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ABSTRACT 


Seismic  modeling  techniques  for  one-  and  two-dimensional  velocity 
models  have  been  developed,  tested  and  applied  to  analysis  of  observed 
seismic  refraction  and  reflection  data  for  the  continental  crust.  The 
reflectivity  method  for  one-dimensional  models  has  proven  to  be  an  efficient 
and  powerful  method  for  interpretation  of  the  amplitude  and  waveform 
of  seismic  record  sections.  The  amplitude  and  waveform  characteristics 
are  shown  to  be  important  to  interpretation  of  fine  structure  of  velocity 
depth  profiles.  Applications  to  modeling  of  a  variety  of  seismic  phases 
are  given.  Two-dimensional  modeling  techniques  utilizing  ray-trace  travel¬ 
time  calculations  and  finite-difference  synthetic  seismogram  calculations 
were  developed.  The  ray- tracing  methods  are  capable  of  accurate  travel - 
time  applications  but  require  modification  for  amplitude  analysis.  The 
finite-difference  method  is  a  powerful  technique  capable  of  modeling 
seismic  data  for  complex  geologic  structures  for  body  and  surface  waves. 
Model  studies  for  simple  one-  and  two-dimensional  velocity  structures 
illustrate  seismic  wave  propagation  including  complex  amplitude  and  waveform 
characteristics  due  to  model  complexity.  The  principal  limitation  of 
the  finite-difference  technique  is  the  large  amount  of  computer  time 
and  storage  required. 
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INTRODUCTION 

During  the  past  few  years  our  research  has  been  aimed  at  developing 
improved  methods  for  modeling  seismograms.  This  effort  has  involved 
both  development  of  new  synthetic  seismogram  techniques  and  implementation 
and  modification  of  existing  computer  programs.  Our  objective  has  been 
not  only  to  provide  adequate  synthetic  seismogram  modeling  techniques, 
which  are  useful  for  a  wide  variety  of  applications,  but  also  to  quanti¬ 
tatively  compare  these  techniques,  ascertain  the  advantages,  disadvantages 
and  range  of  applicability  of  the  various  methods,  and  to  utilize  the 
techniques  to  model  some  real  seismological  data.  Our  applications  have 
largely  been  to  continental  crustal  structure  although  the  techniques 
are  'scale- independent* .  The  methods  can  be  applied  to  a  wide  variety 
of  models  of  geological  and  seismological  interest  including  whole-earth 
seismograms,  crustal  and  upper-mantle  modeling  the  shallow  seismic  refrac¬ 
tion  and  reflection  studies  including  wave  propagation  in  the  ocean  and 
ocean-bottom  region  which  is  of  prime  interest  to  the  Navy.  Both  travel 
time  and  synthetic  seismogram  amplitude  methods  have  been  employed.  Syntheti 
seismogram  methods  utilize  both  ray  theory,  which  is  approximate,  but 
relatively  fast;  and  wave  theoretical  approaches,  which  are  more  exact 
and  include  correct  treatment  of  the  various  phase  conversions  and  wave 
types,  but  may  be  restricted  in  their  range  of  applications  because  of 
the  amount  of  computer  time  required. 

During  the  past  two  years,  our  emphasis  has  been  on  the  development 
of  two-dimensional  synthetic  seismogram  methods.  We  have  utilized  two 
basic  approaches  In  the  development  and  application  of  two-dimensional 
synthetic  seismogram  modeling  techniques.  The  first  involves  a  ray- 
theoretical  method  utilizing  ray  tracing  which  provides  primarily  travel 
times  of  seismic  wave  arrivals.  This  method  can  be  Improved,  however, 
to  produce  reasonably  accurate  synthetic  seismograms  for  seismic  wave 
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propagation  through  complicated  geologic  structures.  Developments  in 
these  areas  will  be  reported  in  a  future  report. 

The  second  method  that  we  have  used  involves  the  numerical  solution 
of  the  two-dimensional  heterogeneous  elastic  wave  equation  using  finite- 
difference  techniques  in  order  to  calculate  complete  synthetic  seismo¬ 
grams  for  complex  two-dimensional  models.  These  techniques  have  the 
advantage  of  providing  complete  seismograms  including  effects  of  wave 
type  conversion,  multiples,  and  inclusion  of  body  waves  as  well  as  surface 
waves.  However,  the  principal  disadvantage  of  the  finite-difference 
techniques  has  been  the  requirement  of  vast  amounts  of  computer  time 
and  storage  in  order  to  apply  these  techniques  on  a  routine  basis. 

In  this  report,  we  summarize  the  status  of  our  synthetic  seismogram 
modeling  development  and  application  efforts.  Previous  reports  relevant 
to  this  research  effort  are  Department  of  Geosciences,  Purdue  University 
reports  0NR-1-80  and  ONR-1-81  prepared  by  Professor  L.W.  Braile.  Addi¬ 
tional  details  of  some  of  the  applications  of  seismic  modeling  performed 
during  the  research  period  are  included  in  reprints  of  published  papers 
which  are  contained  in  the  Appendix  to  this  report.  Recent  developments 
and  improvements  in  finite-difference  modeling  techniques  and  in  an  approxi¬ 
mate  two-dimensional  ray  theoretical  synthetic  seismogram  method,  called 
Disk-Ray  Theory,  will  be  descirbed  in  a  subsequent  report  detailing  research 
peH'ormed  under  contract  number  N00014-82-K-0034. 

SYNTHETIC  SEISMOGRAM  MODELING 

A  summary  of  the  travel  time  and  synthetic  seismogram  modeling  techniques 
that  we  have  utilized  and  have  implemented  on  the  computers  at  Purdue 
University  is  shown  in  Table  1.  This  table  provides  a  summary  of  the 
various  methods  and  the  characteristics  of  these  methods  and  their  associated 
limitations.  For  example,  several  of  the  techniques  are  restricted  to 
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one-dimensional  models  -  that  is,  models  whose  velocities  vary  only  with 
the  depth  parameter.  Other  methods  are  capable  of  calculating  synthetic 
seismograms  for  two-dimensional  models  in  which  lateral  variations  in 
velocity  structure  are  possible. 

Figure  1  is  a  flow  chart  which  illustrates  the  interpretation  procedure 
for  combined  travel-time  and  synthetic  seismogram  modeling  that  we  apply 
to  observed  seismic  data.  A  primary  step  in  this  procedure  is  the  separa¬ 
tion  of  observed  data  into  one-dimensional  and  two-dimensional  classes. 

Once  this  determination  has  been  made,  the  data  are  subject  to  both  travel - 
time  and  synthetic  seismogram  (amplitude  distance)  modeling  procedures 
for  which  the  computer  programs  which  were  described  in  Table  1  are  available. 
In  the  remainder  of  this  report,  we  describe  the  application  of  a  variety 
of  synthetic  seismogram  techniques  which  have  been  developed  or  implemented 
on  our  computers  in  order  to  illustrate  the  capability  of  these  techniques. 
Additional  details  concerning  applications  of  the  synthetic  seismogram 
methods  are  contained  in  the  papers  in  the  Appendix  to  this  report. 

ONE-DIMENSIONAL  SYNTHETIC  SEISMOGRAM  METHODS 

We  have  found  the  Modified  Reflectivity  Method  (Kind,  1978)  to  be 
a  very  powerful  technique  for  synthetic  seismogram  calculation  in  those 
cases  where  the  velocity  structure  of  interest  is  basically  one-dimensional, 
that  is  where  lateral  variations  In  velocity  and  Q  structure  are  not 
significant.  A  variety  of  real  earth  structures  represent  a  close  approxi¬ 
mation  to  this  assumption  and  In  these  cases  the  modified  reflectivity 
method  allows  for  relatively  accurate  and  efficient  modeling  of  observed 
seismic  data.  Some  examples  are  contained  in  the  papers  by  Olsen  and 
Bralle  (1981),  Olsen,  Bralle  and  Johnson  (1980),  Olsen,  Bralle  and  Stewart 
(1982),  and  Banda,  Deichmann,  Bralle  and  Ansorge  (1982)  which  are  given 
in  the  Appendix.  An  additional  example  of  the  application  of  the  synthetic 
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seismogram  technique  to  observe  data  and  the  efficiency  of  calculation 
of  synthetic  seismograms  using  the  modified  reflectivity  method  is  shown 
in  Figure  2.  Our  experience  with  the  modified  reflectivity  method  indicate 
the  importance  of  amplitude  and  waveform  modeling  of  observed  seismic 
data.  Investigation  of  travel-time  relationships  alone  provides  little 
information  as  to  the  detailed  earth  model.  Inclusion  of  amplitude  distance 
variations  and  waveform  characteristics,  which  can  be  treated  using  the 
synthetic  seismogram  technique,  provide  for  much  improved  inference  of 
the  fine  structure  of  the  velocity  depth  curve  as  well  as  inference  of 
Q  structure  (Braile,  1978). 

TWO-DIMENSIONAL  MODELING  -  RAY  THEORY 

Ray- theoretical  techniques  have  been  used  for  modeling  of  travel 
times  of  seismic  waves  in  complex  geologic  structures  for  some  time. 

We  have  implemented  an  efficient  ray-tracing  program  for  accurate  calcula¬ 
tion  of  travel  times  of  refracted  and  reflected  seismic  waves  following 
the  method  of  Cerveny,  Molotkov  and  Psencik  (1978).  This  technique  utilizes 
a  variety  of  numerical  schemes  to  approximate  a  two-dimensional  velocity 
distribution  with  or  without  interfaces.  The  seismic  ray  paths  are  then 
traced  using  iterative  Snell's  Law  application  through  the  velocity  structure. 
Reflections  from  various  Interfaces  may  be  selected  and  refracted  waves 
are  also  Included.  Refractions  simulating  head  waves  from  homogeneous 
layers  can  be  adequately  approximated  by  providing  for  a  small  positive 
velocity  gradient  within  the  homogeneous  media.  The  velocity  gradient 
that  we  normally  use  to  simulate  the  head  wave  is  consistent  with  that 
utilized  In  the  earth  flattening  transformation  and  corresponds  to  approxi¬ 
mately  0.001  km/s/km  for  shallow  layers.  This  technique  provides  adequate 
modeling  of  body-wave  travel  times  in  two-dimensional  structures.  An 
example  of  successful  modeling  of  complex  observed  seismograms  using 
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two-dimensional  ray- trace  modeling  is  shown  in  Figure  3.  One  of  the 
limitations  of  this  technique  is  that  amplitudes  calculated  for  individual 
ray  paths  are  inaccurate  and  thus  the  ray  tracing  technique  is  not  inmediately 
conducive  to  synthetic  seismogram  calculation.  However,  we  have  utilized 
a  modified  ray  theoretical  method  following  the  technique  of  Wiggins 
(1976)  and  McMechan  (1974)  to  produce  accurate  synthetic  seismograms 
from  the  ray  trace  amplitudes.  Progress  in  this  area  will  be  described 
in  a  subsequent  report. 

FINITE-DIFFERENCE  SYNTHETIC  SEISMOGRAM  MODELING 

We  have  applied  the  finite-difference  method  to  numerically  solve 
the  two-dimensional  heterogeneous  elastic  wave  equation  for  the  calcula¬ 
tion  of  synthetic  seismograms  in  velocity  structures  with  arbitrarily 
varying  elastic  properties.  The  idea  of  the  finite-difference  synthetic 
seismogram  method  is  illustrated  in  Figure  4.  A  velocity  structure  is 
approximated  by  a  rectangular  grid  of  points  in  which  the  elastic  properties 
( compress ional  velocity,  shear  velocity,  and  density)  are  tabulated. 
Finite-difference  approximations  to  spatial  and  temporal  derivatives 
allow  solution  of  the  two-dimensional  elastic  wave  equation  for  a  given 
instant  in  time  as  computed  from  displacements  at  previous  times.  Initial 
conditions  due  to  an  implied  elastic  disturbance  (source)  generates  wave 
propagation  which  is  simulated  by  the  finite-difference  calculations 
as  a  function  of  time  for  all  points  in  the  velocity  grid.  Displacement 
time  histories  at  a  variety  of  locations  are  retained  in  the  computer 
and  become  displacement  seismograms  for  the  vertical  and  horizontal  components 
of  particle  velocity.  The  theory  of  finite-difference  synthetic  seismogram 
calculations  has  been  described  by  Boore  (1972),  Alford  et  a^.,  1974; 

Kelly  eta^.,  1976;  Mazzella,  1979;  and  Espindola,  1979. 
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The  two-dimensional  equations  of  motion  for  displacement  in  a  hetero¬ 
geneous  isotropic  elastic  media  are  shown  in  Figure  5.  These  equations 
are  approximated  by  their  finite-difference  approximations  using  an  explicit 
finite-difference  formulation  as  illustrated  in  Figure  6.  Using  these 
finite-difference  approximations,  the  displacement  at  any  time  L  +  1 
for  any  arbitrary  grid  location  is  given  as  a  function  of  times  L  and 
L  -  1  allowing  a  calculation  of  the  displacement  for  all  grid  locations 
at  the  new  time.  Once  this  calculation  is  complete  for  the  entire  veloc 
grid,  it  may  be  repeated  for  as  many  time  steps  as  desired. 

This  finite-difference  application  is  very  straightforward  and  resu. 
in  complete  seismograms  since  it  is  a  numerical  solution  of  the  elastic 
wave  equation.  However,  there  are  a  variety  of  difficulties  which  one 
encounters  in  attempting  to  apply  finite-difference  calculations  to  the 
calculation  of  synthetic  seismograms.  These  difficulties  are  summarized 
in  Figure  7.  The  most  difficult  problems  are  that  of  providing  for  absorbing 
boundary  conditions  at  the  edges  of  the  model  in  which  an  approximate 
absorbing  boundary  condition  due  to  Clayton  and  Engquist  (1977)  is  utilized 
and  the  requirement  for  stability  and  accuracy  which  imply  small  time 
step  and  small  grid  spacing  respectively.  These  requirements  result 
in  very  large  computer  time  and  storage  capacity.  Therefore,  calculation 
of  synthetic  seismograms  for  complicated  realistic  models  of  interest 
for  two-dimensional  geological  structures  may  require  several  tens  of 
minutes  or  even  hours  of  computer  time  utilizing  a  computer  with  a  memory 
capacity  of  a  million  or  more  storage  locations, 

A  finite-difference  synthetic  seismogram  computer  code  has  been 
developed  and  is  described  in  Mazella  (1979).  We  have  applied  this  synthetic 
seismogram  method  to  a  variety  of  models  in  order  to  test  the  accuracy, 
applicability,  and  capability  of  the  technique.  The  velocity  models 
utilized  for  this  testing  are  illustrated  in  Figure  8.  One  of  the  useful 
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by-products  of  the  synthetic  seismogram  calculation  using  the  finite- 
difference  technique  is  a  displacement  field  associated  with  wave  propagation 
in  the  elastic  model  at  any  instant  in  time.  These  displacement  fields 
can  be  stored  for  a  number  of  time  steps  during  the  finite-difference 
calculation  and  the  displacement  fields,  often  called  'snapshots',  are 
illustrative  of  the  wave  propagation  in  the  media.  An  example  of  dis¬ 
placement  time  history  snapshots  for  one  of  the  models  shown  in  Figure 
8  are  illustrated  in  Figure  9. 

In  order  to  verify  the  accuracy  of  the  finite-difference  program, 
we  have  calculated  vertical  component  synthetic  seismograms  for  a  one¬ 
dimensional  model  corresponding  to  a  layer  over  a  half-space  for  which 
synthetic  seismograms  have  also  been  computed  utilizing  the  modified 
reflectivity  method.  Comparison  of  the  seismograms  for  these  two  calcula¬ 
tions  (Figure  10)  illustrates  the  validity  of  the  finite-difference  calculations. 

Calculation  of  synthetics  for  a  complex  two-dimensional  velocity 
structure  (model  INFL  5  illustrated  in  Figure  8)  are  shown  in  Figure  11. 

The  corrected  seismograms  (Figure  11 B)  for  this  complicated  velocity 
structure  illustrate  the  expected  waveform  characteristics  of  propagation 
through  this  model  containing  a  prominent  lateral  velocity  contrast  in 
the  form  of  a  fault.  The  time  delay  due  to  the  fau’t  and  diffracted 
and  headwave  arrivals  from  the  lower  interface  demonstrate  the  capability 
of  this  technique  for  simulating  wave  propagation  and  complex  geologic 
structures. 

The  principal  difficulty  with  the  finite-difference  approach  is 
the  amount  of  computer  time  required  to  compute  synthetic  seismograms 
for  realistic  velocity  structures.  It  is  currently  too  expensive  to 
compute  synthetics  using  this  technique  on  a  routine  basis  for  trial 
and  error  model ing  of  seismic  refraction  and  reflection  data.  However, 
it  is  feasible  to  compute  synthetics  for  a  variety  of  characteristic 
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models  of  laterally  varying  geologic  structures  and  to  use  the  experience 
gained  with  these  models  as  an  aid  to  interpretation  of  observed  seismic 
data.  In  addition,  it  is  highly  desirable  to  have  two-dimensional  synthetic 
seismograms  calculated  by  the  finite-difference  technique  for  certain 
models  in  order  to  use  for  comparison  with  approximate  techniques  as 
verification.  Furthermore,  it  is  possible  that  the  finite-difference 
techniques  could  be  made  to  be  somewhat  more  efficient,  thus  improving 
our  ability  to  utilize  the  finite-difference  synthetic  seismogram  technique 
on  a  routine  basis.  For  example,  we  are  currently  investigating  the 
possibilities  of  using  an  acoustic  formulation  of  the  two-dimensional 
finite-difference  procedure  as  a  preliminary  modeling  method.  In  addition, 
we  are  developing  implicit  finite-difference  schemes  for  synthetic  seismo¬ 
gram  modeling  which  may  be  significantly  faster  than  the  explicit  code 
that  we  are  presently  utilizing. 
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Figure  1 . 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


FIGURE  CAPTIONS 


Flow  chart  illustrating  the  interpretation  procedure  for 
combined  travel- time  and  amplitude  analysis  for  seismic 
refraction  and  reflection  record  sections  for  one-  and  two- 
dimensional  velocity  structures. 

Observed  and  synthetic  seismic  record  sections  for  a  seismic 
refraction  profile  recorded  on  the  eastern  Snake  River  Plain 
(Braile  et  art.,  1982).  Both  record  sections  have  been  ampli¬ 
tude  scaled  with  amplitudes  being  multiplied  by  distance 
to  the  1.5  power  in  order  to  enlarge  the  amplitudes  at  large 
distance  to  make  arrivals  visible.  Travel-time  curves  on 
both  seismograms  are  for  a  plane  layered  velocity  model  which 
was  used  to  calculate  the  synthetic  seismograms  illustrated 
in  the  lower  record  section.  Good  comparison  is  seen  between 
the  travel-time,  amplitude  and  waveform  characteristics  of 
the  primary  phases  for  the  observed  and  synthetic  sections 
indicating  that  the  velocity  and  Q  model  used  for  the  synthetic 
seismogram  calculation  is  substantially  correct. 

Ray- trace  diagram  and  travel  times  compared  with  observed 
seismic  refraction  data  across  the  eastern  Snake  River  Plain 
(from  Sparlin  et  a^. ,  1982).  The  ray  diagrams  show  computer 
plotted  ray  traces  for  refracted  (upper  ray  trace  diagram) 
and  reflected  (lower  ray  trace  diagram)  travel  times  through 
the  complex  velocity  structure.  The  numbers  on  the  velocity 
structure  indicate  compressional  wave  velocity  in  km/s.  The 
travel  times  for  the  various  phases  (designated  A,  B,  C, 
etc.)  are  indicated  on  the  seismic  record  section  at  the 
top  of  the  diagram.  The  observed  travel  times  show  a  good 
match  to  the  theoretical  travel  times  computed  by  the  ray¬ 
tracing  technique. 

Schematic  diagram  illustrating  the  velocity  model  for  finite- 
difference  synthetic  seismogram  calculations.  The  velocity 
structure  in  the  seismic  model  can  vary  arbitrarily  and  is 
specified  by  elastic  properties  at  grid  points  which  are 
distributed  throughout  the  model.  The  source  represents 
initial  conditions  which  are  used  to  calculate  displacements 
as  a  function  of  time  throughout  the  velocity  model.  The 
source  can  be  located  at  any  position  within  the  velocity 
structure.  A  free-structure  boundary  condition  is  applied 
in  the  finite-difference  calculation  and  absorbing  boundary 
conditions  are  utilized  on  the  edges  of  the  velocity  structure 
to  minimize  effects  of  the  fictitious  boundaries.  The  seismo¬ 
meter  locations  can  be  arbitrarily  selected,  but  are  usually 
equally  spaced  along  the  surface  of  the  model  and  represent 
locations  where  the  displacement  time  histories  will  be  stored 
in  the  computer  memory  to  be  plotted  as  seismograms  at  the 
end  of  the  finite-difference  calculation. 

Two-dimensional  equations  of  motion  for  displacement  for 
a  heterogeneous  isotropic  elastic  model. 
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Figure  6.  Form  of  the  finite-difference  approximations  for  second  order 
spatial  derivatives  and  cross  product  spatial  derivatives 
and  second  order  time  derivatives  which  are  approximations 
to  the  two-dimensional  elastic  equation  of  motion  for  displace¬ 
ments  illustrated  in  Figure  5. 

Figure  7.  List  of  factors  which  must  be  considered  in  explicit  finite- 
difference  synthetic  seismogram  calculations. 

Figure  8.  Velocity  depth  models  used  for  finite-difference  synthetic 
seismogram  calculations.  Distance  (X)  and  depth  (Z)  are 
given  for  each  model.  The  lower  boundary  of  each  model  is 
assumed  to  be  an  infinite  half-space,  but  in  fact,  in  the 
program  it  consists  of  a  non- reflecting  artificial  boundary. 

The  numbers  within  the  model  indicate  the  compressional  wave 
velocity  in  km/s  for  each  layer  of  the  models.  Synthetic 
seismogram  record  sections  and  snapshots  of  displacement 
time  histories  are  shown  for  the  velocity  models  indicated 
in  this  figure. 

Figure  9.  Displacement  time  history  snapshots  for  model  INFL2.  The 

Z  component  is  shown  for  times  1,  2,  3,  4,  5,  6  and  7  seconds 
propagation  time.  The  source  is  located  at  i  km  depth  and 
2  km  distance  from  the  origin.  Distance  (X)  and  depth  (Z) 
directions  are  shown.  The  perspective  diagram  indicates 
the  Z  component  of  displacement  at  each  instant  of  time  for 
the  7  different  times  represented  by  the  snapshots.  This 
velocity  model  consists  of  a  single  layer  over  a  half-space. 

The  source  consists  of  a  compressional  point  explosive  source. 

As  the  seismic  waves  propagate  through  the  model,  the  dis¬ 
placement  in  X  and  Z  directions  can  be  visualized  by  these 
snapshots.  Note  that  the  moving-window  process  utilized 
in  the  synthetic  seismogram  program  causes  the  area  of  the 
model  to  be  shifted  to  the  right  (in  the  direction  of  increasing 
X)  for  successive  time  steps.  This  movement  of  the  window 
can  be  seen  by  the  shift  in  the  X  axis  values  shown  at  the 
bottom  of  each  snapshot.  The  efficiency  of  the  non-reflecting 
boundary  can  be  seen  by  the  snapshots,  particularly  at  time 
step  T=2  seconds  in  which  a  compressional  wave  is  impinging 
upon  the  lower  boundary  of  the  window  space  and  no  prominent 
reflection  can  be  seen  from  this  non-reflecting  boundary. 

The  development  of  a  head  wave  due  to  propagation  in  the 
6  km/s  half-space  below  the  4.5  km/s  sedimentary  layer  can 
be  seen  beginning  at  approximately  T=3  seconds  and  becoming 
more  pronounced  with  increasing  time.  The  head  wave  is  clearly 
seen  in  the  T=6  snapshot  in  which  a  strong  compressional 
pulse  is  propagating  at  a  sharp  angle  upwards  and  to  the 
right  through  the  sedimentary  layer  toward  the  surface. 

Figure  10.  Synthetic  seismograms  calculated  for  model  INFL3  (Figure  8) 

which  consisted  of  a  layer  over  a  half-space.  The  Z  component 
seismograms  were  sampled  at  a  5  km  distance  interval  and 
are  plotted  as  a  reduced-time  record  section.  The  strongest 
arrivals  are  the  Rayleigh  waves  which  begin  at  about  1$  seconds 
reduced  time  on  the  5  km  distance  seismogram  but  propagate 
at  a  relatively  slow  apparent  velocity.  The  Rayleigh  waves 
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Figure  11 


are  left  behind  by  the  moving-window  process  and  thus  are 
not  prominent  after  the  first  few  seismograms.  The  direct 
wave  in  the  surface  layer  is  visible  as  the  4.5  km/s  apparent 
velocity  arrival  indicated  on  the  record  section.  The  head 
wave  from  the  6  km/s  half-space  is  clearly  seen  as  a  first 
arrival  from  about  15  km  to  the  maximum  50  km  distance  range 
considered.  Additional  arrivals  consist  of  P-S  conversions 
and  multiple  refractions  and  reflections.  Some  numerical 
noise  exists  in  the  later  sections  of  several  of  the  seismo¬ 
grams.  However,  this  noise  tends  to  be  truncated  by  the 
moving  window  process  as  the  window  proceeds  ahead  of  the 
slowly  propagating  noise  and  thus  the  numerical  noise  does 
not  interfer  with  the  primary  phases  which  propagate  in  the 
early  portion  of  the  record  section.  Also  shown  for  comparison 
is  a  record  section  calculated  for  the  exact  same  model  (INFL  3) 
by  the  modified  reflectivity  method  for  a  similar  source 
wavelet.  The  reflectivity  calculation  included  phase  velo¬ 
cities  of  2.5  -  100  km/s  and  thus  ignores  some  of  the  shear 
and  surface  wave  arrivals.  The  correspondence  between  the 
two  record  sections  is  quite  good  demonstrating  the  accuracy 
of  the  finite-difference  method. 

.  Synthetic  seismogram  record  sections  for  finite-difference 
calculation  for  model  INFL  5  (Figure  8).  This  model  consists 
of  a  two-dimensional  velocity  model  with  a  prominent  fault 
of  5  km  offset  which  occurs  15  km  from  the  source.  The  seismo¬ 
grams  shown  in  the  record  section  indicate  the  beginnings 
of  a  6  km/s  apparent  velocity  head-wave  arrival  in  the  first 
three  seismograms  and  the  delay  due  to  the  fault  structure 
at  15  km  distance  is  seen  beginning  at  approximately  20  km 
on  the  record  section.  However,  full  visualization  of  the 
two-dimensional  effect  of  the  fault  structure  on  the  synthetic 
seismogram  is  prevented  by  the  presence  of  some  peculiar 
numerical  noise  which  occurs  in  the  record  section.  This 
noise  is  caused  by  an  improper  choice  of  the  velocity  of 
the  moving-window  which  is  used  to  speed  up  the  calculations 
of  the  finite-difference  synthetic  seismogram  program.  In 
the  moving-window  calculation,  the  finite-difference  equations 
are  solved  for  only  a  portion  of  the  velocity  model,  called 
a  window,  and  this  window  is  caused  to  move  through  the  velo¬ 
city  model  in  such  a  way  as  to  constantly  include  the  phases 
of  Interest.  The  window  is  shifted  in  discrete  steps  to 
approximately  move  at  the  average  horizontal  velocity  of 
the  seismic  waves.  However,  in  this  particular  example  the 
velocity  of  the  window  was  selected  slightly  too  slow  and 
the  wave  propagation  of  the  first  arrival  (compress ional 
waves)  Impinged  upon  the  leading  edge  of  the  moving  window. 
During  the  next  shift  in  the  window,  this  resulted  in  trunca¬ 
tion  of  the  displacements  at  the  former  window  position  and 
thus  this  truncation  effect  caused  the  numerical  noise  which 
is  seen  in  the  seismograms  from  25  to  45  km  distance  range. 

This  effect  can  certainly  be  avoided  by  rerunning  the  synthetic 
seismogram  calculations  utilizing  a  slightly  faster  velocity 
for  the  moving -window.  Vertical  (Figure  11 B)  and  radial 
(Figure  11 C)  component  seismograms  calculated  by  the  finite- 
difference  method  for  model  INFL  5  (Figure  8)  with  a  corrected 
window  velocity.  The  6.0  km/s  head  wave  is  seen  to  be  delayed 


significantly  and  the  amplitudes  affected  by  the  fault  which 
occurs  at  a  distance  of  15  km  In  the  model.  Much  of  the 
energy  on  the  record  sections  is  due  to  P-S  conversions  and 
Rayleigh  waves  which  are  also  Included  in  the  calculation. 

Note  the  small  numerical  dispersion  of  the  Rayleigh  wave 
train  at  about  3  to  4  seconds  reduced  time  on  the  5  km  seismo¬ 
gram  on  the  vertical  component  record  section.  This  dispersion 
is  due  to  too  large  of  a  grid  spacing  in  the  model  for  the 
slowly  propagating  Rayleigh  waves.  The  spacing  is  however 
adequate  for  the  body  waves  and  thus  only  affects  the  Rayleigh 
wave  which  in  this  case  is  not  of  primary  interest. 


TABLE  1.  SUMMARY  OF  SEISMIC  INTERPRETATION  TECHNIQUES 


METHOD 

CHARACTERISTICS 

LIMITATIONS 

Geometrical  Ray 
theory  calculation 
of  travel-times  of 
reflected  and  refracted 
phases  (Program 

TXCURV) . 

Travel -times  for 
head  waves  and 
reflections  for 
plane  or  dipping 
layered  models; 
exact . 

Model  geometry; 
travel-times  only. 

Ray  Tracing 
(Program  RAY2D) 

Two-dimensional 

Models 

Travel-times  only, 
limited  number  of 
phases ,  ignores  P-S 
conversions ,  guided 
waves  and  surface 

waves. 

Reflectivity 
(Program  SYNCAL) 

Wave-theoretical 
solution  to 
response  of  a  lay¬ 
ered  half  space. 
Includes  P,  SV, 
guided  phases  and 
surface  waves. 

Exact.  Includes  Q. 

Restricted  to  1-D 
models. 

Disk  Ray  Theory 
(modified  RAY 2D 
program) 

2-D  models,  approxi¬ 
mate  solution  for  p 
waves,  can  include 
reflection  coeffi¬ 
cients  and  Q.  Com¬ 
putationally  effi¬ 
cient. 

Inexact  amplitudes  for 
head  waves  and  near- 
critical  arrivals;  ig¬ 
nores  guided  phases, 

P-S  conversions  and 
surface  waves. 

Finite  Difference 

Acoustic  (Program 

FDWVAC) 

2-D  models,  exact 
solution  for  all 
wave  types  in  an 
acoustic  media. 
Significantly  faster 
than  FDWVEQ,  but 
slower  than  RAY 2D 
synthetics. 

Acoustic  rather  than 
elastic .  Large  com¬ 
puter  time. 

Finite  Difference 

Elastic  (Program 

FDWVEQ) 

2-D  models,  exact 
solution  for  all 
wave  types  in  elastic 
media.  Can  be  modi¬ 
fied  to  include  Q. 

Extensive  computer 
time  and  storage . 
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INTERPRETATION  PROCEDURE  FOR  COMBINED 
TRAVEL  -  TIME  AND  AMPLITUDE  ANALYSIS  OF  SEISMIC 
REFRACTION  RECORD  SECTIONS 


Figure  1. 


TWO  DIMENSIONAL  EQUATIONS  OF  MOTION 
FOR  DISPLACEMENT 
( HETEROGENEOUS,  ISOTROPIC  MEDIA  ) 
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IN  x  AND  2  DIRECTIONS 


a 


X  +  Zp. 
P 


IS  THE  COMPRESSIONAL 
VELOCITY 


0 


IS  THE  SHEAR  VELOCITY 


P 


IS  THE  DENSITY 


Figure  5. 
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FINITE  DIFFERENCE  APPROXIMATIONS 
C  EXPLICIT  CASE  ) 


SECOND  ORDER  SPATIAL  DERIVATIVES 


•gjj-  j«2  (x,  z)g^-  u(x,  z,  t)j  =  |  oz(m+-j,  n)  ^u(m+  I ,  n,  I)  -  u(m,  n,  l)J 
-a2(m-j,  (0  ^u(m,  n,  D-  u(m-1,  n,  l)j  y/(&H)z 


CROSS  PRODUCT  SPATIAL  DERIVATIVES 


g^-£a2(x,  z)  gj|-uCx,  z,  t)j  =  |  a2(m,  n+  1 )  ^u(ni+ 1 ,  n+1 ,  I)  -  u(m- 1 ,  n+1 ,  l)j 
-  «2<m,  n- 1 )  ^u(m+ 1 ,  n- 1 , 1)  -  u(m- t ,  n- 1 ,  oj  y^CAK)2 


SECOND  ORDER  TIME  DERIVATIVE 
u(w,  n.  1+1)  *  2u(m,  n,  I)  -  u(m,  n,  1-1)  +  f  ^Spatial  derivatives 

WHERE:  az<m±y,  n)  =  [a2(m  +  1,  n)  +  a2(m,n)j  j 2 

m  AND  n  REFER  TO  THE  x  AND  z  DIRECTION 
GRID  POINTS  AND  I  IS  THE  TIME  STEP 


Figure  6. 
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PRACTICAL  ASPECTS  OF  EXPLICIT  FINITE  DIFFERENCE 
SYNTHETIC  SEISMOGRAM  CALCULATION 

-  TWO  DIMENSIONAL  APPLICATION 

-  SOURCE  GENERATION 

-  FREE  SURFACE  BOUNDARY  CONDITION 

-  ABSORBING  BOUNDARIES  AT  EDGES  OF  MODEL 

-  STABILITY  CONDITION  (SMALL  TIME  STEP) 

-  ACCURACY  CONDITION  (SMALL  GRID  SPACING) 

-  NUMBER  OF  TIME  STEPS  (LENGTH  OF  SEISMOGRAMS) 

-  NUMBER  OF  GRID  POINTS 

-  LARGE  COMPUTER  TIME  AND  STORAGE  REQUIREMENTS 


Figure  7. 
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Figure  9A. 


Figure  9B. 


Figure  9E. 
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Figure  9F. 
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INFL5 

2-COMPONENT 


Figure  HA. 


FINITE  DIFFERENCE  SYNTHETIC  SEISMOGRAMS 


Figure  11B. 
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Abstract.  Average  up  par  nantla  Pn  velocities 
and  hast  flow  wara  computed  within  contlnantal 
physiographic  provlncaa  In  North  Aaarlca  fro* 
publish ad  data.  F#  velocity  and  haat  flow  data 
display  an  lnvaraa  ra la t Ions hip  and  wara  found 
to  corralata  with  a  statistically  significant 
•  corrals tlon  coaff lclant.  Tsaparaturas  at  tbs 
rrnar  want  Is  boundary  wara  estimated  frow  tha 
haat  flow  valuas,  and  thasa  wara  uaad  to  damon- 
strata  a  corralatlon  batwaan  Fa  velocity  and 
temperature  of  upper -man tie  notarial.  Tha 
valus  of  (3Fn/3T)p  thus  obtained  (-4.4  x  10~* 
to  -8.1  x  10"*  ka  s-1  *C_l)  la  within  tha  range 
.of  taaperature  derivatives  determined  from 
laboratory  studies  of  ultraaaf lc  rocks.  This 
dependency  of  Pn  velocity  on  taaperature  im¬ 
plies  that  one  possible  explanation  for  tha 
observed  geographical  distribution  in  upper- 
nan  tie  aelsalc  velocity  la  that  tha  variation 
in  PQ  velocity  la  prlaarlly  a  taaperature 
affect.  Contained  with  tha  relationship  ba¬ 
twaan  haat  flow  and  crustal  age  for  continents, 
the  Pn  versus  haat  flow  correlation  suggests 
a  relationship  between  Pn  velocity  and  crustal 
age,  probably  due  to  progressiva  cooling  of 
tha  contlnantal  lithosphere  after  a  thermo- 
tec tonic  event. 

Introduction 

Seismologists  have  bean  aware  of  variations 
of  contlnantal  upper-mantle  compress lone 1-wave 
velocity  (Pn,  velocity  at  tha  Mohorovi2l£  dis¬ 
continuity)  on  a  regional  scale  since  about 
1960  (Herrin,  1969].  Tbs  Implications  of  this 
observation  are  not  dearly  understood,  however, 
and  several  possible  explanations  have  bean 
suggested.  Horal  and  Slnaona  [1968]  found  a 
correlation  between  travd-tlae  anomalies  for 
tha  Longshot  nuclear  explosion  and  heat  flow. 
They  could  not  demonstrate  a  similar  corres¬ 
pondence  with  gravity.  They  concluded  that 
thermal  anomalies  in  the  upper  mantle  (perhaps 
30- ka  deep)  were  the  probable  cause  of  travel- 
time  anomalies.  Pakiser  and  Stelnhart  [1964], 
and  Wgrrcn  and  Haaly  [1973],  on  the  other  hand, 
related  Pa  and  mean  crustal  vdocitles  to  den¬ 
sity  variations.  Chung  [1977]  attributed  varia¬ 
tions  In  P  to  compositional  differences  In 
the  upper  mantle.  Fagernes  and  Xanestrom  [1973] 
found  that  variations  In  the  ratio  of  Pa  velo¬ 
city  to  S-  velocity  within  a  region  were  too 
large  to  be  explained  solely  by  temperature 
variations,  and  proposed  density  variations 
as  the  cause.  Pakiaer  [1963],  noting  that  Pn 
and  cruatal  thicknaas  are  related,  found  that 
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the  assumption  of  iaoetasy  requires  variations 
of  density  In  the  upper  mantle  and  hehca  a 
relationship  between  Fa  and  density.  Ha  notes, 
however,  that  '...elevated  tmeperstures  nay 
affect  velocity  more  than  density'  [Pakiser. 

1963,  p.  3734]. 

This  study  was  undertaken  In  an  attempt  to 
quantitatively  determine  the  relationship  be¬ 
tween  continental  Pa  velocities  and  heat  flow, 
and  to  determine  whether  the  <«pHeH  tempera¬ 
ture  variations  alone  are  sufficient  to  be  a 
viable  explanation  for  the  variation  In  upper 
mantle  seismic  velocity.  Tha  continental  crust 
and  upper  mantle  of  North  Aaarlca  was  selected 
for  study  because  of  the  large  n usher  of  seis¬ 
mic  refraction  profiles  and  heat  flow  measure¬ 
ments  available  and  because  of  tbs  range  of 
haat  flow  valuea  present.  Because  the  loca¬ 
tions  of  seismic  and  thermal  observations  do 
not  coincide,  and  because  of  large  uncertainties 
which  are  possible  in  any  given  observation, 
the  data  were  grouped  by  provinces  and  mean 
values  for  each  province  analysed. 

SelMlc  Data 

One-hundred  thirty-one  seismic  refraction  pro¬ 
files  for  the  continental  crust  and  upper  mantle 
of  North  America  were  compiled  from- published  lit¬ 
erature  and  other  sources  (Table  1).  Observed 
Pa  velocities  and  crustal  thickness  were  tabu¬ 
lated  from  the  profiles  resulting  In  133  observa¬ 
tions  of  tha  continental  Pn  velocity.  The  data 
were  sorted  according  to  physiographic  province 
(Fennenan  [1946]  for  the  United  States  and 
Douglas  and  Pries  [1972]  for  Canada)  and  naan 
Pa  velocity  and  crustal  thickness  were  calculated 
for  each  province  (Table  2).  The  locations 
of  the  seismic  profiles  and  the  physiographic 
provinces  are  shown  in  Figure  1.  Estimated 
errors  for  the  mean  Pn  velocity  and  crustal 
thickness  were  taken  to  be  the  estimated  stan¬ 
dard  deviation  of  the  mean,  S  *  a/&.  where 
o  la  the  standard  deviation  of  the  data  (Pn 
or  for  a  given  province)  and  N  is  the  number 
of  observations.  Tha  estinated  standard  devia¬ 
tion  of  the  mean  provides  a  reasonable  choice 
for  estimating  the  uncertainty  of  the  mean  P„ 
velocity  for  each  province  because  possible 
errors  of  individual  P  observations  are  gener¬ 
ally  not  given.  However ,  S  la  likely  to  be 
a  poor  estimate  of  the  uncertainty  of  the  mean 
Pa  value  for  provinces  containing  only  a  small 
number  of  observations  such  as  tha  Cascade 
Range,  middle  Rocky  Mountains,  southern  Rocky 
Mountains,  and  Columbia  Plateau  provinces. 

Two  possible  sources  of  errors  In  the  P„ 
observations  have  not  been  analysed  In  detail, 
but  are  expected  to  be  negligible  for  tha  pur¬ 
pose  of  correlation  of  Pa  velocity  and  heat 
flow  by  province.  They  are  velocity  anisotropy 
and  dlffsrsncss  in  the  apparent  and  true  Pn  velo¬ 
city  due  to  the  earth's  sphericity.  Although 
velocity  anlsotrophy  In  tha  upper  mantle  of 
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TABU  1.  Seismic  Refraction  Profiles  for  North  America 


Profile 

Beference 

Profile 

Beference 

1 

Richards  and  Welker  [1959] 

178 

Hales  and  Ration  [1973] 

2.3 

Johnson  and  Couch  [1970] 

179 

Hill  [1972] 

4,3,6, 7 

Ewing  et  al.  [1966] 

184 

Jackson  and  Fakiser  [1965] 

6.9 

Stewart  [1968a] 

185 

Brails  et  al.  [1974] 

10 

Boiler  [1965] 

193 

Keller  et  al.  [1975] 

11 

Dlaent  et  al.  [1961] 

194 

Hales  and  Ration  [1973] 

U 

Boiler  [1964];  Prodehl  [1979] 

195 

Olaen  at  al.  [1979] 

11 

Langaton  and  Belabor gar  [1974] 

198 

Barr  [1971] 

12 

Hamilton  at  al.  [1964] 

199 

Mereu  and  Hunter  [1969] 

13 

%erg  et  al.  [1966] 

200 

Hodgson  [1953] 

14 

Johnson  [1965] 

201,202 

Hall  and  Hajnal  [1973] 

13,16 

Hit.  hall  and  Landlaaan  [1971] 

204,205 

Green  at  al.  [1980] 

17 

Topposada  and  Sanford  [1976] 

206 

Lyons  et  al.  [1980] 

18.19 

Haaly  [1963] 

207 

White  and  Savage  [1965] 

20 

Boiler  and  Jackson  [1966] 

208 

Warren  et  al.  [1972] 

21,22,23 

Bats  [1954] 

212 

Hill  and  Fakiser  [1966] 

28,29 

Baton  [1963] 

216 

Carder  [1973] 

31 

Gibbs  and  Boiler  [1966] 

217 

Forsyth  at  al.  [1974] 

32 

Carder  et  al.  [1970] 

218,219 

Hobson  at  al.  [1967] 

33 

Ryall  and  Stuart  [1963] 

220 

Giah  at  al.  [1981] 

33 

Berry  and  West  [1966] 

221 

Saieh  ae  al.  [1982] 

36,37 

Cohan  and  Meyer  [1966] 

222 

Sparlln  et  al.  [1982] 

38 

Craa  [1961] 

223 

Baldwin  [1980] 

39 

Doraan  et  al.  [1972] 

225 

Hall  and  Hajnal  [1969] 

40 

MeCany  and  Meyer  [1966] 

226 

Slnno  et  al.  [1981] 

41 

Merkel  end  Alexander  [1969] 

227,228 

Shor  at  al.  [1968] 

42 

Warren  et  al.  [1966] 

229 

Keller  et  al.  [1975]; 

43,44 

Warren  [1969] 

Mueller  and  Landlaaan  [1971] 

45 

Hales  et  al.  [1970] 

230 

Clee  at  al.  [1974] 

46 

Willden  [1963] 

231 

Martin  [1978] 

47 

Jackson  et  al.  [1963] 

232 

Mereu  and  Jobldon  [1971] 

48 

Jackson  and  Fakiser  [1965] ; 

233,234, 

Prodehl  and  Fakiser  [1980] 

235.236 

Berry  and  Forsyth  [1975] 

49 

Chandra  and  Ci—lnga  [1972] 

237,238 

Eaton  [1966] 

50 

Bennett  et  al.  [1975] 

239.240 

Stewart  [1968b] 

51,32 

Johnson  at  al.  [1972] 

241 

Brails  at  al.  [1982] 

53,54.55 

Berry  and  Fuchs  [1973] 

242.243, 

56.57.58. 

244.245, 

59,60,61, 

246.247, 

62,63 

MeCany  and  Meyer  [1964] 

248 

Prodehl  [1979] 

64 

Boiler  and  Haaly  [1963] 

249,250. 

66 

Ewing  et  al.  [1955] 

251 

Warren  and  Jackson  [1968] 

68.69 

Herasy  et  al.  [1959] 

252.253 

Bonney  at  al.  [1962] 

72 

Mereu  et  al.  [1976] 

254.255 

Janas  at  al.  [1968] 

75.76 

Stauber  end  Boors  [1978] 

256 

Bates  and  Ball  [1975] 

86,87 

Barrett  et  al.  [1964] 

257 

Tuva  [1951];  Stelnhart  and  Mayer  [1961] 

154 

Shor  [1962] 

238 

Stelnhart  and  Meyer  [1961] 

137 

Stelnhart  and  Meyer  [1961] 

259 

Tuve  [1953];  Stelnhart  and  Mayer  [1961] 

158 

Slichtsr  (1951; 

260 

Tuve  [1954];  Stelnhart  and  Meyer  [1961] 

Stelnhart  and  Meyar  [1961] 

261 

Warren  [1968] 

163 

Stelnhart  and  Meyer  [1961] 

262 

Stelnhart  et  al.  [1964] 

167 

Press  [I960] 

263 

Hales  et  al.  [1968] 

173,174 

Slichtsr  [1951]; 

Stelnhart  and  Meyer  [1961] 

175 

Berg  et  al.  [1960] 

continent*  ha*  boon  reported  [Sanford,  1973], 
an  analyele  of  western  United  Statee  data  by 
Banford  at  al.  [1979]  found  a  relatively  snail 
and  poorly  doflnad  anlaotropy  of  about  3Z.  Ba¬ 
ca  vis*  the  azimuths  of  the  aelaalc  refraction 
profile*  wi thin  each  province  are  readonly  dis¬ 
tributed,  any  anisotropy  present  would  con¬ 
tribute  to  a  larger  estlaated  standard  deviation 
of  the  aean  for  the  Fn  neaaurenanta  for  the 


province  but  would  not  be  expected  to  signifi¬ 
cantly  affect  tha  naan  value  of  Fa.  Of  couree, 
for  province a  that  have  only  a  few  observations, 
the  effect  could  be  nor*  pronounced.  Possible 
errors  due  to  neglect  of  consideration  of  the 
spherical  earth  can  b*  evaluated  and  ahown  to 
be  anall.  The  errors  are  due  to  the  fact  that 
■any  analyses  of  selsnic  velocities  utilise 
a  flat-earth  fornulatlon  in  which  apparent  valo- 


£  ± 


TABLE  2.  But  Flow.  F„  Velocity.  Estimated  Temperature  at  the  Mohorovicic 
Discontinuity  and  Crustal  Thickness  Averaged  by  Province  for  North  America 


Includes  Superior,  Grenville,  Churchill,  Slave,  Bear  and  Hudson  Platform  provinces 
Includes  Great  Plains,  Osark  and  Interior  Platfons  provinces. 
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Fig.  1.  Locations  of  refraction  profiles  for  the  continental  crust  and  upper  aantle  of  North 
America.  References  for  the  profiles  are  given  in  Table  1.  Province  boundaries  are  fro*  Fenn email 
[1946]  for  the  United  States  and  fro*  Douglas  and  Price  [1972]  for  Canada. 


dtp  along  tha  surface  of  the  P„  arrivals  in 
the  absence  of  dip  on  the  Hobo  is  equal  to  the 
upper  mantle  velocity.  Consideration  of  a 
spherical  earth  indicates  that  the  true  velo¬ 
city  of  the  upper  aantle  is  slightly  lower  than 
the  apparent  velocity.  The  difference  is  depen¬ 
dent  on  the  velocity  and  the  depth  of  the  re¬ 
fractor  (the  crustal  thickness)  but  varies  from 
-0.03  to  -0.06  km/ a  for  the  range  of  velocities 
and  thicknesses  Included  here.  These  correc¬ 
tions  have  not  been  applied  to  the  PB  data  used 
here  because  publications  that  report  the  Pa 
data  generally  do  not  Include  an  indication 
as  to  whether  spherical-earth  or  flat-earth 
calculations  were  utilised.  Including  this 
correction  would  have  no  significant  effect 
on  the  correlation  analysis  discussed  below. 

Variations  in  the  P„  velocities  due  to  dif¬ 
ferences  in  crustal  thickness  were  corrected 
for  pressure  by  recalculstlng  each  to  a  common 
depth  of  3S  ka.  To  do  this,  the  pressure  gra¬ 
dient  was  assumed  to  be  llthoatatic,  l.e. . 
dP/dZ  *  gp  where  P  is  pressure,  Z  is  depth, 
p  is  density,  and  g  -  980  ca/az,  the  accelera¬ 


tion  of  gravity.  The  variation  of  velocity 
with  depth  will  then  be  dV/dZ  «  (3V/3P)T 
(dP/dZ)  *  (>V/3P)y  gp.  On  the  basis  of  experi¬ 
mental  studies  [Christensen,  1974;  Manghnani 
et  al..  1974],  (3V/3P)T  was  taken  to  be  0.013  ka 
s-1  kbar-1.  Tha  average  crustal  density  was 
taken  as  2.9  ga/ca3  [Verhoogen  et  al.,  1970, 
p.  617].  The  pressure  correction  is  then  4.2  x 
10"3  ka  s"1  ka"1  X  (35-H,.  ka)  where  Hq  Is  the 
average  crustal  thickness  for  the  province  cal¬ 
culated  froa  the  seisaic  data.  Varying  the 
pressure  correction  over  a  fairly  wide  range 
(or  neglecting  the  correction  completely)  does 
not  significantly  affect  the  correlation  analysis 
discussed  below. 

Heat  Plow  Data 

Hast  flow  values  for  the  United  States  pub¬ 
lished  through  1976  were  compiled  by  Sass  et 
al.  [1976].  These  values  were  suppleasnted 
by  data  taken  froa  Henyey  and  Lee  [1976], 

Pollack  and  Chapaan  [1977],  and  Costain  et  al. 
[I960]. 
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TABLZ  3.  Regression  Parameters 

for  P_  ■ 
n 

A  +  Bqo 

and  P_  - 

D 

A  ♦  BT^ 

Case 

Regression 

Variables 

Method 

Province* 

Excluded 

N 

A 

°A 

B(xl04) 

Oj(xl04) 

r 

.  o 

1 

’o  '  Pn 

RMA 

- 

13 

8.420 

0.123 

-68.8 

15.9 

-0.550 

0.10 

2 

«o  '  Pn 

LSQC 

- 

13 

8.386 

0.121 

-68.0 

19.4 

-0.550 

0.10 

3 

%  ‘  pn 

RMA 

SH 

12 

8.487 

0.131 

-75.6 

16.0 

-0.679 

0.05 

4 

S-Pn 

LSQC 

SM 

12 

8.516 

0.113 

-82.9 

17.4 

-0.679 

0.05 

5 

T.  Pn 

RMA 

- 

13 

8.271 

0.097 

-4.11 

0.97 

-0.519 

0.10 

6 

T.-Pn 

LSQC 

13 

8.490 

0.135 

-8.46 

2.15 

-0.519 

0.10 

7 

T.  *Pn 

RMA 

SM 

12 

8.312 

0.102 

-4.40 

1.00 

-0.619 

0.05 

8 

T.-pn 

LSQC 

SN 

12 

8.456 

0.112 

-7.29 

1.70 

-0.619 

0.05 

9 

T.-Pn 

RMA 

SN.SR 

11 

8.370 

0.119 

-5.53 

1.26 

-0.653 

0.05 

10 

T.-Pn 

LSQC 

SN.SR 

11 

8.488 

0.117 

-8.08 

1.88 

-0.653 

0.05 

11 

<o"Pn 

LSQ 

- 

13 

8.223 

0.115 

-37.8 

17.3 

-0.550 

0.10 

12 

qo  '  Pn 

LSQ 

SN 

12 

8.327 

0.120 

-51.3 

17.5 

-0.679 

0.05 

13 

T.  Pn 

LSQ 

- 

13 

8.133 

0.082 

-2.13 

1.06 

-0.519 

0.10 

14 

T.-P» 

LSQ 

SR 

12 

8.189 

0.088 

-2.72 

1.09 

-0.619 

0.05 

15 

T.-Pn 

LSQ 

SN.SR 

11 

8.240 

0.101 

-3.61 

1.40 

-0.653 

0.05 

q  -  boat  flow  averaged  by  province;  PQ  -  upper  mantle  compresalonal  wave  velocity 
averaged  by  province  <!■/•) ;  Ta  -  temperature  eatlaate  at  the  Moho  (*C);  UU  -  reduced 
major  axle  method;  LSQC  -  leaat  equarea  cubic  method;  LSQ  -  atandard  laaat  square* 
linear  regresalon;  SN  -  Sierra  Nevada  province;  SB  -  southern  Rocky  Mountains  province; 
M  -  number  of  data  points;  A  -  Intercept;  oA  -  estimated  atandard  deviation  of .A;  B  - 
slope;  o.  -  estimated  standard  deviation  of  B;  r  -  correlation  coefficient;  a  -  level 
of  significance  for  r. 


Heat  flow  values  were  sorted  by  province, 
and  the  mean  and  estimated  atandard  deviation 
of  the  mean  were  calculated  for  each  province 
(Table  2).  Value*  of  beat  flow  greater  than 
200  sH/m*  were  considered  anomalous  and  ellml-  - 
ns ted.  Such  a  high  value  for  observed  beat 
flow  generally  lapllas  the  action  of  hydro- 
thermal  system*  or  other  near-surface  conditions 
not  representative  of  the  regional  heat  flow. 

Temperatures 

In  order  to  evaluate  the  possible  dependence 
of  the  P„  velocity  on  temperature,  the  tempera¬ 
ture  at  the  Mohorovl£i&  discontinuity  (Moho) 
for  each  province  must  be  calculated.  Such 
thermal  calculations  would  require  additional 
Information  on  the  thermal  conductivity  and 
beat  generation  values  for  each  province.  In 
addition,  assumptions  pertaining  to  thermal 
equilibria*,  the  distribution  of  heat  producing 
elements  and  tha  possible  contributions  of  ther¬ 
mal  convection  to  the  heat  flow  would  have  to 
be  made  for  each  province.  Such  an  analysis 
Is  beyond  the  scop*  of  this  study  and,  therefore, 
we  have  estimated  the  temperature  at  the  Moho' 
beneath  each  province  by  utilising  the  geo¬ 
thermal  gradient  curves  presented  by  Lachenbruch 
end  Sees  (1977 J.  The  temperature  estimates 
were  Interpolated  from  the  curves  of  Lachenbruch 
and  Sesa  (19771,  using  the  average  heat  flow 
and  cruetal  thickness  calculated  here  (Table  2) . 
Temperature  error  estlmetea  are  based  on  the 


estimated  uncertalntiee  In  beat  flow  and  crustal 
thickness  for  each  province  (Table  2). 

Correlation 

The  correlations  of  continental  upper  mantle 
seismic  velocity  (PB)  and  heat  flow  (q0)  as  well 
as  ?B  velocity  and  temperature  estimated  at  the 
Moho  (Tb)  were  evaluated  by  least-squares  linear 
regression  of  the  form  y  ■  A  +  Bx  where  y  Is 
the  ?„  variable  and  x  Is  q.  or  T_.  The  coef¬ 
ficients  A  (Intercept),  B  (slope),  and  the  cor¬ 
relation  coefficient  (r),  as  well  as  standard 
deviations  estimates  of  A  (a*)  and  B  (og),  were 
calculated  and  are  given  In  Table  3.  It  Is 
well-known  that  standard  linear  regression 
methods  (x  considered  the  Independent  variable 
and  y  considered  the  dependent  variable)  yield 
unreliable  estimate*  of  the  A  and  B  coefficient* 
for  the  situation  In  which  both  tb*  x  and  y 
variables  are  subject  to  error.  This  can  be 
easily  demonstrated  by  comparing  the  results 
of  the  regressions  for  y  on  x  and  for  x  on  y. 
Because  both  the  x(q„  or  Ta)  and  y(Pn)  variables 
analysed  her*  are  subject  to  error,  the  're¬ 
duced  major  axis'  (RMA)  and  the  'least-squares 
cubic'  (LSQC)  methods  of  linear  regression  were 
utilised  In  order  to  evaluate  the  correlation 
of  P„  versus  q0  and  Pq  versus  Ta.  In  the  RMA 
method  [Karmack  and  Haldane,  1930],  the  x  and 
y  variables  era  scaled  by  standard  deviations 
of  the  data  and  the  perpendicular  distances 
of  points  to  the  least  square*  line  are  minimised. 
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Fig.  2.  Fn  velocity  and  beat  flow  data  averaged  by  province  for  North  America.  The  pro¬ 
vince  abbreviations  and  the  data  ara  given  in  Table  2.  Error  bars  correspond  to  plus  and 
minus  two  estimated  standard  error*  of  tha  mean  for  tha  P„  and  heat  flow  observations. 

The  best  fitting  regression  line  for  the  case  of  the  Sierra  Nevada  province  excluded  is 
shown  along  with  the  regression  parameter a  (case  4,  Table  3).  Data  from  Australia  [Cull 
and  Denham,  1979]  are  shown  by  triangles  with  error  bar*  (for  hast  flow  only)  indicating 
plus  and  minus  one  standard  deviation.  Each  of  the  Australian  data  points  corresponds 
to  an  individual  geological  province.  The  data  from  Australia  were  not  used  in  the  re¬ 
gression  line  calculation. 


In  the  LSQC  method  [York,  1966,  1967],  each 
observation  la  weighted  according  to  ita  esti¬ 
mated  error  In  both  the  x  and  y  directions. 

The  solution  of  a  cubic  equation  yields  the 
coefficients  A  and  I  and  estimates  of  their 
standard  deviations.  Because  each  data  point 
Is  individually  weighted,  and  because  errors 
in  both  variables  are  considered,  this  method 
has  the  advantage  that  data  points  with  the 
smallest  errors  have  the  greatest  effect  on 
the  regression  line.  The  regression  paranstara 
for  the  correlations  of  Pn  velocity  and  heat 
flow  and  for  P„  velocity  and  temperature  are 


given  in  Table  3.  Regression  parameters  for 
standard  least  squares  fit  to  the  q^F-  and 
T^-Pjj  data  are  shown  for  comparison  with  the 
IMA  and  LSQC  methods  as  cases  11  to  15  In  Table 
3.  The  correlation  coefficient  r  and  the  asso¬ 
ciated  level  of  significance  a  are  determined 
from  the  linear  correlation  of  the  observation 
without  regard  to  estimated  errors  in  the  data 
and,  therefore,  are  the  same  for  the  three 
methods . 

Tbs  PB  velocity  data  are  plotted  versus  heat 
flow  and  estimated  temperature  at  the  Mo  ho  in 
Figures  2  and  3.  Error  bars  for  the  mean  F_ 
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Fig.  3.  F.  velocity  and  estimated  temperature  at  the  Mo ho  data  averaged  by  province 
for  North  America.  The  province  abbreviations  and  data  ara  given  in  Tabl*  2.  Error 


bars  correspond  to  plus  and  minus  two  estimated  standard  deviations  of  tha  mean  cal¬ 
culated  for  the  PB  data  and  estimated  error  limits  for  the  temperatures.  Tha  best 
fitting  regression  line  for  the  case  of  the  Sierra  Nevada  province  excluded  ia  shown 
along  with  the  regression  parameters  (case  8,  Table  3). 
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TABLE  4.  Comparison  of  Temperature  Darlvatlvas  Inferred  Prom 
P  -  T  Regression  with  Laboratory  Determinations  of 
(3Vp/JT)p  for  Dltramaflc  Rocks 


*p  4 

(3T>p  *  10 

km/a/*C 


Rsfsrenca 


-4.12 

-4.94 

-5.6 

-6.1 

-4.4 

-6.5 

-6.7 

-5.6 

-7.1 

-6.7 

-6.2 

-6.9 


Regression  of  Pn  and  T  [this  study] 

Dunite  (Kara  and  Richter.  1981] 

Paridotlta  (Kara  and  Richter,  1981] 

Dunite  (Christensen,  1979] 

Dunite  [Ramananantoandro  and  Manghnanl,  1978] 
Perldotite,  data  from  Flelltz  [1976] 

Harcburgita  [Peselnick  and  Nicolas,  1978] 

Har  zburgite  [Peselnick  and  Nicolas,  1978] 

Har zburgite  [Peselnick  and  Nicolas,  1978] 

Lhercolite  [Peselnick  and  Nicolas,  1978] 
Lherzolite  [Peselnick  and  Nicolas,  1978] 
Harsburgice  [Peselnick  at  al.,  1977]  (corrected 
values  by  Peselnick  and  Nicolas  [1978]) 

Her zburgite  [Peselnick  at  al.,  1977]  (corrected 
values  by  Peselnick  and  Nicholas  [1978]) 


velocity  and  heat  flow  values  are  scaled  to 
plus  and  minus  two  estimated  standard  devia¬ 
tions  of  the  aean  that  correspond  approximately 
to  95Z  confidence  Intervals.  The  error  bars 
for  temperature  at  the  No ho  are  estimated  by 
the  range  of  temperatures  read  from  the  geo¬ 
thermal  gradient  curves  of  Lachenbruch  and  Sass 
[1977],  corresponding  to  the  9SZ  confidence 
Interval  ranges  of  observed  beat  flow  and 
crustal  thickness.  The  best  fitting  regres¬ 
sion  lines  for  the  case  of  the  Sierra  Nevada 
province  excluded  are  also  shown  In  Figures 
2  and  3. 

If  the  Sierra  Nevada  data  point  Is  excluded 
from  the  least  squares  calculations,  the  cor¬ 
relations  between  heat  flow  and  Pn  velocity 
and  between  temperature  and  Pn  velocity  are 
statistically  significant  at  the  a  -  0.05  level 
of  algnlflcance  (cases  3,  4,  7,  8,  9  and  10, 
Table  3).  The  Slerre  Nevada  data  point  is  one 
of  the  most  discrepant  points  on  the  plots 
shown  In  Figures  2  and  3  and  Its  exclusion  from 
the  regression  calculations  is  based  on  the 
fact  that  tha  estimated  mantle  heat  flow  for 
the  SlerTa  Nevada  Is  anomalously  low  [Blackwell, 
1971]  and  thus  the  estimated  temperature  at  the 
Mo ho  baaed  on  average  surface  heat  flow  nay 
be  substantially  in  error.  The  estimated 
temperature  at  the  Mo ho  for  the  Southern  Rocky  . 
Mountains  (SR)  province  also  appears  highly 
anomalous  (Figure  3).  The  high  value  of  In¬ 
ferred  te^ersture  Is  due  to  the  large  average 
crustal  thickness  of  the  SR  Province  and  the 
relatively  high  observed  heat  flow  (Table  2). 
Such  a  high  temperature  should  not  exist  at 
the  Mo ho  because  it  Is  above  the  melting  points 
for  likely  mantle  and  lower  crustal  rocks. 
However,  because  of  the  large  error  bars  asso¬ 
ciated  with  the  SR  data  point.  Its  Inclusion 
In  the  regression  analysis  has  little  effect 
on  the  resulting  slope  estimates  (compere  cases 
7  and  9  or  8  and  10  from  Table  3).  Although 


the  correlations  of  q,,  and  P„  or  T_  and  P„  are 
statistically  algnlficant.  the  differences 
between  the  intercept  and  alope  coefficients 
for  the  various  regressions  shown  in  Table  3 
(cases  1-4  for  q  and  P  ;  cases  5-10  for  T 
and  P  )  are  not  significant.  ® 

U 

Laboratory  Data 

To  determine  whether  the  regional  variations 
In  te^ierature  calculated  from  heat  flow  are 
sufficient  to  explain  the  entire  variation  In 
continental  Pn  velocities,  experimental  values 
of  (3Pp/3T)p  must  be  considered.  The  regres- 
sion  line  values  found  In  this  study  -4.4  x  10~* 
to  -8.1  x  10”*  km/a/*C  compare  well  with  the 
range  of  experimental  measurements  of  the  tem¬ 
perature  derivatives  of  compressional  velocity 
for  ultraaaflc  rocks  as  shown  In  Table  4.  Al¬ 
though  other  factors  could  also  have  a  signifi¬ 
cant  affect,  these  results  are  consistent  with 
tha  hypothesis  that  temperature  effects  could 
explain  tha  entire  observed  variation  In  ?n 
velocity  between  provinces. 

Discussion 

This  study  has  demonstrated  a  statistically 
significant  relationship  between  continental 
Fn  velocity  and  heat  flow  on  a  regional  scale 
for  North  America.  Heat  flow  has  also  been 
related  to  crustal  aga  [Polyak  and  Smirnov, 

1968]  and  to  thicknsHS  of  the  lithosphere 
[Crough  and  Thompson,  1976;  Pollack  and 
Chapman,  1977;  Kono  and  Amano,  1978].  Figure 
4  shows  the  observed  relationship  between  Pn 
velocity,  crustal  age,  and  heat  flow  which 
was  obtained  by  using  the  relationship  between 
heat  flow  and  age  determined  by  Polyak  and 
Smirnov  [1968]  (as  plotted  by  Chapman  and 
Pollack  [1975])  and  the  relationship  between 
heat  flow  and  Pn  velocity  found  In  this  study. 
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Fig.  4.  Relationship  of  obsarvad  haat  flow,  cruatal  aga  and  Fn  valoclty.  Heat  flow 
aga  data  are  fro*  Polyak  and  Snlrnov  {1968]  aa  plotted  by  Chapman  and  Pollack  [1975] 
for  continents.  Pn  velocity  -  haat  flow  relation  la  Inferred  by  the  linear  regression 
for  North  American  data  shown  In  Figure  2.  Crustal  age  Is  Interpreted  as  tine  since 
the  last  t hansel  event  affecting  the  continental  lithosphere.  Lower  schematic  dia¬ 
gram  Illustrates  the  evolution  of  the  continental  lithosphere  and  upper  mantle  velo¬ 
city  as  a  function  of  time  (cooling  of  the  lithosphere) .  Pn  values  are  In  ka/s  Infer¬ 
red  fron  the  q0  “  *n  -  ags  curve  above  and  the  Ta  estimates  (In  *C)  are  Inferred  froa 
the  Pn  -  Ta  relationship  Illustrated  In  Figure  3.  The  depth  to  the  base  of  the  litho¬ 
sphere  is  shown  baaed  on  the  depth  to  the  estlaated  teaperstura  of  partial  melting  of 
upper  aantle  materials  froa  the  model  of  Pollack  and  Chapman  [1977],  The  dotted  line 
In  the  lower  part  of  the  figure  shows  the  llthosphere-asthenorphere  boundary  according 
to  the  relation  determined  by  Kono  and  Aaano  [1978]. 


Although  correlations  of  beat  flow  age  and 
beat  flow  P  velocity  do  not  necessarily  Imply 
a  causal  relationship  between  age  and  PQ  velo¬ 
city,  the  Pn  versus  age  estimates  for  provinces 
In  North  America  shown  In  Figura  4  suggest 
that  this  relationship  Is  valid.  Lithospheric 
thickness  estimates  from  Pollack  and  Chapman 
[1977]  and  Kono  and  Amano  [1978]  are  also  shown 
In  Figure  4.  These  results  suggest  that  thick¬ 
ening  of  the  lithosphere  and  aa  Increase  In 
upper  aantle  P„  velocity  are  related  processes 
caused  by  cooling  of  the  continental  litho¬ 
sphere  with  time  after  a  thermo-tectonic  event. 
The  observed  range  of  continental  upper  mantle 
velocities  («7.6  -  8.3  ka/s)  can  be  explained 
aa  primarily  the  results  of  regional  differences 
In  temperature  at  the  Moho  discontinuity. 

Alternative  explanations  for  the  beat  flow, 

P_  and  crustal  age  data  are  (1)  that  Pn  velo¬ 
city  Is  a  function  of  composition  that  may  vary 
with  age  (for  example,  by  loss  of  volatiles  or 


differentiation  by  magmatic  processes)  and,  (2) 
that  variations  In  Pn  velocity  are  due  to  aniso¬ 
tropy  In  which  the  degree  of  orientation  of 
anisotropic  mineral  grains  varies  vlth  tlae. 
Rybach  and  Buntebarth  [1982]  have  recently  pre¬ 
sented  data  to  indicate  that  seismic  velocity 
and  heat  generation  in  crystalline  rocks  are 
Inversely  related  and  thus  variations  In  com¬ 
position  (and,  therefore,  heat  generation) 
could  be  a  controlling  factor  in  variations 
In  upper  aantle  seismic  velocity.  However, 
the  agreement  of  (3P„/3T)p,  estimated  fron 
the  regression  of  Ta  and  Pn,  with  experimentally 
determined  temperature  derivatives  for  ultra- 
mafic  rocks  suggests  that  the  entire  variation 
may  be  due  to  teaperature  offsets. 

Data  from  North  America  were  selected  for 
the  analysis  of  the  relation  between  PQ  velo¬ 
city  and  heat  flow  because  the  volume  and 
geographical  distribution  of  data  were  suf¬ 
ficient  to  permit  statistical  analysis.  A 
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similar  PB  versus  beat  flow  relationship  (Fig¬ 
ure  2)  was  suggested  by  Cull  and  Oenhaa  [ 1979] 
for  Australia.  Although  no  detailed  analysis 
has  been  performed ■  a  brief  examinee Ion  of  the 
PQ  and  heat  flow  data  for  Europe  and  Asia  sug¬ 
gests  that  a  similar  relationship  will  be  found 
for  these  continental  regions. 
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ABSTRACT.  Seismic  energy  propagating  through  vertically  and 
laterally  varying  structures  of  the  earth's  crust  and  lower 
lithosphere-uppermost  mantle  is  responsible  for  the  numerous  and 
complex  seismic  phases  observed  on  short-period  seismograms  at 
regional  distance  ranges  (100  to  2000  km).  Recent  advances  in 
techniques  for  computing  synthetic  seismograms  make  it  practical 
to  calculate  complete  seismograms  Chat  realistically  model  many 
features  of  regional  phases.  A  modified  reflectivity  method 
program  is  used  to  interpret  some  details  of  record  sections  of 
Nevada  Test  Site  (NTS)  underground  explosions  that  were  observed 
700  to  800  km  from  the  sources. 


I.  INTRODUCTION 

Regional  seismic  phases  recorded  by  high-gain,  short-period  or 
broadband  instruments  are  likely  to  play  an  increasingly  im¬ 
portant  role  in  seismic  source  location  and  identification  as 
acceptable  magnitude  thresholds  are  pushed  to  lower  levels. 
From  the  standpoint  of  complexity  of  seismograms,  the  epicentral 
distance  range  between  *200  km  and  the  transition  to  simpler 
teleseismic  waveforms  around  2000  km  presents  many  challenges  to 
the  seismic  analyst.  In  this  range,  propagation  paths  can 
traverse  the  crust,  the  lower  lithosphere,  and  the  uppermost 
mantle  where  both  vertical  and  lateral  heterogeneities  strongly 
influence  waveform  characteristics.  Good  observational  data  are 
rare  for  testing  analysis  techniques  developed  for  regional 
problems.  In  contrast  to  the  numerous  detailed  crustal  refrac¬ 
tion/reflection  profiles  that  have  been  obtained  from  many  parts 
of  the  world  out  to  distances  -200  km,  relatively  few  long-range 
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profiles  exisc  where  station  spacing  is  sufficiently  tight  to 
facilitate  a  clear  interpretation  of  the  onset,  development,  and 
amplitude  vs.  distance  behavior  of  the  many  observable  phases. 
Thus,  although  signals  from  sources  of  interest  may  be  easily 
observable  at  regional  distances,  derivation  of  source  para¬ 
meters  from  observations  at  sparsely  located  observatories  or 
arrays  will  require  careful  analysis  and  modeling  of  the  in¬ 
tricacies  of  wave  propagation  at  these  scales. 

Phases  of  interest  in  regional  identification  studies  fall 
into  two  main  categories*.  large  amplitude,  long  duration,  but 
somewhat  indistinct  wave  groups  such  as  Lg  and  P;  and  body  waves 
(mainly  compress ional)  that  appear  either  as  first  arrivals  or 
closely  following  as  possible  wide  angle  reflections/near- 
critical  refractions  from  interfaces  and/or  steep  velocity 
gradients  in  the  deep_  crust,  lower  lithosphere,  and  uppermost 
mantle.  The  Lg  and  P  phases  are  often  the  largest  amplitude 
features  on  regional  short_-period  seismograms,  but  a  clear 
explanation  of  how  Lg  and  P  propagate  is  still  lacking  [1]; 
this  lack  perhaps  is  reflected  in  the  fact  that  seismologists 
frequently  use  the  notations  P  or  Pg  interchangeably  in 
reference  to  a  broad,  large  amplitude  phase  following  Pn.  We 
adopt  the  P  notation  here.  The  phase  in  question  propagates 

very  well  in  the  western  United  States,  but  attenuates  rapidly 
in  the  eastern  U.S.  A  group  velocity  around  6  km/s  implies  P 
propagates  as  coapressional  waves  multiply  reflected  within  the 
crust — which  may  thus  act  as  a  waveguide.  Similarly,  the  -3.5 
km/s  group  velocity  for  Lg  suggests  shear  waves  multiply  re¬ 
flecting  within  the  crustal  layers.  Some  authors  [2]  prefer  to 
treat  Lg  as  a  superposition  of  higher  mode  Love  and  Rayleigh 
waves  propagating  in  a  nearly  laterally  homogeneous,  vertically 
layered  crust.  In  any  case,  the  propagation  physics  is  com¬ 
plicated  and  will  require  quite  sophisticated  synthetic  seismo¬ 
gram  codes  to  properly  model  and  interpret  observed  waveforms. 

Record  sections  of  long-range  seismic  refraction  profiles 
often  show  one  or  more  nearly  parallel  travel  time  (T)  vs.  dis¬ 
tance  (d)  branches  following  within  several  seconds  of  first 

arrivals  (3,  4,  5].  Each  secondary  branch  may  be  traceable  only 
over  a  distance  interval  of  SO  to  200  km  before  being  replaced 
in  a  "shingle-like"  fashion  with  another  branch  or  set  of  ar¬ 
rivals  (5,  6,  19].  These  are  usually  interpreted  as  parts  of 
cusp  phases  arising  from  critical  refractions  and/or  wide-angle 
reflections  from  first  order  discontinuities  or  steep  velocity 
gradients  in  the  upper  mantle.  Archambeau  et  al .  [7]  and 

Burdick  and  Relmberger  [8],  for  example,  have  derived  velocity 
vs.  depth  tndels  for  the  major  features  of  the  upper  mantle 

beneath  the  U.S.  by  a  joint  analysis  of  travel  times,  amplitude 
vs.  distance  variations,  and  waveform  fitting  of  the  first  few 
compress ional  arrivals  .observed  at  widely  separated  seismograph 


SEISMOGRAMS  OF  EXPLOSIONS  AT  REGIONAL  DISTANCES  IN  THE  US. 


4SS 


stations  throughout  tha  U.S.  Theae  and  similar  model*  by  othera 
are  moat  valid  for  depths  graatar  than  about  250  km.  Although 
thaaa  analyaaa  auggeat  that  tha  main  faaturea  of  mantle  atructura 
at  daptha  below  about  300  ka  (corresponding  to  coapreasional 
first  arrivals  at  epicentral  ranges  beyond  -1500  km)- may  be  more 
uniform  over  a  global  scale  [8],  it  is-  known  that  significant 
lateral  variations  in  lower  lithosphere  and  uppermost  mantle 
properties  occur  beneath  the  continents  on  regional  and  perhaps 
even  finer  scales  (8,  9,  10,  111.  In  the  depth  range  between 
the  Moho  and  -300  ka,  several  types  of  structural  variations 
have  been  suggested  in  the  literature  that  would  give  rise  to 
wide  angle  reflections,  converted  phases,  and  similar  closely 
spaced  arrivals  on  seismograms  at  regional  ranges.  /  These  in¬ 
clude  the  presence  or  absence  of  the  S-wave  and/or  the  P-wave 
low  velocity  cone  (LVZ)  in  the  asthenosphere,  high  velocity 
mantle  lids  (12,  13],  alternating  lamellae  of  positive  and 

negative  velocity  gradients  (6,  19],  etc.  These  early  arriving 
phases  often  have  better  defined  onsets  than  the  P  and  Lg  phases 
and,  since  they  are  observed  at  distances  beyond  that  where  a 
true  head  wave  Pn  arrival  can  be  expected,  they  may  be  useful  in 
regional  source  location  and  identification.  In  order  to  make 
use  of  the  information  contained  in  these  arrivals  (especially 
the  amplitude  vs.  distance  behavior  for  particular  paths  of 
interest),  it  will  be  necessary  to  use  modern  sophisticated 
synthetic  seismogram  techniques  to  derive  localised  fine  scale 
details  frcm  generalised  crust-mantle  models. 

The  purpose  of  this  paper  is  to  explore  a  few  of  the  prob¬ 
lems  in  modeling  regional  short-period  seismograms  by  means  of  a 
modified  reflectivity  method  [14]  computer  program  developed  by 
R.  Kind  [15].  This  numerical  program  accounts  for  the  effects 
of  a  buried  source  and  is  thus  capable  of  computing  ’complete' 
seismograms — including  refracted  waves,  surface  reflected  body 
waves  such  as  the  pP  phase,  and  surface  waves.  The  effects  of 
anelastic  attenuation  (Q)  for  each  layer  are  included  as  an 
integral  part  of  the  method  [15].  The  most  severe  limitation  of 
the  technique  for  studies  of  regional  seismograoia  is  the  assump¬ 
tion  of  lateral  homogeneity  (this  is  also  a  limitation  for  nor¬ 
mal  modes  summation  techniques).  An  item  of  interest  will  be 
the  extent  synthetics  can  be  siade  to  match  observed  wsve forms 
under  this  restriction. 

Two  problems  are  considered.  The  first,  labeled  the  B-3 
model  for  brevity,  employs  a  simple  model  consisting  of  three 
layers  in  the  crust  without  velocity  gradients  and  an  almost 
uniform  velocity  mantle.  A  large  range  of  apparent  surface 
phase  velocities  is  used  in  order  to  display  S  phases  and  sur¬ 
face  waves.  The  second  calculation,  the  A-10  model,  treats  the 
mantle  structure  in  detail,  but  confines  sttention  to  compres¬ 
sions!  phases  near  ther  start  of  the  seismogram.  The  more 
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Pig.  1.  (a)  Location  map  of  the  western  United  States  with 
relative  positions  of  the  Nevada  Test  Site  and  the  Y-ESRP 
recording  line.  (b)  Enlargement  showing  positions  of  stations 
that  recorded  the  27  September  1978  RUMMY  explosion.  Asterisk 
denotes  approximate  area  for  mantle  ray  turning  points  from  NTS 
explosions. 


important  conclusions  of  the  A-10  model  are  summarised  here — a 
fuller  discussion  of  this  calculation  and  the  impl icat ions  for 
uppermost  mantle  structure  beneath  the  western  U.S.  can  be  found 
in  a  previous  publication  [16] 1 

A  comparison  of  the  synthetic  seismogram  calculations  has 
been  made  with  a  100-km-long  record  section  of  short-period  ver¬ 
tical  component  seismograms  obtained  in  eastern  Idaho  during  the 
1978  Yellowstone-Eastern  Snake  River  Plains  (Y-ESRP)  seismic 
profiling  experiment,  for  these  observations,  the  sources  were 
underground  nuclear  explosions  at  the  Nevada  Test  Site  (NTS)  at 
distances  between  720  and  820  km  from  the  nearly  radially  ori¬ 
ented  linear  station  array  (Fig.  1).  Only  the  records  from  the 
largest  NTS  explosion,  the  %  “  5.7  ROMfY  event  at  1720»00.076 
OfT,  27  September  1978,  are  reproduced  here  since  they  have  the 
best  signai-to-noise  ratio  of  the  three  NTS  explosions  observed 
during  the  experiment.  Additional  details  of  the  Y-ESRP  instru¬ 
mentation,  experiment,  and  data  can  be  found  elsewhere  [16]. 
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2 .  COMPUTAT IONAL  TECHNIQUE 

As  discussed  by  Kind  [15]  snd  by  Fuchs  and  Muller  (14),  the 
reflection  coefficient  and  tine  shift  calculations  in  the  re¬ 
flectivity  method  are  carried  out  in  the  frequency  domain  and 
then  Fourier  transformed  to  plot  seismograms.  We  included 
Miller's  [17]  earth  flattening  approximation  in  both  of  our 
problems  to  account  for  earth  curvature  effects.  Both  P  and  S 
velocities  are  independently  specified  in  all  calculations, 
since  the  reflection  coefficients  are  functions  of  both  P  and  S 
velocity  contrasts  at  non-normal  incidence  angles  and  are 
required  even  when  only  computing  P  phases  over  a  narrow  time 
window.  In  the  A-10  calculation,  for  example,  the  departure  of 
the  P/S  velocity  ratio  in  a  layer  from  that  given  by  /Poisson's 
ratio  -  1/4  is  an  important  factor  in  our  interpretation  [16], 
Densities  are  given  by  a  Birch's  Law  relation  (density  *  0.252  + 
0.3788*P  velocity).  The  attenuation  factor  Qa  for  P  waves  was 
chosen  as  25  in  the  source  layers,  200  in  the  upper  crust,  and 
1000  in  the  lower  crust  and  the  uppermost  mantle  layers;  for  the 
LVZ  modeling  of  the  A-10  model,  Qa  in  the  asthenospheric 
layers  was  adjusted  as  part  of  the  fitting  procedure  (see 
Fig.  5).  The  attenuation  factor  for  S  waves  was  always  assumed 
to  be  4Qa/9  [20].  The  explosive  source  algorithm  [16]  was  used 
with  the  source  buried  at  a  depth  of  0.640  km  in  a  layer  of  P 
velocity  *  3.55  km/s.  These  were  close  to  actual  field  values 
for  the  NTS  RUMMY  explosion.  Time  intervals,  number  of  samples, 
and  computed  lengths  of  seismograms  were  chosen  so  that  the 
dominant  frequency  of  the  source  spectrum  was  1.6  Hs  for  the 
A-10  calculation — again  close  to  the  observed  value.  In  order 
to  save  computer  time  for  the  extended  duration  B-3  seismogram 
sections,  the  parameters  were  chosen  so  that  the  dominant  fre¬ 
quency  of  the  source  was  shifted  to  0.25  Hz;  although  this  was 
lou  compared  to  observed  frequencies,  we  felt  it  was  adequate 
for  the  puposes  of  this  initial  study.  To  avoid  long  computer 
runs,  the  wave  field  was  only  computed  within  a  limited  phase 
velocity  window;  1  km/s  to  20  km/s  for  B-3,  and  6.5  km/s  to 
1000  km/s  for  A-10.  These  integration  limits  sometimes  in¬ 
troduced  spurious  single  cycle  "phases"  at  these  apparent  velo¬ 
cities  in  the  computed  record  sections.  The  limit  velocities 
were  chosen  so  as  to  not  overlap  or  interfere  with  arrivals  of 
interest  in  the  observations.  In  the  record  section  plots,  the 
amplitudes  of  each  trace  have  been  multiplied  by  station  dis¬ 
tance  to  maintain  a  convenient  scaling  of  the  amplitudes  of  the 
phases  which  are  subject  to  geometrical  spreading  and  attenua¬ 
tion  due  to  anelasticity. 
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Fig.  2.  Vertical  component  low  time  resolution  seismic  record 
section  o£  the  KUMfY  explosion  as  recorded  at  Snake  River  Plains 
stations.  The  time  scale  is  compressed  to  show  envelope  be¬ 
havior;  individual  waveforms  not  readily  seen.  The  P  phase  is 
the  broad  feature  at  reduced  times  between  30  and  60  » .  Upward 
ground  motion  to  che  left. 


3.  DISCUSSION 

3.1  The  Extended  Time  Seismograms;  B-3  Model 

Figure  2  is  a  true  relative  amplitude  vertical  component  record 
section  of  the  RUMfY  explosion  recorded  on  ten  matched  short- 
period  (l*Hx  natural  frequency)  instruments  deployed  in  the 
eastern  Snake  River  Plains  (Fig.  1).  Although  Che  time  scale  is 
Coo  compressed  to  reveal  many  details  of  che  waveforms,  several 
important  overall  features  can  be  noted.  The  broad  (-40-second- 
long)  envelope  of  the  ?  phase  appears  at  reduced  times  between 
approximately  30  to  60+  seconds,  and  is  the  largest  amplitude 
feature  on  the  record.  In  contrast,  the  Lg  phase  expected  at 
reduced  times  of  -130+  seconds  (an  average  velocity  of  about  3.5 
km/s)  is  poorly  developed  on  these  unfiltered  records;  it  is 
only  obvious  at  the  770-km  station.  A  few  impulsive  arrivals 
can  be  seen  (such  as  the  first  arrivals  at  reduced  time  -10 
seconds,  idtich  will  be  discussed  in  Sec.  3.2,  and  perhaps  an 
Sn  (?)  phase  at  trt4*80  seconds  and  -780  km),  but  the 
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impression  one  gets  by  viewing  this  observed  section  is  thet  the 
correlations  seem  to  be  better  described  as  broad  energy 
correlations  rather  than  phase  correlations.  A  similar  con¬ 
clusion  is  suggested  by  seismograms  from  central  Asia  shown  in 
the  paper  of  Ruzaikin  et  al.  (1].  A  coherent  structure  in  the  P 
and  Lg  phases  is  difficult  to  trace  from  station  to  station  even 
though  the  stations  are  only  separated  by  8  km  on  the  average. 

The  results  of  an  attempt  to  model  late  time  arrivals  over 
a  regional  distance  range  is  shown  in  Fig.  3.  A  rudimentary, 
almost  trivial,  crust/mantle  velocity  structure  was  assumed  that 
consisted  of  three  constant  velocity  layers  in  the  crust  over¬ 
laying  a  nearly  constant  velocity  half space.  (A  slight  negative 
gradient  in  P  velocity  was  introduced  just  below  the  Moho  in 
order  to  suppress  the  Pn  amplitudes  as  required  by  the  observa¬ 
tions;  see  Sec.  3.2.)  We  note  several  points. 

(a)  The  seismogram  section  from  100  to  900  km  and  the  enlarged 
individual  record  for  800  km  shows  a  surprising  amount  of 
complexity  at  times  beyond  the  first  arrivals  even  though 
an  extremely  simple  earth  model  and  source  function  is 
used.  Groups  corresponding  to  the  P  and  Lg  phases  can  be 
identified. 

(b)  There  appears  to  be  a  considerable  amount  of  S-wave  energy 
although  none  is  present  in  the  explosion  source  algorithm. 
This  is  probably  due  to  P-to-S  and  S-to-P,  etc.,  conver¬ 
sions  at  interfaces  and  to  multiples  which  the  program 
adequately  includes. 

(c)  The  calculated  dispersed  fundamental  mode  Rayleigh  wave  is 
very  large.  There  are  at  least  two  reasons  this  Rayleigh 
wave  is  not  representative  of  the  observations.  First,  no 
corrections  for  the  short-period  bandpass  response  of  the 
seismometers  were  included  in  the  synthetics.  Second,  the 
assumed  source  spectrum  has  too  much  energy  at  the  longer 
periods  as  compared  with  a  near  point-source  representative 
of  a  NTS  explosion,  thus  over  enhancing  the  Rayleigh  waves. 
Long-period  Rayleigh  waves  from  actual  underground  explo¬ 
sions  are  probably  generated  or  modified  and  enhanced  by 
mechanisms  such  as  spall  closure  and/or  tectonic  strain 
release;  these  mechanisms  are  not  adequately  treated  by  the 
explosion  algorithm  used  for  the  present  calculation. 

(d)  Because  the  calculated  seismogram  sections  are  quite  com¬ 
plicated  even  for  this  simple  earth  model,  they  give  the 
impression  that  broad  "packets  of  energy"  can  be  more 
readily  corjelated  than  any  well  defined-  phases—for  at 
least  the  P  and  Lg  phases.  This  was  the  case  with  the 
observations  in  Fig.  2.  In  order  to  better  understand  the 
gross  behavior  of  these  phases  with  distance  and  to 
identify  the  origin  of  obscure  features,  it  will  be 
necessary  to  include  calculations  of  the  horizontal 
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Fig.  3.  (a)  Synthetic  seisnograa  vertical  component  record 

eection  calculated  froa  the  P  and  S  velocity  va.  depth  structure 
(Model  8-3)  shewn  in  (b).  (c)  Expanded  plot  of  the  synthetic 

seismogram  at  the  800-ka  distance.  Approxiaate  arrival  tiae  and 
average  velocity  windows  for  different  phases  or  groups  are 
indicated;  the  phase  velocities  of  the  different  wave  types  are 
equal  to  or  slightly  greater  than  tha  average  velocities.  The 
Rayleigh  waves  on  plot  (a)  are  arbitrarily  clipped  in  plotting 
to  avoid  large  overlays  in  the  seisaograas . 


SEISMOGRAMS  OF  EXPLOSIONS  AT  REGIONAL  DISTANCES  IN  THE  U.S. 


461 

(radial)  component  and  to  perform  calculations  at  small 
station  separation  to  increase  recognizability  of  phase 
correlations. 

These  results  suggest  that  the  modified  reflectivity  method, 
even  with  the  restrictive  assumption  of  lateral  homogeneity,  can 
be  a  useful  technique  in  understanding  the  intricacies  of  Lg  and 
P  phases  and  the  types  of  earth  structures  that  most  affect  them. 
In  addition,  these  studies  suggest  that  observations  of  complex 
and  apparently-incoherent  seismic  phase  arrivals — even  over  short 
distances — do  not  necessarily  imply  strong  lateral  heterogeneity 
in  crustal  structure.  Parameter  studies  would  help  identify 
those  aspects  where  refinements  due  to  lateral  heterogeneity 
and/or  scattering  need  to  be  considered  in  order  to  better  match 
observat ions . 

3.2  Early  Time  Arrivals;  A-10  Model 

Figures  4a  and  4b  are  enlarged  portions  of  the  first  few 
seconds  of  the  digitized  RUMMY  vertical  component  seismograms 
(see  also  Fig.  2)  that  show  details  of  the  earliest  arrivals. 
We  have  interpreted  [16]  this  record  section  in  terms  of  three 
different  compressions!  phases,  all  having  apparent  velocities 
close  to  8  km/s:  (a)  an  extremely  weak  leading  arrival  labeled 
Pn,  which  was  lost  in  the  background  noise  for  the  two  other, 
lower  yield,  NTS  shots  that  were  also  recorded  during  the  Y-ESRP 
experiments;  (b)  a  stronger  phase  labeled  follows  Pn  by 

about  two  or  three  seconds  for  epicentral  distances  between  700 
and  780  km;  (c)  beyond  780  km,  the  Piid  phase  appears  Co  be 
overtaken  and  overwhelmed  by  a  low-frequency  phase,  P;  ,  whose 
amplitude  increases  rapidly  with  distance  out  to  at  least  the 
farthest  station  of  the  linear  array.  The  detailed  reasons  for 
these  labels  and  identifications  are  discussed  in  [16];  they  can 
be  sunsnarized  as  follows. 

The  phase  labeled  Pn  could  be  a  wide  angle  reflection  from 
a  weak  P-velocity  contrast  in  the  lower  lithosphere  below  the 
Moho  rather  than  a  true  headwave  (in  the  strict  sense  of  the 
mathematical  definition)  that  travels  along  the  M-discontinuity 
interface  over  the  entire  800-km  path.  However,  the  sub-Moho 
P  velocity  (7.7  to  7.9  km/s)  in  this  region  of  the  Great  Basin 
is  known  to  be  close  to  both  the  average  and  the  apparent  ve¬ 
locity  observed  in  Figs.  2  and  4.  This,  plus  the  fact  that 
other  travel  time  arguments  [16]  suggest  there  is  no  evidence 
for  mantle  lids  or  other  thin  but  fairly  high  gradient  zones 
down  to  a  depth  of  about  100  km,  argues  that  the  most  straight¬ 
forward  explanation  for  this  arrival  is  that  it  is  a  Pn-type 
phase.  We  calculate  that  the  energy  at  800  km  is  greatly  re¬ 
duced  because  the  wave  travels  in  a  region  beneath  the  Moho  that 
has  a  slight  negative  velocity  gradient. 
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Fig.  4.  U)  True  relative  amplitude  record  section  of  early 
compress ional  arrivals  from  the  RUMMY  explosion.  (b)  Same  as 
(a)  with  increased  amplitudes  to  show  weak  Pn  phase.  Upward 
motion  to  the  left.  All  traces  low  pass  filtered  at  3  Hz. 


The  sudden  onset  at  about  780  km  and  subsequent  rapid 
amplitude  growth  of  the  ?i  phase  indicates  it  is  the  cusp  of 
the  critically  refracted  P-waves  from  the  steep  velocity  gradi¬ 
ent  at  the  base  of  the  asthenospher ic  low  velocity  zone.  The 
observed  dominant  low  frequency  content  is  Chen  explained  by  the 
attenuation  of  the  high  frequency  components  as  the  energy 
travels  first  downward  and  then  back  up  through  Che  very  lov-Q 
region  of  the  LVZ.  The  notation  of  P}  for  this  phase  follows 
the  convention  established  by  Archambeau  et  al .  (7). 
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The  travel  times,  moderate  amplitudes,  and  relatively  high 
frequency  content  imply  the  phase  identified  as  Pj^  is  a  vide 
angle  reflection  from  a  discontinuity  near  the  base  of  the 
mantle  lid  (*  top  of  LVZ)  in  this  area. 

The  conclusions  concerning  these  three  early  arriving  coa- 
pressional  phases  summarized  above  were  confirmed  by  using  the 
modified  reflectivity  program  to  quantitatively  model  the  arri¬ 
val  times,  amplitudes,  and  waveforms  in  the  first  IS  seconds  of 
the  record  sections.  The  procedure  was  to  begin  with  a  generic 
P-velocity  vs.  depth  model  for  the  western  U.S.  (the  T-7  model) 
derived  from  a  wider  data  set  by  Burdick  and  Helmberger  [8]  and 
then  to  perturb  the  model  to  achieve  a  better  fit  [16].  Because 
of  the  influence  of  S-velocity  contrasts  on  the  P-w|ve  re¬ 
flectivity  calculations,  an  S-velocity  vs.  depth  model ''derived 
by  Priestly  and  Brune  (18]  from  an  analysis  of  Rayleigh  and  Love 
wave  dispersion  on  paths  crossing  the  area  of  interest  in  the 
Creat  Basin  of  Eastern  Nevada  was  incorporated  into  the  syn¬ 
thetic  seismogram  modeling.  The  starting  T-7  and  Priestly-Brune 
(P/B)  velocity  models  are  shown  by  dotted  lines  in  Fig.  5.  The 
generic  T-7  P-wave  model  has  a  pronounced  mantle  lid  with  a 
strong  positive  P-velocity  gradient  beneath  the  Moho  for  depths 
from  33  to  65  km.  Calculation  of  synthetics  for  this  lid 
structure  gave  very  large  amplitudes  for  the  "Pn"  arrival, 
which  was  superimposed  on  a  strong  reflection  from  the  base  of 
the  lid  at  65  km  (16).  Thus,  the  T-7  ♦  P/B  starting  model  gave 
results  very  different  from  observations.  However,  as  seen  in 
Tig.  5,  only  small  changes  to  the  initial  model  were  necessary 
to  siatch  the  observations.  To  bring  the  calculated  synthetic 
seismograms  into  agreement  with  observations,  the  gradient  at 
the  base  of  the  LVZ  had  to  be  raised  to  shallower  depths  and  the 
positive  gradient  lid  replaced  with  a  smooth  but  gradual  nega¬ 
tive  gradient  starting  at  the  M-discontinuity.  The  final  model, 
A.-10,  that  matches  observations  is  shown  by  the  solid  lines  in 
Fig.  5.  Figure  6  is  the  comparison  between  the  observed  and 
synthetic  record  sections.  Interestingly,  no  discontinuity  in 
P-velocity  ia  necessary  to  explain  the  ?xid  reflections;  the 
reflections  can  be  adequately  modeled  by  a  small  negative  step 
in  S  velocities  at  a  depth  of  about  100  km.  The  synthetics, 
however,  do  not  seem  to  adequately  model  the  long  oscillatory 
trains  following  the  Pj  phase  onset.  This  is  probably  due  to 
interference  effects  caused  by  fine  structure  in  the  lower  LVZ 
velocity  gradient  that  we  have  not  yet  modeled  by  thin  enough 
layers  in  the  calculation  [16]. 

These  calculations  illustrate  that  synthetic  modeling  tech¬ 
niques  can  be  helpful  in  phase  identification  and  in  quantitative 
calculations  of  amplitude  vs.  distance  behavior  and  waveform 
characteristics.  With  a  sophisticated  reflectivity  method  cal¬ 
culation  we  were  able  to  model  several  important  features  of 
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Fig.  5.  P-velocity  (a)  and  S-velocity  (Q)  vs.  depth  plots  for 
the  T-7/Priestly-Brune  and  A-10  models.  Assumed  Q  structure  at 
left;  a  (dimensionless)  is  Poisson's  ratio. 


regional  short-period  seismograms.  The  technique  appears 
promising  in  advancing  knowledge  of  wave  propagation  and  source 
identification  at  regional  distance  ranges. 


ACKNOWLEDGMENTS 

We  especially  thank  Rainer  Kind  for  making  available  to  us  the 
atodified  reflectivity  method  computer  program  that  we  have 
adapted  for  our  analyses.  Paul  A.  Johnson  was  responsible  for 
the  computer  runs.  We  thank  Terry  C.  Wallace  and  Mike  Shore  for 
discussions  and  for  reviewing  the  manuscript.  The  calculational 
and  data  reduction  efforts  for  this  research  were  supported  by 
Che  U.S.  Department  of  Energy  and  partially  by  ONR  Earth  Physics 
Program  grant  N00014-75-C-0972  to  L.W.B.  The  Snake  River  Plains 
data  were  collected  during  research  partially  funded  by  the  U.S. 
National  Science  Foundation  (grant  EAR-77-23707  to  the  Univer¬ 
sity  of  Utah  and  EAR-77-23357  to  Purdue  University)  and  by  the 
U.S.G.S.  Geothermal  Research  Program  grant  14-08-000 l -G- 5 32  to 
Purdue . 


t 


66 


466  K.H.  OLSEN  AND  L  W  BRAILt 

REFERENCES 

l.  A.  I.  Ruzaikin,  I.  L.  Nersesov,  V.  I.  Khalturin,  and  P. 

Molnar,  J.  Geophys.  Res.  82,  pp.  307-316,  1977.' 

2  L.  Knopoff,  F.  Schwab,  K.  Nakanishi,  and  F.  Chang,  Geophys. 
J.  R.  Ascr.  Soc.  39,  pp.  41-70,  1974. 

3.  V.  Z.  Ryaboi,  Izv.  (Bull.)  Acad.  Sci.  USSR,  Geophys.  Ser., 
AGU  Trans.  3,  pp.  177-184,  1966. 

4.  R.  P.  Masse,  Bull.  Seism.  Soc.  Am.  63,  pp.  911-935,  1973. 

5.  A.  Hirn,  L.  Steinmetz,  R.  Kind,  and  K..  Fuchs,  Z.  Geophys. 
39,  pp.  363-384,  1973. 

6.  R.  Kind,  J.  Geophys.  40,  pp.  189-202,  1974. 

7.  C.  B.  Archambeau,  E.  A.  Flinn,  and  D.  G.  Lambert,  J. 

Geophys.  Res.  74,  pp.  5825-5865,  1969. 

8.  L.  J.  Burdick  and  D.  V.  Helmberger,  J.  Geophys.  Res.  83, 

pp.  1699-1712,  1978. 

9.  J.  E.  York  and  D.  V.  Helmberger,  J.  Geophys.  Res.  78,  pp. 
1883-1886,  1973. 

10.  M.  Cara,  Geophys  J.  R.  Astr.  Soc.  57,  pp.  649-670,  1979. 

11.  B.  A.  Romanowicz  and  M.  Cara,  Geophys.  Res.  Lett.  7,  pp. 

417-420,  1980. 

12.  A.  L.  Hales,  Earth  and  Planet.  Sci.  Lett.  7,  pp.  44-46, 

1969. 

13.  D.  P.  Hill,  Geol.  Soc.  Aa.  Bull.  83,  pp.  1639-1648,  1972. 

’4.  K.  Fuchs  and  G.  MUller,  Geophys.  J.  R.  Astr.  Soc.  23,  pp. 

417-433,  1971. 

15.  R.  Kind,  J.  Geophys.  44,  pp.  603-612,  1978. 

16.  K.  H.  Olsen,  L.  W.  Braile,  and  P.  A.  Johnson,  Geophys.  Res. 
Lett.  7,  pp.  1029-1032,  1980. 

17.  G.  MUller,  J.  Geophys.  42,  pp.  429-436,  1977. 

18.  K.  Priestly  and  J.  Brune ,  J.  Geophys.  Res.  83,  pp. 
2265-2272. 

19.  K.  Fuchs,  Tectono physics  56,  pp.  1-15,  1979. 

20.  L.  Knopoff,  Rev.  Geophys.  2,  pp.  625-660,  1964. 

/ 


GEOPHYSICAL  RESEARCH  LETTERS,  VOL.  7,  NO.  12,  PAGES  1029-1032,  DECEMBER  1980 


67 


SEISMIC  VELOCITY  AND  Q-STRUCTURE  OF  THE  UPPER  MANTLE  LID 
AND  LOW  VELOCITY  ZONE  FOR  THE  EASTERN  GREAT  BASIN 

K.  H.  Olsen,1  L.  W.  Braile,2  and  P.  A.  Johnson1 

1Geoseiences  Division,  Los  Alamos  Scientific  Laboratory,  Los  Alamos,  New  Mexico  87S45 
^Department  of  Geosciences,  Purdue  University,  West  Lafayette,  Indiana  47907 


Abstract.  A  100-km-long  record  section  of  NTS 
explosions  recorded  in  the  eastern  Snake  River  Plains 
(7°<  a  <  8°)  shows  the  cusp  of  critical  refractions 
from  the  steepened  P  velocity  gradient  at  the  bottom 
of  the  upper  mantle  LVZ.  Synthetic  seismograms 
calculated  with  a  modified  reflectivity  program  have 
heen  used  to  derive  a  regional  velocity  model  of  the 
upper  mantle  beneath  the  eastern  Great  Basin.  The 
model  suggests  that  observed  very  weak  P»  arrivals 
are  due  to  a  slight  negative  velocity  gradient  below  the 
Moho  and  that  no  high  velocity  mantle  lid  exists  in  this 
region. 

Introduction 

The  seismic  velocity  versus  depth  structure  of  the 
upper  mantle  and  lower  crust  beneath  tectonically 
active  areas  of  the  western  United  States  has  been 
studied  extensively  for  nearly  20  years.  This  has  been 
possible  because  Nevada  Test  Site  (NTS)  underground 
explosions  and  western  U.S.  and  Mexican  earthquakes 
provide  frequent  seismic  sources  in  an  area  well 
covered  by  seismograph  stations.  Compressional 
velocity  distributions  have  mainly  been  determined  by 
integrating  the  slope  of  the  travel  time  curve,  dT/di, 
using  the  Herglotz-Wiechert  method.  The  required 
travel  time  (T)  versus  distance  (a)  date  have  been 
analysed  from  short-period  recordings  obtained  along 
long-range  profiles  (Archambeau  et  al.,  1969;  Masse  et 
al.,  1972)  and/or  from  apparent  velocities  measured 
directly  across  large  seismic  arrays  (Johnson,  1967). 
Recently,  availablility  of  high  speed  computers  and 
development  of  sophisticated  synthetic  seismogram 
modeling  techniques  make  it  practical  to  fit  the  travel 
time  and  amplitude  date  by  a  trial  and  error  procedure 
(Burdick  and  Helmberger,  1978;  Wiggins  and 
Helmberger,  1973).  The  important  advantage  of  the 
synthetic  seismogram  method  is  that  it  makes  optimum 
use  of  amplitude  date  and  detailed  waveform  fitting  to 
derive  P  velocity  structure. 

Many  compressional  and\shear  wave  studies  show 
that  a  major  feature  of  the  mantle  structure  beneath 
the  western  U.S.  is  a  low  velocity  zone  (LVZ)  in  the 
depth  range  between  60  and  300  km.  It  is  well  known 
that  significant  lateral  variations  in  LVZ  properties 
(thickness,  depth,  values  of  minimum  S  and  P  velocities, 
presence  or  absence  of  a  lithospheric  "lid,"  etc.)  occur 
over  distances  of  several  hundred  kilometers  and 
perhaps  to  even  finer  scales  (Burdick  and  Helmberger, 
1978;  York  and  Helmberger,  1973;  Romanowicz  and 
Cara,  1980).  On  the  other  hand,  Burdick  and 
Helmberger  (1978)  suggest  mantle  structure  deeper 
than  about  300  km  is  more  uniform  over  a  global  scale 
and  therefore  amenable  to  modeling  using  widely 
spaced  sources  and  seismograph  stations  if  emphasis  is 
placed  on  long  period  body  wave  arrivals  at  distances 
beyond  10°.  Here  we  report  on  a  record  section  of 
NTS  explosions  taken  with  matched  short-period 
instruments  having  a  sufficiently  small  station  spacing 
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(8  km)  and  yet  long  enough  MOO  km)  to  identify  at  least 
three  distinct  (T,  a)  branches  for  P  waves  whose 
raypaths  -  bottom  in  the  uppermost  mantle  beneath  a 
small  area  in  east-central  Nevada.  Modeling  of  the 
arrival  times,  amplitudes,  and  waveforms  usii«  a 
reflectivity  method  synthetic  seismogram  program 
(Kind,  1978;  Fuchs  and  Muller,  1971)  enables  us  to 
perturb  the  generic  western  U.S.  models  into  a 
crust-upper  mantle  model  which  gives  fine  details  of 
the  LVZ  transition  in  this  region. 

}  Observations 

Our  observations  are  recordings  of  two  NTS  nuclear 
explosions  obtained  while  our  equipment  was  deployed 
in  eastern  Idaho  during  the  Yellowstone-Snake  River 
Plains  (Y-SRP)  cooperative  seismic  profiling  experi¬ 
ment  (Braile  et  al.,  1979).  Twelve  special  high-explosive 
shots  plus  blasts  at  two  quarries  were  used  as  sources 
for  crustal  profiles  in  eastern  Idaho  and  Yellowstone 
Park.  Figure  la  shows  the  area  of  the  Y-SRP 
experiment;  Figure  lb  indicates  those  stations  that 
were  recorded  on  an  approximate  radial  line  to  two  NTS 
explosions  on  September  27, 1978  (Table  1).  Because 
RUMMY  and  DRAUGHTS  explosion  sites  were  within 
3  km  of  each  other,  our  observed  record  sections  are 
nearly  identical  except  DRAUGHTS  amplitudes  are 
about  1/4  RUMMY  amplitudes.  We  discuss  only  the 
better  signal -to- noise  RUMMY  seismograms. 

Instrumentation  consisted  of  13  vertical  component 
short-period  (1  Hz  natural  frequency)  seismometers. 
Ten  of  these  were  telemetered  to  a  centrally  located 
site  and  recorded  on  analog  magnetic  tape;  the  three 
southernmost  instruments  were  recorded  on  portable 
smoked  paper  units  and  FM  tape  recorders.  All  records 
were  digitized  at  100  samples  per  second  and  filtered 
(0-3  Hz)  for  this  analysis. 

The  reduced-time,  true  relative  amplitude  record 
section  for  the  RUMMY  explosion  is  displayed  in 
Figure  2.  Three  separate  compressional  phases  within 
the  first  four  seconds  are  marked  on  Figure  2a;  our 
reasoning  in  so  identifying  these  arrivals  is  as  follows: 

(1)  The  very  first  arrivals  with  an  apparent  velocity 
of  7.8-7. 9  km/s  are  so  weak  that  they  could  easily  be 
missed  on  initial  inspection.  From  the  Y-SRP  refrac¬ 
tion  date,  we  determined  that  the  M-discontinuity  is 
40  km  below  these  stations  and  the  mantle  Pn 
velocity  is  close  to  7.9  km/s.  An  enlarged  view  of  the 
first  12  seconds  is  shown  in  Figure  2b  where  the 
consistency  of  the  Pn  arrivals  across  the  spread  is 
more  apparent.  These  Snake  River  Plains  seismic 
stations  had  quite  low  background  noise  so  the 
implication  is  that  a  true  headwave  P„  arrival  will 
rarely  be  seen  at  distances  beyond  600km  in  the 
western  U.S.,  except  from  events  of  mp  >6.  In  these 
SRP  seismograms  the  ratio  of.  the  amplitudes  of  the 
Pn  arrivals  to  those  of  the  P  phase  is  smaller  than 
0.005.  (The  F  energy  arrives  at  reduced  times  greater 
than  32  seconds  so  is  not  shown  in  Figure  2a).  Other 
investigators  (e.g.,  Hill,  1972,  1973)  have  commented 
that  Pn  energy  at  these  distances  is  probably  very 
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Figure  1.  (a)  Location  map  of  the  western  UA  with 
relative  positions  of  the  Nevada  Teat  Site  and  the 
Y-SRP  refraction  line,  (b)  Enlargement  showing  seismic 
station  in  Idaho  used  for  the  September  27 
observations.  Asterisk  denotes  the  approximate  area  in 
the  Great  Basin  for  mantle  ray  turning  points  from  NTS 
explosions. 


weak  and  that  care  must  be  taken  when  attempting  to 
extend  the  Pn  branch  during  long  range  refraction 
profiling. 

(2)  Following  Pn  by  about  three  secontb  are 
stronger  arrivals  also  having  apparent  velocities  close 
to  8  km/s.  A  striking  feature  of  Figure  2a  is  that 
beyond  780  km  the  amplitude  of  the  second  arriving 
phase  increases  rapidly  with  distance  and  the  dominant 
frequency  is  noticably  lower  (■'•0.8  He)  than  the 
frequencies  for  A<780km  and  for  the  Pn  phase  (both 
■'•1.6  Hz).  This  qualitative  observation  strongly  suggests 
that  the  rapid  increase  of  the  low  frequency  phase  for 
A<  780  km  is  a  manifestation  of  a  critical  distance 
effect  and  that  the  high  frequency  energy  has  been 
attenuated  alor*  the  travel  path.  The  obvious  place  for 
thfa  to  occur  is  during  the  two-way  transit  of  energy 
through  the  mantle  LVZ  (which  also  has  a  high 
anelasticity,  i.e.,  low  Q).  These  critical  refractions  are 
shown  schematically  in  the  ray  diagram  of  Figure  3;  we 
follow  the  convention  of  Arc  ham  beau  et  al.  (1968)  in 
labeling  this  cusp  phase  Pi- 

(3)  The  higher  frequency  second  arrivals  for 
A<  780  tan  we  attribute  to  large  angle  reflections  from 
an  interface  lying  mainly  above  the  LVZ.  These 
reflections  are  overtaken  and  overwhelmed  by  Pi  for 
a<780.  Because  of  the  apparent  velocity  near  8  km/s, 
the  high  frequency  content,  and  travel  time  just  longer 
than  Pn,  the  synthetic  seismogram  modeling  dis¬ 
cussed  below  suggests  this  reflection  occurs  at  the  base 
of  the  mantle  lid  and  hence  our  notation  of  Pu<j. 

Modeling 

Our  technique  in  modeling  the  record  section  was  to 
first  use  a  fast  asymptotic  ray  theory  computer 
program  (Cerveny,  1879)  to  fit  travel  times  and 
approximate  amplitudes.  For  more  exact  modeling  we 
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then  used  a  reflectivity  method  program  developed  by 
Kind  (1978),  which  properly  accounts  for  the  effects  of 
a  buried  source  and  thus  allows  computation  of 
complete  seismograms.  One  advantage  of  the 
reflectivity  method  over  Cagniard-de  Hoop  techniques 
(Helmberger,  1973;  Helmberger  and  Burdick,  1979)  is 
that  Q  values  can  be  individually  assigned  to  each 
model  layer  rather  than  distributed  over  the  entire  path 
as  part  of  a  linear  operator. 

Burdick  and  Helmberger's  (1978)  T-7  model  was 
adopted  as  the  starting  model  for  compressions!  wave 
velocities  in  the  crust  and  mantle.  The  generic  T-7 
model -was  constructed  mainly  from  long  period  data  to 
the  NW  and  SE  of  NTS— with  emphasis  on  velocity 
structure  below  200  km  (arrivals  for  A  >  10°).  Since 
our  observations  are  in  the  range  7°<  A  <8°,  we 
perturbed  the  initial  model  only  at  depths  above 
250  km.  The  T-7  model  has  a  Pn  velocity  of 
7.95  km/s  with  a  positive  gradient  below  the 
M-diseontinuity  to  8.05  km/sec  at  the  bottom  of  the  lid 
at  65  km.  A  substantial  LVZ  for  P-velocities  is 
ineludedfbelow  65  km  (Figure  5). 

We  did  not  calculate  synthetics  for  shear  wave 
phases  but  were  required  to  include  a  realistic 
S- velocity  structure  because  shear  wave  velocity 
contrasts  can  have  a  major  influence  on  P-wave 
reflection  coefficients— especially  for  large  angles  of 
incidence.  Priestly  and  Brune  (1978)  used  dispersion  of 
fundamental  mode  Rayleigh  and  Love  waves  to  derive  a 
shear  velocity  model  in  the  eastern  Great  Basin  very 
dose  to  the  area  of  the  mantle  turning  points  of  this 
study.  The  combined  P-  and  S-velocity  model,  T-7/PB, 
is  shown  in  Figures  along  with  Poisson's  ratio  (o) 
calculated  from  the  tabulated  velocities. 

The  modified  reflectivity  synthetics  for  the  T-7/PB 
velocity  model  are  shown  in  Figure  4a  tor  an  extended 
range  from  600  km  to  960  km.  The  Pi  phase  can  be 
seen  only  for  distances  beyond  840  km  and  reduced 
times  greater  than  13  seconds.  In  order  to  bring  the 
synthetic  Pi  phase  into  agreement  with  the  observed 
arrival  times  and  to  shift  the  cusp  from  ■'•820  km  back 
to  ■'•780  km  it  was  necessary  to  bring  the  gradient  at  the 
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Figure  2.  (a)  True  relative  amplitude  record  section  of 
P-wave  arrivals  from  RUMMY.  (b)Same  as  (a)  with 
increased  amplitudes  to  show  weak  Pn  phase. 
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Figure  3.  Schematic  ray  diagram  showing  the  head 
wave  phase  Pn>  the  P^  phase  critically  refracted 
from  the  gradient  near  the  bottom  of  the  LVZ,  and  the 
I  Pjjd  phase  reflected  from  the  base  of  the  mantle  lid. 


bottom  of  the  LVZ  to  shallower  depths  and  to  make 
slight  adjustments  to  its  curvature.  Also  note  from 
Figure  4a  that  the  T-7/PB  model  gives  too  large 
amplitudes  for  P„  arrivals,  and  the  Pud  phase  from 
the  discontinuity  at  65  km  depth  arrives  about  2.5 
seconds  too  early  so  is  superimposed  on  the  Pn  phase 
throughout  much  of  the  600-960  km  range.  The  result 
of  perturbing  the  T-7/PB  velocity  model  to  better 
match  the  observations  of  Figure  2  is  model  A-10  shown 
in  Figure  5.  The  A-10  synthetics  are  compared  with  the 
original  T-7/PB  model  in  Figure  4b  and  with^the  more 
limited  distance  range  observations  in  Figure  6. 

Discussion 

We  can  summarise  the  nature  and  reasons  for  the 
various  model  perturbations  as  follows: 

(1)  The  steepened  positive  P-velocity  gradient  in  the 
lower  part  of  the  LVZ  has  been  raised  in  order  to  fit 
the  arrival  times  and  cusp  distance  of  the  P^  phase. 

(2)  bi  order  to  match  the  observed  weak  Pn,  a 
slight  negative  gradient  just  below  the  M-d  is  continuity 
is  required  instead  of  the  positive  gradient  of  the 
generic  T-7  model.  In  fact,  our  A-10  P-veloelty  model 
suggests  that,  beneath  this  part  of  the  Great  Basin,  the 
LVZ  may  be  in  contact  with  the  crust  at  the 
M-d  is  continuity.  Similar  indications  of  the  absence  of  a 
high  velocity  lid  in  parts  of  the  western  U.S.  have  been 
cited  by  Archambeau  et  al.  (1969)  (especially  for  their 
SHOAL-FALLON  SE  profile). 

(3)  The  large  negative  discontinuity  at  a  depth  of 
65  km  present  in  both  the  T-7  P- velocities  and  the 
Pries  uy-Bneie  model  appears  to  be  too  shallow  to 
properly  match  the  observed  Pu<r^n  trtv*l  time 
delay.  The  observed  delay  is  better  reproduced  if  the 
discontinuity  is  at  about  100  km  depth. 

(4)  The  calculated  amplitudes  of  the  Pjjd  re¬ 
flection  are  much  too  large  if  the  T-7/PB  velocity 


Figure  4.  Synthetic  seismograms  (Z-component)  of 
early  compreasional  phase  arrivals  from  (a)  the  generic 
T-7/PB  model  and  (b)  the  A-10  mantle  model  (Figure  5) 
which  match  the  observations  in  the  720  to  820-km 
range.  Travel  time  curves  calculated  from  the  Cerveny 
program.  Amplitude  multiplied  by  distance  for 
convenient  plotting. 
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Figure  5.  P-velocity  (a)  and  S-velocity  (B)  versus  depth 
plots  for  the  T-7/PB  and  A-10  models.  Assumed 
Q-structure  for  both  models  shown  at  left,  o  is 
Poisson's  ratio. 


contrasts  (-U3  km/s  for  P,-0.39  km/s  for  S)  are  used  at 
a  depth  of  MOO  km.  In  fact,  by  including  the  effect  of 
S-velocity  contrasts  on  the  P-wave  reflection  co¬ 
efficients,  we  can  match  observed  Pyd  amplitudes  by 
keeping  the  P-velocity  contrast  at  zero  and  relyiiv 
entirely  on  an  S  contrast  of  -0.15  km/s  to  produce  the 
effect  on  these  wide  angle  reflections.  In  moviig  the 
bottom  of  the  S-velocity  lid  from  a  depth  of  65  km 
proposed  by  Priestly  and  Brune  to  the  M00  km  required 
in  our  A-10  model,  we  introduced  a  negative  gradient  in 
the  S- velocities  between  these  two  depths;  a  similar 
negative  gradient  can  be  seen  in  the  higher  mode 
inversions  by  Cara  (1979),  but  these  are  not  (dotted  in 
Figure  5.  HOI  (1972)  and  Hales  (1969)  previously  re¬ 
ported  arrivals  following  Pn  by  two  to  three  seconds 
at  distances  of  ■vSOOkm  in  sections  from  long  range 
refraction  experiments  in  the  Columbia  Plateau 
(EDZOE  experiment)  and  the  central  U.S.  (EARLY 
RISE),  respectively.  Their  interpretations— using  travel 
time  information  only— suggest  a  thin  (M0  km),  sharp, 
but  high  velocity  (8.0  to  8.4  km/s)  lid  at  depths  of 
90-100  km  is  present  in  those  regions.  Our  Great  Basin 
data  agrees  in  placing  a  discontinuity  (which  is  perhaps 
the  "boundary"  between  the  lithosphere  and  the 
athenosphere)  at  M00  km  but  our  P-velocity  contrast 
cannot  be  as  pronounced  as  those  implied  by  Hill  and 
Hales  and  still  give  rise  to  the  comparitively  weak 
amplitudes  that  we  observe  in  the  750  km  raiye. 

(5)  the  Q-structure  (for  P-waves)  used  for  the  A-10 
model  was  a  generalization  of  proposed  values  that 
have  appeared  in  recent  literature,  the  most  important 
segment  is  the  low  value  centered  in  the  LVZ.  A  Qq 
value  in  the  range  between  50  and  100  appears  to 
adequately  attenuate  the  higher  frequency  components 


Figure  6.  Comparison  of  observed  (a)  and  synthetic  (b) 
seismogram  record  sections  for  the  720  to  820-km 
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of  the  Pi  phase.  Q„  »  100  in  the  LVZ  does  not 
attenuate  Pi  enough,  whereas  Qa  ~  25  completely 
obliterates  the  Pi  phase  in  the  synthetics.  Our  value 
of  50<Qo<100  is  of  the  same  order  as  that  deduced  by 
Helmberger  (1973)  from  Cagniard-de  Hoop  techniques. 

(6)  One  possible  shortcoming  of  the  A-10  model  is 
the  failure  to  reproduce  details  of  the  oscillations  of 
the  observed  Pj  phase.  We  do  not  believe  this  to  be  a 
result  of  an  inadequately  detailed  source  spectrum, 
since  a  comparison  of  the  explosion  source  spectrum 
algorithm  used  in  the  modified  reflectivity  code  is 
reasonably  represented  by  the  source  spectrum  plus 
instrument  response  function  calculated  from  known 
physical  parameters  of  these  explosions  (Mueller  and 
Murphy,  1971).  Archambeau  et  al.  (1969)  observed 
compressional  wave  energy  spread  out  in  long,  rather 
complicated  oscillatory  wave  trains  near  caustics  and 
attributed  this  to  interference  between  refracted  and 
reflected  components  near  the  cusp.  Our  model  layer 
thicknesses  (^5  km)  in  the  region  of  the  lower  depths  of 
the  LVZ  (120-150  km)  are  of  the  same  order  as  the 
wavelengths  ( MO  km)  of  the  dominant  short  period 
energy.  Thus,  we  believe  the  oscillatory  Pi  trains 
may  be  due  to  small  details  of  fine  structure  in  the 
transition  zone  which  we  have  not  yet  attempted  to 
model  at  the  required  resolution. 

In  conclusion,  relatively  minor  adjustments  in  the 
T-7/PB  model  for  the  western  U.S.  yield  an  uppermost 
mantle  structure  that  reproduces  in  detail  the  upper 
mantle  arrivals  and  very  weak  Pn  observed  in  the 
Great  Basin. 
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Ttx  short-period  seismic  phases  knows  as  P  and  Lg  are  ofte^  recorded  at  distances  of  200-1000  km  on  loog-ragt 
rafraetion  profiles  and  are  usually  the  larpett-amplitudc  feature*  on  record  uctiotu  for  this  nape.  P  and  Lg 

propagate  a*  multiply  reflectad  rrsnprniirmil  and  shear  waves  in  a  crustal  waveguide  whose  principal  boundaries  an 
the  Mobo  and  the  tree  surface  Equivalently,  they  can  be  interpreted  as  the  interference  pattern  produced  by  a 
superposition  of  higher-mode  t.  SV  tad  SH  wave  propgating  is  a  leaky  waveguide  For  comprasrional  waves,  the 
waveguide  efficiency  is  a  strong  function  of  frequency  and  depends  on  the  pretence  or  abeenee  of  low  velocity  layers 
within  a  law  kilometers  of  the  surface,  such  as  deep  sedimentary  sections  commonly  found  is  active  -"-w  areas.  Such 
low-velocity  surface  layers  create  constructive  interference  effects  far  upcoming  P  waves  inodent  at  near  gracing  angles 
at  the  baa  surface  and  land  to  efficient  P  propagation.  Several  good  examples  of  strong  P  phases  can  be  found  on 
long-range  refraction  profiles  for  the  tactooicaDy  active  centra  United  States;  the  550  km  profile  eastward  from 
SHOAL  to  Delta.  UT  is  analyzed  hoe.  We  have  used  a  modified  reflectivity-method  computer  program  to  model 
crustal  phases  for  the  SHOAL- Della  profile.  The  reflectivity  technique  accounts  for  all  body  and  surface  waves 
contributing  to  the  short-period  seismogram*.  It  is  found  that  the  synthetic  waveforms  realistically  model  the  obmrved  P 
characteristic*.  In  this  ease,  the  decay  of  P  amplitudes  with  distance  appears  to  be  by  nrfaee-reflectian 

t— t-g-  from  the  svaveguide  rather  than  by  inelastic  attenuation  due  to  Q  of  crustal  rocks. 


1.  Introduction 

Until  recently,  the  velocity  structure  of  the 
Earth’s  crust  was  inferred  mainly  from  travel-time 
studies  of  easily  correlatable  seismic  phases  (usu¬ 
ally  first  arrivals)  recorded  in  refraction  profiling 
experiments.  Amplitude  and  waveform  informa¬ 
tion  oo  scunic-record  sections  was  used  mainly  in 
a  qualitative  way  in  phase  identifications  and  as 
an  indication  of  velocity  gradients,  transition  zones, 
and  attenuation  properties.  Since  1970,  the  in¬ 
creased  availability  of  high-speed  computers  and 
the  development  of  sophisticated  synthetic-seismo¬ 
gram  modeling  techniques  has  made  it  possible  to 
use  amplitude  and  waveform  data  optimally  to 
derive  detailed  velocity  and  Q  structure  (MQDer 
and  Fuchs,  1976;  Braile,  1977).  In  particular,  we 
employ  here  a  modified  reflectivity-method  pro¬ 
gram  developed  by  Kind  (1978)  that  is  capable  of 


computing  “complete**  seismograms,  all 

body  and  surface  waves.  Our  objective  is  to  in¬ 
vestigate  the  possible  utility  of  modeling  rather 
complex  short-period  crustal  phases  (e.g,  P  and 
Lg)  as  an  aid  to  the  interpretation  of  crustal 
properties  and  structures.  In  this  paper  we  attempt 
to  describe  the  more  important  characteristics  of 
the  P  and  Lg  phases,  investigate  conditions  for  the 
generation  and  propagation  of  these  phases,  and 
determine  which  properties  of  the  Earth  structure 
influence  this  propagation.  Because  of  the  com¬ 
plexity  and  uncertainties  about  the  theoretical  de¬ 
tails  of  the  generation  and  propagation  of  these 
crustal  phases,  they  have  heretofore  been  little 
used  in  the  interpretation  of  seismic  refraction 
profiles  —  even  though  they  are  often  the  largest- 
amplitude  features  on  record  sections.  We  show 
one  aspect  of  how  modeling  crustal  phases  with 
complete  synthetic  seismograms  can  improve  our 
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knowledge  of  Earth  structure.  Because  relatively 
few  complete  seismic-record  sections  for  regional 
distances  (200-2000  km)  showing  P  and  Lg  phases 
are  available  from  observed  data,  we  rdy  heavily 
on  the  analysis  of  P  and  Lg  waveforms  by  means 
of  synthetic  seismograms  calculated  using  the 
modified  reflectivity  method  for  a  variety  of  Earth 
models. 

For  clarity  in  both  typography  and  context,  we 
shall  henceforth  write  “P-bar"  in  place  of  the  P 
notation. 

2Jtegkmal  Crustal  Phases  (P-bar,  Lg) 

The  P-bar  and  Lg  phases  are  usually  weD  re¬ 
corded  on  short-period  instruments  at  regional 
distance  ranges  (100-2000  ltm).  They  are  .com¬ 
monly  seen  as  broad  (A/>20  s),  complex  wave 
trains  having  group  velocities  of'*-  6.0  km  s~ 1  for 


P-bar  and  -  3J  km  s~‘  for  Lg  (Olsen  and  Braile. 
1981). 

Both  P-bar  and  Lg  phases  propagate  efficiently 
in  the  continental  crust  and  are  often  the  largest- 
amplitude  arrivals  on  short-period  seismograms  at 
regional  distances  (Fig.  1).  In  addition  to  the  com¬ 
plex  and  long-duration  character  of  both  P-bar 
and  Lg  phases,  these  arrivals  also  display  little 
coherence  in  waveform  or  even  in  the  envelope  of 
the  complex  wave  packet  over  relatively  short  dis¬ 
tances.  These  characteristics  of  P-bar  and  Lg  wave 
propagation  can  be  modeled  qualitatively  using 
synthetic-seismogram  techniques.  An  example  of  a 
vertical-component  synthetic-seismogram  calcula¬ 
tion  resulting  from  a  crustal  model  appropriate  to 
the  western  United  Sutes  or  other  tectonic  areas  is 
sbfwn  in  Fig.  2.  Both  P-bar  and  Lg  phases  are 
prominent  arrivals  on  the  record  sections,  having 
durations  of  tens  of  seconds  and  relatively  large 
amplitudes  in  the  short-period  range.  Unlike  the 


T-X/7  fs) 


Fig.  I.  Veaal  eepo— I  record  roctioo  recorded  vanward  from  the  ScbaffarviDe  ibotpoiat  n  the  Superior  Province  of  the 
Canadian  ihtaid  ahoariag  exemplar  of  P-bar  cad  Lg  enuul  pharos  The  broad  P-bar  and  Lg  phaaoe  approximately  folio*  the 
average  velocity  baa*  6-2  ka  »-1  aad  3.6  ka  tape cthoty.  for  ranger  X  beyond  —  ISO  km.  Ampbtuda  have  bam  muhtpbed  by 

dtiuw  for  convenient  reeling  Adapted  from  Berry  aad  Ftacb*  (1973). 
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primary  waveforms,  such  as  Pn  and  So,  which 
display  waveform  coherence  over  long  distances, 
the  P-bar  and  Lg  waveform  characteristics  change 
considerably  with  distance,  and  phase  correlations 
are  usually  impossible.  However,  correlation  of  the 
energy  represented  by  the  group-velocity  envelope 
is  usually  possible.  An  enlargement  of  one  of  the 
seismograms  for  the  record  section  shown  in  Fig.  2 


is  illustrated  in  Fig.  3,  indicating  the  waveform 
characteristics  of  the  complex  P-bar  and  Lg 
as  a  function  of  reduced  time  as  well  as  their 
average  or  group  velocity.  It  is  seen  for  this  model 
that  the  P-bar  phase  has  average  velocities  in  the 
range  of  6  km  s-1  and  the  Lg  phase  has  group 
velocities  near  3.6  km  s~'. 

Observations  and  synthetic-seismogram  asaly- 
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Pig.  3.  EaiergwawM  of  Uw  Jf  *900  km  synthetic  aremmgram  from  Fig.  2  showing  more  details  of  individual  phases. 
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Fig.  4.  VenicaJ-compcocni  synthetic  uiimogum  record  MCtioo  mlmUtod  usin|  the  modified  reflectivity  method  for  doeely  spaced 
(10  km)  Mumographs  u  disteww  between  MO  and  950  km  is  explosion  tourer.  Now  how  rapidly  the  waveform  oohrrroor  verier 
with  distance  m  the  P-ber  window  (reduced  ttmet  40-90 1).  The  P- velocity  crustal  model  is  lhosre  at  the  left. 


sis  presented  here  and  by  Bouchon  (1982)  suggest 
that  both  P-bar  and  Lg  propagate  as  multiply 
reflected  compressions!  and  shear  waves  within 
the  crust,  or,  equivalently,  they  are  the  interference 
patterns  produced  by  the  superposition  of  a  large 
number  of  leaky  P-  and  S-modes  respectively 
(Haskell.  1966).  We  have  performed  a  parameter 
study  of  short-period  propagation  from  explosion 
sources  in  continental  crustal  models,  which  shows 
that  the  modified  reflectivity  technique  can  repro¬ 
duce  many  of  the  observed  characteristics  of  these 
phases,  such  as  envelope  modulation,  rapidly  vary¬ 
ing  lateral-waveform  coherence,  etc.  (Olsen  w  -1 
’9813 

(  fat  synthetic  modeling  of  Lg  using  shallow 
Explosion  sources,  the  SV  energy  in  the  Lg  group 


; 

window  arises  from  P  to  S  conversion.  Once  con¬ 
verted,  SV  energy  is  very  efficiently  confined  in  a 
waveguide  whose  principal  boundaries  are  discon¬ 
tinuities  or  sleep  velocity  gradient*  at  tb.  Moho 
and  at  the  free  surface.  Although  we  have  not  here 
modeled  the  SH  component  of  motion  of  the  Lg 
phase,  we  expect  that  it  also  can  be  represented  as 
a  guided  wave  with  the  near  surface  and  Moho  as 
waveguide  boundaries,  as  Bouchon  (1982)  has  also 
suggested.  The  total  channeled  energy  is  relatively 
insensitive  to  fine-scale  detail  (a  few  km)  of  verti¬ 
cal  velocity  gradients  near  the  waveguide 
boundaries,  but  the  waveform  coherence  at  surface 
seismometers  can  still  vary  appreciably  over  dis¬ 
tances  of  — 10  km  (Figs.  4  and  5).  Undercritical 
and  wide-angle  reflections  at  the  waveguide 


1NTER-STRT1QN  aXWCE  FO*  **(L6I 


m 


Fig.  5.  Synthetic  record  rection  Qhmretmg  rapid  lawn]  chutes  in  waveform  coharcocr  in  the  L*  window  (reduced  tunes  130-190  s) 

for  the  same  model  and  rente*  as  m  Fig  4. 


'Z 


V  * 


boundaries  are  nearly  total,  so  attenuation  and  the 
frequency  spectra  depend  almost  entirely  on  Q  of 
crustal  materials.  Thus  Lg  is  a  good  measure  of  Q. 

There  is  frequently  a  confusion  or  imprecision 
in  the  literature  concerning  the  nomenclature  for 
the  P-bar  phase,  with  the  term  Pg  being  applied  to 
any  phase  following  Pn  and  having  an  apparent 
velocity  of  —  6.0  km  s~'.  In  the  western  United 
Sutes  where  the  complex  P-bar  phase  can  be  very 
prominent  at  distances  of  200-2000  km,  the  Pg 
phase  —  a  bead  wave  traveling  along  the  top  of 
the  “granitic"  basement  —  attenuates  extremely 
rapidly  and  is  usually  unobservable  beyond  —  ISO 
1cm  (Ryall  and  Stuart.  1963).  So,  although  P-bar 
onsets  fall  on  the  extension  of  the  Pg  travel-time 
branch,  the  two  phases  have  quite  different  prop¬ 
agation  and  attenuation  characteristics.  Similarly 
to  Lg.  the  “true"  P-bar  phase  propagates  by  multi¬ 
ple  wide-angle  reflection  of  PxnP,  but  the  effi¬ 
ciency  of  the  upper  waveguide  boundary  is  a  com¬ 
plex  function  of  frequency  and  is  very  sensitive  to 
the  presence  of  low-velocity  layers  at  or  near  the 
surface.  Our  parameter  studies  (e.g.,  Figs.  6  and  7) 
confirm  that  a  low-velocity  surface  layer  —  such 
as  the  thick  sedimentary  section  in  tectonic  areas 
such  as  the  Basin  and  Range  and  Colorado  Plateau 
provinces  of  the  western  United  Sutes  —  is  the 
key  factor  for  efficient  propagation  of  P-bar.  The 
controlling  factor  is  the  nature  of  the  ftee-surface 
reflection  of  upcoming  P-waves  previously  re¬ 
flected  from  the  Moho  and  midcrustal  discontinui¬ 
ties.  la  the  absence  of  low-velocity  surface  layers 
(Lc.,  in  shield  areas)  nearly  all  upward-traveling  P 
energy  is  converted  to  SV  upon  free- surface  reflec¬ 
tion.  Note  the  very  small  values  of  the  P-P 
energy-reflection  coefficient  (R0PPa<  0)  in  Fig. 
6(A);  the  energy  fraction  converted  to  SV  at  the 
free  surface  is  Ron™  1  —  Ron  (HaskeU,  1966). 
Thus,  P-SV  conversion  is  nearly  complete  and 
P-bar  does  not  propagate  efficiently  in  shields. 

On  the  other  hand,  when  a  low-velocity 
near-surface  layer  is  present,  Ron  is  a  strong 
function  of  frequency  (owing  to  constructive  and 
destructive  interference  effects)  but  on  the  avenge 
is  large  (Figs.  6(B)  and  6(C).  A  relatively  large 
amount  of  compressional  energy  is  reflected  back 
into  the  crust,  and  P-bar  propagates  with  mod¬ 
erate  efficiency.  Figure  7  illustrates  the  effect  of  a 
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Fig. «.  Eaergy-ftfWctioo  coefficient  as  a  function  of  angle 
of  incidence  aad  frequency  for  lb*  thm  Basis  and  Range 
crustal  medals  shown  in  Fig  14.  is  the  fraction  of  modem 

F-wava  energy  reflected  ss  P  from  the  free  surface.  Now  that, 
for  the  free  surface,  the  fraction  of  F-wtve  energy  convened  to 
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low-velocity  surface  layer  on  both  P-bar  and  Lg 
propagation.  Figure  7(B)  illustrates  that  as  the 
thickness  of  the  low-velocity  surface  layer  is  in¬ 
creased,  the  efficiency  of  propagation  of  the  P-bar 
phase  increases.  With  no  low- velocity  sedimentary 
layer  (Model  TIOCANS1),  the  P-bar  arrivals  are 
extremely  small  —  having  amplitudes  of  the 
order  as  those  of  primary  arrivals  such  as  Pn  and 
P*.  With  the  presence  of  a  thin  low-velodty 
surfacelayer.  large-amplitude  P-bar  arrivals  are 
present.  Varying  the  thickness  of  the  low-velocity 
surface  layer  does  not  change  thepeak  amplitude 
of  the  P-bar  phase  appreciably,  but  does  signifi¬ 
cantly  affect  the  complex  waveform  characteristics 
of  the  phase.  The  exact  waveform  character  of  the 
P-bar  phase  cannot  be  predicted  at  any  given 
distance,  even  for  these  laterally  homogeneous 
modbls,  because  the  lateral  coherency  (Fig  4)  varies 
rapidly  owing  to  the  constructive  and  destructive 
interference  effects  inherent  in  the  propagation  of 
P-bar.  Barker  et  aL  (1981)  have  emphasized  that 
the  P-bar  and  Lg  phases  are  signficantly  affected 
by  local  geological  structure  (site- response  condi¬ 
tions).  They  have  shown  that  variations  in  ampli¬ 
tude,  by  a  factor  of  as  much  as  10,  and  significant 
differences  in  the  character  and  waveform  of  both 
P-bar  and  Lg  may  be  caused  by  local  geological 
conditions  in  the  vicinity  of  the  recording  seismo¬ 
graph.  However,  these  local  effects  must  be  dis¬ 
tinguished  from  the  waveform  complexity  which  is 
due  to  constructive  and  destructive  interference 
effects  in  the  propagation  of  the  P-bar  and  Lg 
phases.  For  example,  in  the  very  simple  and  later¬ 
ally  homogeneous  models  that  we  have  shown  here 
(Figs.  2,  4,  5,  7(B),  (Q,  the  detailed  characteristics 
of  the  P-bar  and  Lg  phases  are  complicated  be¬ 
cause  of  propagation  effects  even  though  the  pos¬ 
sible  contributions  of  local  geological  conditions 
have  not  been  included. 

Haskell  (1966)  has  discussed  the  attenuation  of 
P-bar  amplitudes  with  distance  due  to  P  to  SV 
conversion  by  multiple  surface  reflections,  and 


CANS!)  low- velocity  turtle*  layers;  (B)  synthetic  teumopmaa 
at  900  km  ranga  showing  c  jges  in  amplitude  and  waveform 
characteristic*  m  the  P-bar  time  windo w,  (Q  same  as  (B) 
except  for  the  Lg  tune  window.  Only  vertical  coopoecnu  are 
ihoem. 


concluded  that  this  leakage  attenuation  usually 
dominates  over  Q-effects  even  in  the  western 
United  Sutes  where  P-bar  is  often  very  strong  on 
short-period  seismograms.  Our  work  on  the 
SHOAL- Delta  record  section  discussed  below  is 
in  agreement  that  P-bar  attenuation  is  a  poor 
measure  of  Q  for  crustal  rocks. 

Shurbet  (1960)  observed  P-bar  phases  from 
Nevada  Test  Site  (NTS)  explosions  at  Lubbock, 
TX  (—1400  km  range)  and  suggested  (Shurbet, 
1969)  that  a  low-velocity  (V-.  —  5.5  km  s“ ')  chan¬ 
nel  at  depths  between  5  and  10  km  would  explain 
several  features  of  the  P-bar  phase  on  the  Lubbock 
seismograms.  Our  reflectivity-method  calculations 
show  that  such  a  low-velocity  crustal  channel  is 
not  uniquely  required  for  P-bar  propagation  — 
low-velocity  layers  at  the  surface  are  sufficient  to 
form  a  good  P-bar  waveguide.  However,  price 
P-bar  is  essentially  an  interference  pattern 
produced  by  layers  less  than  one  seismic  wave¬ 
length  in  thickness  (analogous  to  thin. film  coat¬ 
ings  in  optics),  shallow  low-velocity  channels  would 
affect  the  details  of  the  P-bar  coda.  Thus,  certain 
details  of  P-bar  phases  may  in  some  ««««•*«  lend 
support  to  suggested  low-velocity  channel*  in  the 
crust. 

Strong  P-bar  phases  are  often  found  on  long- 
range  refraction  profiles  for  other  parts  of  the 
world,  as  well  as  for  the  western  United  States. 
Good  examples  are  the  Eschelohe  NW  profile  in 
southern  Germany  (Mueller,  1977),  the  900  km 
Brest  to  Toulon  profile  in  France  (Him  et  aL, 
1973;  Kind,  1974)  and  the  800  km  EDZOE  profile 
in  the  Canadian  Rockies  and  foreland  of  Alberta 
(Menu  et  aL,  1977).  The  Superior  Province 
Canadian  shield  profiles  discussed  by  Berry  and 
Fuchs  (1973)  show  fine  examples  of  Lg  as  well  as 
intermittent  P-bar  phases  (Fig.1).  Peterschmitt 
(quoted  by  Mueller,  1977)  has  analyzed  P-bar  on 
several  European  profiles  and  also  concluded  that 
P-bar  consists  of  a  superposition  of  multiply  re¬ 
flected  PmP  waves,  which  are  recorded  only  if  the 
upper  reflection  point  atthe  surface  lies  in  sedi¬ 
mentary  terrain. 

As  suggested  by  Haskell  (1966)  and  Bouchoa 
(1982),  the  lower  boundary  of  tbe  waveguide  for 
P-bar  and  Lg  propagation  is  the  Moho  (with  possi¬ 
ble  contributions  from  other  discontinuities  within 


the  lower  crust).  Three  observations  provide  evi¬ 
dence  for  the  efficiency  of  P-bar  and  Lg  propaga¬ 
tion  being  due  to  a  waveguide  effect  in  which  the 
lower  boundary  of  the  waveguide  is  the  Moho 
discontinuity.  First,  Haskell  (1966)  has  shown  that 
the  Medio  velocity-transition  zone  (or  discontinu¬ 
ity)  is  an  effective  reflector  for  wide-angle  P  and 
SV  reflection  (large  reflection  coefficients).  Sec¬ 
ondly,  observed  seismic  refraction  and  wide- 
angle-reflection  record  sections  for  both  thjrfH  and 
tectonic  areas  (Figs.  8-10)  indicate  that  P-wave 
reflection  from  the  Moho  (phase  PmP)  not  only  is 
large,  but  multiple  PmP  arrivals  can  also  be  dis¬ 
tinguished  on  tbe  seismic-record  section.  These 
multiple  PmP  reflections  (as  well  as  other,  more 
complicated  paths)  contribute  to  high-amplitude 
arrivals  which  have  a  complicated  and  long-dura¬ 
tion^  wave  character  but  which  travel  within  an 
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Fit. 1  Record  taction  from  tbe  Canadian  thiald  experiment 
(Barry  and  Fuchx,  1973)  iUustratins  bow  multiple  reflectiont  of 
PmP  (PmPPmP)  merge  at  rantat  beyond  400  km  to  term  part 
of  the  P-bar  wave  group. 


Fi(.  9.  V«nic«l  mmpotm  raoord  mctioo  from  NTS  to  Elko,  NV.  ilhutntmg  mcrpng  of  multiple  PmP  at  beyond  300  km  to 

form  pan  of  P-6ar. 


average  velocity  of  —6  1cm  s~'  across  the  record 
sections  at  ranges  of  several  hundred  kilometers  or 
greater.  Thirdly,  our  synthetic-seismogram  model 
studies,  as  shown  above  and  in  the  synthetic  wave¬ 
form*  calculated  for  the  SHOAL-Delta  profile 
discussed  below,  as  well  as  the  results  of  Bouchon 
(1982),  also  indicate  that  the  Moho  discontinuity 
is  an  effective  waveguide  boundary  for  P-bar  and 
Lg  propagation.  Thus,  we  view  the  generation  and 
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propagation  of  the  P-bar  phase  as  illustrated 
schematically  in  Fig.  11.  Multiple  PmP  reflections, 
as  well  as  more-complex  multiples  wihin  the  low- 
velocity  surface  layer  and  possible  P  to  S  conver¬ 
sions  and  multiples  at  the  Conrad  discontinuity, 
contribute  to  wide-angle  reflected  energy  having  a 
complex  interference  waveform  propagating  at  —6 
km  s-'  at  distances  of  several  hundred  kilometers 
from  a  source.  A  low-velocity  surface  layer  is 


Fij.  11.  Schematic  diagram  indicating  how  multiple  reflection*  of  rmF  merge  to  font  the  crustal  r-bar  phase  for  ranges  beyond  ->300 
km.  Except  for  details  of  interference  effects  in  the  Deer-surf  see  low-velocity  lsyert,  the  tune  fenerxl  physical  picture  pertains  to  the 
SV  coenpooent  of  Lf  propagation.  with  PmP  replaced  by  SmS,  etc. 


necessary  for  efficient  propagation  of  the  P>bar 
phase.  If  the  surface  low-velocity  layer  is  not  pre¬ 
sent,  the  majority  of  multiple  PmP  energy  is  even¬ 
tually  convened  to  SV  motion  as  indicated  by  the 
reflection  coefficients  shown  in  Fig.  6  and  the  re¬ 


sults  of  synthetic-seismogram  calculations  il¬ 
lustrated  in  Fig.  7.  We  expect  that  both  the  SV  and 
SH  components  of  Lg  wave  propagation  ere  anal¬ 
ogous  to  those  shown  in  Pig.  11  tor  P-  bur,  except 
that  the  low-velodty  surface  layer  is  not  required 
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in  order  to  produce  an  efficient  waveguide.  The 
free  surface  and  the  Moho  and  other  crustal  dis¬ 
continuities  are  adequate  for  producing  a  wave* 
guide  for  the  Lg  phase. 


3.  Application  to  SHOAL- Delta  profile 

Figure  12  is  a  map  of  the  western  United  Slates 
showing  the  locations  of  several  refraction  profiles 
that  display  examples  of  the  P-bar  phase.  This 
paper  discusses  only  the  interpretation  of  the  P-bar 
phase  along  the  550  km  profile  extending  eastward 
from  the  nuclear  explosion  SHOAL  to  Delta,  UT. 


*r>*. _ 


The  long-range  SHOAL  profile  was  part  of  a 
network  of  refraction-profile  recordings  mainly  in 
California,  Nevada  and  Utah  carried  out  between 
1961  and  1963  by  field  parties  from  the  US. 
Geological  Survey  (ProdehL  1979).  Approximately 
31  different  explosion  and  earthquake  source  sites 
were  used;  there  were  several  explosions  within  the 
Nevada  Test  Site  (shown  as  only  one  source  loca¬ 
tion  in  Fig.  12).  Many  of  the  crustal  models  de¬ 
rived  by  Prodehl  (1979)  from  the  1961-1963  pro¬ 
files  were  obtained  from  200-  300  km  profiles  using 
chemical  explosions  in  the  Pacific  Ocean,  in 
or  in  drill  boles.  Two  of  these  segmented 
shorter-range  reversed  profiles  (NTS- Boise.  ID, 
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Fig.  IX  Locations  of  hoc  loeg-rmnge  eeumic  profiles  (heavy  boat)  in  lb*  western  Uni  tad  Slates  which  exhibit  strong  P-bar  phases 
Asterisks  are  nuclear  explosion  sites;  dots  are  rhsmirel  shotpoiais.  Dashed  lines  are  intersecting  shorter- range  profiles  which  wen 
used  for  detailed  crustal- structure  studies.  Nuclear  shotpoints.  S.  SHOAL;  O,  GASBUOOY;  NTS,  Nevada  Test  Site.  Othr- 
shotpoinu:  O.  Delta,  UT;  Eu,  Eureka.  NV;  EL  Elko,  NV;  B.  Botes,  U>;  H,  Hiko.  NV;  F,  FaOoo.  NV  < earthquake).  Simplified  from 
Prodehl  (1979). 
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and  Fallon,  NV-Delta,  UT)  intersect  near  Eureka, 
NV,  which  is  near  the  midpoint  of  the  SHOAL- 
Delta  long-range  profile. 

The  observed  reduced-time  vertical-component 
record  section  from  the  SHOAL  shotpoint  is  given 
in  Fig.  13(A)  (ProdehL  1979),  where  the  onset 
times  of  the  Pn,  PmP  and  P-bar  phases  have  been 
marked.  For  the  —  200  km  profiles  used  for  de¬ 
tailed  crustal-structure  studies,  ten  —  li  km  (six- 
seismometer)  spreads  were  separated  by  distances 
of  — 10-15  km  (ProdehL  1979).  In  contrast,  as 
shown  in  Fig.  13(A),  the  long-range  SHOAL  pro¬ 
file  had  —  50  km  station  separations,  so  the  inter- 


OBSERVED  UNO  SYNTHETIC  SEISnOGMRS 

SHOWING  ftwa  PHASE  A 


Fig.  13.  Comparison  of  observed  (upper)  and  synthetic  (Model 
SD9)  (b»*r)  seismic  record  Met ioes  (or  the  SHOAL-Dcha 
k»g-raa|c  refraction  profile.  Approximate  onset  tune*  lor  Pn, 
PmP  end  P-bar  are  indicated  by  dashed  lines. 


station  coherence  (or  lack  thereof)  of  the  different 
compression*]  phases  is  not  apparent  Our  syn¬ 
thetic-seismogram  modeling  shows  rapid  changes 
in  interstation  coherence  for  the  P-bar  phase  over 
distance  intervals  sometimes  less  than  10  km.  Such 
rapid  coherence  variation  is  one  of  the  more  dis¬ 
tinctive  characteristics  of  P-bar  and  Lg  inter- 
ference-pattem-type  phases,  as  discussed  above. 
However,  the  P-bar  phases  from  SHOAL  (Fig  13) 
well  illustrate  another  characteristic  of  these  phases 
—  the  long  oscillatory  trains  following  emergent 
onsets.  In  Fig  13(A),  the  P-bar  trains  become 
identifiable  at  ranges  between  300  and  400  km. 

Our  procedure  for  modeling  the  SHOAL- Delta 
seismograms  was  to  accept  the  basic  Basin  and 
Range  crust/upper-mantle  velocity  model  inter¬ 
preted  by  Prod ehl  (1979)  from  several  short  (150— 
250  km)  refraction  profiles  in  the  vicinity  of 
Eureka,  NV  (Fig  12)  involving  shotpoints  at 
Eureka,  NV;  Delta,  UT;  Hiko,  NV;  Elko,  NV; 
and  NTS.  This  is  the  crustal  model  given  in  Pro- 
dehl's  (1979)  Table2  and  shown  here  in  Fig  14. 


SHOAL  DELTA  (EUREKA) 


p-vuocmr  Ikm/tl 

Fig.  14.  P-velodty  and  Q  structure  models  near  Eureka,  NV. 
Tbe  SD9/SD10/SDI6  Models  used  for  the  synthetic  iismo- 
grasn  calculations  differ  only  is  tbe  top  3  km.  Assumed  explo¬ 
sion  source  at  0  J  km  depth  for  all  calculations, 
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We  then  adjusted  the  details  of  the  P-velocaty 
model  only  for  the  uppermost  3  km  (Models  SD10, 
SD9  and  SD16,  shown  in  Fig.  14).  Our  purpose 
here  is  not  to  make  a  complete  reinterpretation  of 
the  SHOAL-Dclta  profile,  but  only  to  demon¬ 
strate  how  the  presence  or  absence  of  low-velocity 
surface  layers  determines  the  corresponding  pres¬ 
ence  or  absence  of  P-bar  phases  on  Basin  and 
Range  record  sections.  Our  synthetic  modeling 
indicates  that  we  could  continue  improving  the  fit 
to  the  relatively  sparse  SHOAL-Dclta  observa¬ 
tions  by  adjusting  deeper  crustal  velocity  discon¬ 
tinuities  and/or  gradients  as  well  as  the  near¬ 
surface  velocity  structure.  However,  we  feel  that 
these  other  model  adjustments  might  unduly  con¬ 
fuse  the  key  physical  issue  of  the  P-bar /low-veloc¬ 
ity  surface-layer  relation  that  we  want  to  em¬ 
phasize  here.  For  all  SD- series  Models  the.  source 
depth  was  held  constant  at  0.3  km  —  the  ap¬ 
proximate  SHOAL  explosion  depth.  Poisson’s  rate 
was  assumed  to  be  0.25  (Fs  =  Vr/  v/3 ),  and  a 
Nafe-Drake  velocity-density  relationship  was  as¬ 
sumed  (Olsen  and  Braile,  1981).  The  assumed  vari¬ 
ation  of  Q  with  depth  is  also  shown  in  Fig.  14.  The 
assumed  peak  frequency  for  the  source  spectrum 
with  U  Hz. 

The  crustal  model  that  gives  the  better  agree¬ 
ment  with  the  observed  record  section  is  SD9, 
which  has  a  pronounced  low-velocity  layer  (Vr  — 
4.0  km  s"\  3  km  thick)  at  the  surface.  The  ob¬ 
served  and  synthetic  record  sections  are  compared 
in  Fig.  13. 


4.  Discussion 

The  key  role  of  near-surface  layers  in  control¬ 
ling  the  amplitude  and  character  of  the  P-^ar 
phase  can  be  followed  by  comparing  Figs.  6,  13 
and  15.  Figure  6  is  a  comparison  of  the  reflection 
coefficients  (square  root  of  energy)  for  near-surface 
P  to  P  reflection  (Rorr) loT  near-critical  to  grazing 
angles  of  incidence  of  upward-traveling  P-waves. 
These  plane-wave  reflection  functions  were  calcu¬ 
lated  using  a  Thomson- Haskell  layer  matrix 
method  (Haskell,  1966),  which  is  the  same  algo¬ 
rithm  used  in  the  reflectivity  program  (Kind,  1978). 
Although  the  details  of  the  interference  fringes  are 


somewhat  different,  both  the  SD9  Model  (uniform 
3  km  layer  of  4  km  s-1  sediment)  and  the  SD16 
Model  (uniform  gradient  from  Vf  =  4.0  km  s~ 1  at 
0.6  km  depth  to  Vf  —  5.7  km  s' 1  at  3.0  km)  show 
appreciable  P  energy  returned  to  the  crust  for 
energy  in  the  —60s  incidence  angle  and  0.5- 1.5 
Hz  range,  which  are  the  ranges  in  which  most  of 
the  energy  is  transmitted  for  these  models.  Figures 
15(B)  and  (C)  show  that  the  integrated  result  of 
.  such  multiple  reflections  is  a  well-developed  P-bar 
phase  in  both  cases.  Conversely,  the  SD10  Model 
with  no  low-velocity  surface  layer  (Fig.  6(A))  has 
nearly  zero  P-P  conversion  and  completely 
supresses  P-bar  propagation  (Fig.  15(A)).  Since 
'there  are,  undoubtedly,  lateral  variations  in  the 
near-surface  layer  structure  and  velocity,  waves 
propagating  in  the  real  Earth  will  “average’’  minor 
variations,  and  a  detailed  matching  of  P-bar  wave¬ 
form  and  modulation  characteristics  is  not  justi¬ 
fied.  However,  the  total  energy  in  the  P-bar  phase 
as  represented  by  the  RMS  amplitudes  ovei  the 
several-second  duration  of  the  P-bar  phase  is  a 
quantity  that  depends  on  the  presence  of  low- 
velocity  surface  layers  over  the  propagation  path, 
but  this  energy  is  relatively  insensitive  to  assumed 
fine  details  in  the  synthetic  models!  so  long  as 
some  low-velocity  layers  are  present  (see  also 
Fig.  7(B)). 

In  summary,  the  short-period  regional  crustal 
phases,  commonly  known  as  P-bar  and  Lg,  propa¬ 
gate  as  multiple  P-wave  and  S-wave  reflections  in 
a  crustal  waveguide  whose  principal  boundaries 
are  the  Moho  and  the  free  surface,  as  shown 
schematically  in  Fig.  11.  Near-surface  low-velocity 
layers  create  constructive  interference  effects  for 
upcoming  P-wave  energy  at  wide  angles  of  inci¬ 
dence,  which  results  in  a  large  fraction  of  the 
P-wave  energy  being  retained  in  the  waveguide 
and  hence  in  efficient  propagation  of  the  P-bar 
phase.  For  no  surface  low-velocity  layers,  nearly 
all  upcoming  P  energy  converts  to  SV  and  is  lost 
to  further  P-bar  propagation.  Although  the  Moho 
and  the  free  surface  are  also  the  principal  wave¬ 
guide  boundaries  for  Lg.  Lg  propagation  effi¬ 
ciency  is  little  affected  by  near-surface  low  veloci¬ 
ties. 

Comparison  of  our  synthetic  modeling  with  the 
observations  discussed  by  Barket  et  al.  (1982)  sug- 
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gests  that  there  are  two  important  factors  is  the 
observed  complexity  of  P-bar  and  Lg  phases.  The 
first  is  the  complicated  interference  patterns  of 
propagation,  even  for  very  simple  laterally  homo¬ 
geneous  crustal  layers  such  as  we  have  discussed 
here,  and  the  second  is  the  local  geology  or  site-re¬ 
sponse  effects.  The  first-order  effect  is  that  of 
laterally  homogeneous  layering;  propagation 
within  this  structure  is  sufficient  to  produce  very 
complex  seismograms  which  display  little  lateral 
coherency,  as  we  have  shown.  Consideration  of  the 
effects  of  local  geological  conditions,  as  discussed 
by  Barker  et  al.  (1982),  leads  to  further  waveform 
and  amplitude  complexities  which  are  difficult  to 
disentangle  in  single-station  seismograms. 
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Abstract.  The  amplitude  of  the  Pg  phase,  as  recorded  in 
explosion  seismology  studies,  is  analyzed  with  the  aid  of 
synthetic  seismograms.  Parameters  such  as  source  fre¬ 
quency,  low-velocity  cover  above  the  crust  (sediments  or 
weathered  layer),  low-velocity  layers  within  the  upper  crust, 
velocity  gradients,  thickness  of  the  gradient  zone,  attenua¬ 
tion  and  Poisson's  ratio  strongly  influence  the  amplitude- 
distance  pattern  of  the  Pg  phase.  A  systematic  study  clearly 
shows  that  different  models  of  the  continental  upper  crust 
display  distinct  amplitude-distance  characteristics.  These 
models  could  not  be  distinguished  by  travel-time  interpreta¬ 
tion  alone. 

In  the  presence  of  gradient  zones  the  amplitude-distance 
curve  shows  different  patterns  depending  on  the  source  fre¬ 
quency.  The  higher  the  frequency,  the  more  pronounced 
are  the  relative  maxima  in  the  amplitudes.  The  presence 
of  a  low-velocity  cover  at  the  surface  accentuates  the  char¬ 
acter  of  the  amplitude-distance  curves  even  if  the  cover 
is  thin  (a  few  hundred  meters).  Moreover,  a  low-velocity 
cover  produces  P  to  S  conversions  and  multiples  following 
the  Pg  which  obscure  possible  secondary  crustal  phases. 
The  thickness  of  the  velocity  gradient  zone  influences  the 
amplitude  decay  and  the  width  of  the  relative  maxima. 
Low- velocity  layers  within  the  upper  crust  cause  a  faster 
drop-off  of  the  amplitudes  than  would  be  expected  from 
ray  theory.  Detailed  Pg  amplitude  studies  are  thus  useful 
in  improving  the  knowledge  of  the  physical  properties  of 
the  upper  continental  crust.  The  application  of  the  derived 
criteria  to  two  sets  of  real  data  allow  us  to  determine  fine 
details  of  the  velocity-depth  function  which  are  of  great 
importance  for  the  understanding  of  the  earth's  crust. 

Key  words:  Explosion  seismology  -  Upper  continental  crust 
-  Seismic  amplitude  -  Source  frequency  -  Low-velocity 
layer  -  Velocity  gradient. 


Introduction 

Understanding  the  fine  structure  and  physical  properties 
of  the  continental  crust  is  one  of  the  principle  goals  of 
explosion  seismology.  Interpretation  of  seismic  refraction 
data  using  traditional  travel-time  methods  gives  only  a 
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rough  picture  of  the  velocity-depth  structure.  The  increasing 
use  of  the  ray  tracing  interpretation  techniques  and  meth¬ 
ods  derived  from  the  Herglotz-Wiechert  travel-time  inver¬ 
sion  (including  r-p  methods)  result  in  velocity  models  of 
the  upper  crust  that  include  velocity  gradients  which  are 
very  often  poorly  defined.  In  fact,  layers  with  constant  ve¬ 
locity,  the  simplest  model,  fit  the  travel-time  data  equally 
well  in  most  cases.  Healy  (1963)  has  shown  that  very  differ¬ 
ent  velocity-depth  models,  from  homogeneous  layers  to 
continuous  gradient  zones  can  fit  equally  well  the  travel- 
time  data  for  a  particular  phase.  However,  as  will  be  shown 
below,  the  amplitude-distance  character  of  these  models 
may  be  significantly  different.  At  present,  either  very  de¬ 
tailed  travel-time  information  or  amplitude  studies  are  the 
only  techniques  available  for  accurate  determination  of 
velocity  gradients  in  the  earth’s  crust. 

In  the  following  discussion,  we  use  the  term  upper  crust 
for  the  upper  10-15  km  of  t'.e  crystalline  continental  base¬ 
ment  lying  immediately  beneath  the  surface  sediments  and 
above  the  lower  crustal  layer.  The  P- wave  velocity  of  the 
upper  crust  normally  varies  between  5.7  and  6.3  km/s. 
These  velocities  are  characteristic  of  sialic  rocks  at  the  ap¬ 
propriate  upper  crustal  temperature  and  pressures.  The 
compressional  seismic  wave  critically  refracted  (head  wave) 
in  the  upper  crust  is  usually  called  Pg.  In  explosion  seismic 
studies,  the  Pg  phase  is  normally  recognized  as  the  first 
arrival  in  the  distance  range  of  about  10-100  km.  The  nota¬ 
tion  Pg  is  also  used  in  earthquake  seismology  studies  but 
in  these  cases  it  often  refers  to  a  different  phase  at  larger 
distances. 

Studies  of  the  Pg  phase  recorded  in  refraction  profiles 
(Muller  and  Fuchs,  1976;  Muller  and  Mueller.  1979;  Banda 
and  Ansorge,  1980;  Braile  et  al.,  1982)  have  shown  how 
amplitude  information  can  be  used  to  greatly  reduce  the 
range  of  models  fitting  the  travel-time  interpretation.  As 
refraction  surveys  are  becoming  more  detailed,  with  closely 
spaced  recordings  and  improved  amplitude  control,  we  feel 
that  the  qualitative  comparisons  attempted  so  far  in  most 
amplitude  studies  are  not  enough,  and  that  understanding 
the  variation  of  Pg  amplitudes  as  a  function  of  various 
parameters  will  provide  further  insight  into  the  velocity 
structure  of  the  upper  crust.  In  turn,  this  information  will 
serve  as  a  basis  for  comparison  with  laboratory  measure¬ 
ments  and  petrological  studies  leading  to  a  better  knowl¬ 
edge  of  the  physical  properties  of  the  basement. 

In  this  paper  we  discuss,  on  the  basis  of  synthetic  se¬ 
ismograms,  the  Pg  amplitude-distance  curves  and  their  vari- 
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ation  with  parameters  such  as  frequency  content  of  the 
source,  presence  of  low-velocity  cover  and  of  low-velocity 
layers  within  the  upper  crust,  velocity  gradients,  thickness 
of  the  gradient  zone,  attenuation  (Q~  *)  and  Poisson’s  ratio. 
We  also  present  examples  of  amplitude-distance  modelling 
for  observed  Pg  amplitude  data  from  central  Europe  and 
western  North  America. 

Methods  of  Computation 

Synthetic  seismograms  can  at  present  be  calculated  by  a 
variety  of  methods.  For  our  purposes  we  have  used  the 
reflectivity  method  developed  by  Fuchs  (1968)  and  Fuchs 
and  Muller  (1971)  with  the  fundamental  modification  by 
Kind  (1978)  and  the  asymptotic  ray  method  described  by 
terveny  et  al.  (1977).  These  two  methods  and  others  have 
been  extensively  discussed  in  the  literature  (see  Chapman, 
1978  and  Spudich  and  Orcutt,  1980  for  reviews). 

Since  the  objective  of  this  paper  is  to  describe  in  detail 
the  amplitude-distance  behaviour  of  the  Pg  wave  and  to 
understand  related  later  arrivals  (multiply  reflected,  re¬ 
fracted  and/or  converted),  we  have  used  mainly  the  reflecti¬ 
vity  method.  The  modification  by  Kind  (1978)  includes 
these  secondary  effects  by  taking  into  account  the  free  sur¬ 
face  and  placing  the  source  in  the  reflectivity  zone.  This 
technique  requires  velocity  models  consisting  of  laterally 
homogeneous  and  isotropic  layers.  Velocity  gradients  are 
approximated  by  a  stack  of  thin  layers  with  small  velocity 
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contrasts.  For  our  computations  we  have  used  thicknesses 
corresponding  to  less  than  one  wavelength  of  the  frequency 
range  considered. 

Record  sections  of  the  vertical  component  of  ground 
displacement  have  been  calculated  for  all  the  models  with 
a  distance  interval  of  5  km  between  0  and  100  km  (Fig.  1). 
Except  where  stated,  the  VfjVs  ratio  was  assumed  to  be 
)/T  (Poisson  ratio  =  0.25)  and  the  Qr  values  were  fixed  at 
100  for  the  sediments  or  weathered  layer  and  at  500  for 
the  basement  with  Qs  equal  to  4Q,/9.  The  depth  of  the 
source  was  fixed  at  100  m  depth. 

All  record  sections  were  plotted  with  a  reduction  veloci¬ 
ty  of  6.0  km/s  and  the  amplitudes  of  each  trace  were  multi¬ 
plied  by  the  distance  for  a  more  accurate  reading  of  the 
amplitudes  (Fig.  1).  Amplitudes  were  read  taking  the  maxi¬ 
mum  (peak-to-peak)  of  the  first  cycle  (reading  the  first  pulse 
has  given  identical  results  for  theoretical  seismograms).  Fi¬ 
nally  the  amplitude  readings  were  plotted  as  function  of 
distance,  as  shown  in  Fig.  1. 

The  phase  velocity  interval  used  to  determine  the  range 
of  angles  of  incidence  over  which  the  reflectivity  program 
integrates  was  0. 2-0.4  km/s  below  the  minimum  velocity 
in  the  models  and  1. 0-2.0  km/s  above  the  maximum.  This 
range  includes  all  of  the  compressional  and  shear  waves 
of  interest  propagating  in  the  upper  crust.  The  large  phase 
velocity  range  leads,  especially  for  high  frequencies,  to  very 
long  computer  runs  (e.g.  for  the  model  in  Fig.  1,  23.3  min 
CPU  time  on  a  Cyber  M0722A).  However,  it  is  worthwhile 


Fig.  I.  Example  of  synthetic  record-section  calculated  using  the  reflectivity  method  with  a  dominant  frequency  of  4  Hz  for  model 
PG!  IK.  displayed  in  the  left  inset.  The  amplitude-distance  curve  isolid  line)  is  shown  together  with  the  curves  for  l/.v  and  Use*  decays 
in  the  right  inset.  Note  the  prominent  secondary  phases:  Pm.  whispering  gallery  (see  terveny  et  al.,  1977),  PgiS).  direct  Pg  converted 
to  S  at  the  base  of  the  sediments,  PgPg.  Pg  reflected  once  at  the  free  surface,  PgPg(S).  reflected  Pg  converted  to  5 
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Table  1.  Model  parameters 


Model 

Lower 
limit 
of  layer 
(km) 

Layer 

thick¬ 

ness 

(km) 

y, 

(km/s) 

Number  Gradient 
of 

layers 

(lun/s/km) 

PG1 

2.0 

2.0 

4.0 

2 

0.0 

15.0 

13.0 

6.0 

5 

0.0 

20.0 

5.0 

6.5 

2 

0.0 

PG3 

2.0 

2.0 

4.0 

2 

0.0 

15.0 

13.0 

5.9-6.1 

5 

0.015 

20.0 

5.0 

6.1 

2 

0.0 

PG6 

2.0 

2.0 

4.0 

2 

0.0 

15.0 

13.0 

5.75-6.25 

5 

0.038 

2t  )  • 

5.0 

6.25 

2 

0.0 

PGll 

2.0 

2.0 

4.0 

2 

0.0 

15.0 

13.0 

5. 5-6.5 

5 

0.077 

20.0 

5.0 

6.5 

2 

0.0 

PG11K 

2.0 

2.0 

4.0 

2 

0.0 

7.0 

5.0 

5.5-5.885 

5 

0.077 

20.0 

13.0 

5.885 

3 

0.0 

PGll  It  K. 

2.0 

2.0 

4.0 

2 

0.0 

10.0 

f  0 

5.5-6.12 

10 

0.077 

20.0 

10.0 

6.12 

2  ’ 

0.0 

PG1IL 

2.0 

2.0 

4.0 

2 

0.0 

7.0 

5.0 

5.5-5.885 

5 

0.077 

12.0 

5.0 

5.5 

1 

0.0 

20.0 

8.0 

6.3 

1 

0.0 

PG11LL 

2.0 

2.0 

4.0 

2 

0.0 

10.0 

8.0 

5.5-6.12 

8 

0.077 

15.0 

5.0 

5.5 

I 

0.0 

20.0 

5.0 

6.3 

1 

0.0 

PG12 

2.0 

2.0 

4.0 

2 

0.0 

15.0 

13.0 

6.1-5.9 

5 

-0.015 

20.0 

5.0 

6.5 

2 

0.0 

PG22 

2.0 

2.0 

4.0 

2 

0.0 

15.0 

13.0 

5.3-6.7 

5 

0.108 

20.0 

5.0 

6.7 

2 

0.0 

SULZ3 

J  0 

1.0 

4.5 

2 

0.0 

1.3 

0.3 

4.8-56 

1 

2.67 

5  0 

3.7 

5.6-S.9 

5 

0.081 

10.0 

5.0 

5.9-6.1 

7 

0.04 

21.0 

11.0 

6.1 

2 

0.0 

SULZ4 

1.0 

1.0 

4.3 

2 

0.0 

1.3 

0.3 

4.8-5.6 

1 

2.67 

6.3 

5.0 

5.6-6.0 

6 

0.08 

17.0 

10.7 

6.0 

2 

0.0 

SULZ6 

1.0 

1.0 

4.6 

2 

0.0 

5.7 

4.7 

5.6-5.95 

7 

0.074 

8.2 

2.5 

5.95-6.0 

4 

0.02 

17.0 

8.8 

6.0 

1 

0.0 

SULZ7 

1.0 

1.0 

4.6 

2 

0.0 

5.7 

4.7 

5.6-5.95 

7 

0.074 

7.0 

1.3 

5.95-5.98 

3 

0.023 

17.0 

10.0 

5.98 

1 

0.0 

performing  these  computations  to  reveal  the  influence  of 
the  sediments.  Many  tests  were  run  before  choosing  the 
best  illustrative  models  for  the  purpose  of  this  paper,  which 
are  listed  in  Table  1. 

To  compare  the  performance  of  the  reflectivity  method 
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Fig.  2.  Amplitude-distance  curves  (peak-peak)  of  Pn  phase  for 
model  EC4  (Cerveny,  1979)  computed  using  reflectivity  ( continuous 
line)  and  using  ray-method  (dotted  line) 


with  the  asymptotic  ray  theory  method,  we  started  by  deter¬ 
mining  the  amplitude-distance  curves  for  one  of  the  models 
published  by  Cerveny  (1979)  for  which  he  used  both  meth¬ 
ods  (Fig.  2).  Although  this  case  deals  with  a  gradient  zone 
in  the  upper  mantle,  the  situation  is  analogous  to  a  gradient 
zone  in  the  upper  crust  with  a  sedimentary  layer  above. 
The  ray  method  obviously  sharpens  the  peak  of  the  ampli¬ 
tude-distance  curve.  This  is  because  the  asymtotic  ray  meth¬ 
od  represents  a  “high  frequency  approximation”  to  the 
wave  equation  (Cerveny  et  al.,  1977;  Chapman,  1978).  As 
is  shown  below,  similar  results  were  obtained  from  other 
computations. in  this  study.  The  ray  method  is  inexpensive 
and  for  much  of  the  available  data  and  some  crustal  models, 
which  can  also  include  lateral  inhomogeneities,  this  approx¬ 
imation  is  accurate  enough.  For  more  detailed  studies,  in 
which  the  models  can  be  approximated  by  flat  homoge¬ 
neous  layers,  the  reflectivity  method  is  more  appropriate 
and  has  therefore  been  used  in  this  paper. 

The  synthetic  seismograms  computed  in  this  study  rep¬ 
resent  ground  displacements  instead  of  ground  velocities, 
as  measured  in  observed  seismograms.  As  shown  in  the 
example  in  Fig.  3,  for  which  both  velocity  and  displacement 
were  calculated,  the  difference  in  the  amplitude-distance 
behaviour  is  not  significant.  At  least  for  the  models  pre¬ 
sented  here,  the  results  from  displacement  can  thus  be 
applied  directly  to  observed  velocity  data.  Moreover,  it 
should  be  noted  that  the  dominant  frequency  used  in  dis¬ 
placement  computations  is  increased  when  the  displacement 
seismograms  are  differentiated  to  obtain  velocity  (see 
spectra  in  Fig.  3). 
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Fig.  3.  Comparison  of  amplitude-distance  curves  of  displacement  Fig.  4.  Amplitude-distance  curves  for  model  PGllK  [upper  nghi 
( continuous  line)  and  velocity  ( dotted  line )  seismograms  for  model  inset)  for  frequencies  of  2.5,  4.  5  and  8  Hz  (Founer  spectra  in 

PG 1  IKK.  (upper  right  inset).  The  corresponding  spectra  show  the  the  lower  inset)  computed  using  reflectivity  method  ( commous 

frequency  shift  between  displacement  and  velocity  signals.  The  lines).  Dotted  line  corresponds  to  ray-method  computation  with 

bottom  insets  show  the  signals  (at  35  km)  for  displacement  and  program  SH1S4 

velocity;  the  intervals  used  to  calculate  the  spectra  are  marked 
by  vertical  bars 


Discussion  of  Model  Parameters 

Source  Frequency 

The  frequency  content  of  seismic  refraction  data  varies  from 
about  2  Hz  to  more  than  20  Hz  depending  on  the  shooting 
technique,  charge  size,  frequency  response  of  the  instrument 
and  local  geological  environment.  Comparison  of  the  shape 
of  the  amplitude-distance  curves  for  low  frequency  source 
signals  ( —  2  Hz)  with  those  for  high  frequency  ( *  8  Hz) 
shows  significant  differences  when  velocity  gradients  are 
truncated  at  shallow  depth  or  low-velocity  layers  are  pres¬ 
ent  in  the  model.  For  that  reason  we  have  computed  most 
of  the  models  for  frequencies  2-8  Hz. 

Figure  4  shows  an  example  of  the  results  for  model 
PG1  IK  computed  for  2.5, 4,  6  and  8  Hz  dominant  frequen¬ 
cies.  The  fact  that  the  velocity  structure  of  the  upper  crust 
changes,  from  a  positive  gradient  (0.077  km/s/km)  to  zero 
gradient  at  7  km  depth  produces  a  different  response  of 
the  medium  depending  on  the  source  frequency.  The  varia¬ 
tion  of  amplitude  with  distance  is  more  pronounced  at  the 
higher  frequencies.  As  discussed  in  more  detail  below,  this 
is  due  to  the  fact  that  shorter  wavelengths  are  affected 
mainly  by  the  focusing  effect  of  the  gradient  zone,  thus 
producing  strong  relative  amplitude-maxima  at  a  distance 
of  45  km.  At  the  same  distance,  longer  wavelengths  are 
already  affected  by  the  homogeneous  layer  beneath  the  gra¬ 
dient  zone.  Higher  frequencies  are  thus  more  informative 
in  amplitude  studies. 


For  the  same  model,  amplitudes  were  computed  with 
the  asymptotic  ray  theory  and  are  also  displayed  in  Fig.  4. 
As  was  to  be  expected,  the  ray  method  works  reasonably 
well  for  high  frequencies,  although  the  slope  of  the  ampli¬ 
tude  decay  is  somewhat  enhanced.  Asymptotic  ray  theory 
accounts  for  the  influence  of  different  frequencies  on  the 
amplitudes  of  waves  reflected  at  first  order  discontinuities 
(terveny  et  al.,  1977)  but  not  for  waves  refracted  from  a 
gradient  zone  (Banda,  1979).  Therefore,  if  we  are  dealing 
with  good  quality  data,  suitable  for  amplitude  studies,  it 
is  of  fundamental  importance  to  compute  the  theoretical 
seismograms  using  the  reflectivity,  or  other  wave  theory 
method,  with  a  source  which  has  a  dominant  frequency 
similar  to  that  of  the  experimental  data. 

Thickness  and  Velocity  Structure  of  Sediments 

Qualitative  differences  in  the  character  of  the  wave  field 
due  to  the  presence  of  sediments  are  well  known.  Their 
influence  was  studied  quantitatively  with  models  which  are 
identical  except  for  the  velocity  structure  of  the  sediments. 

Figure  5  shows  the  results  for  models  with  0,  0.2,  2  and 
5  km  of  sediments  overlying  an  upper  crust  containing  a 
velocity  gradient  of  0.077  km  s/km.  Beyond  about  60  km 
the  shapes  of  the  amplitude-distance  curves  are  essentially 
identical.  The  most  important  difference  is  that  even  for 
models  with  as  little  as  0.2  km  of  sediments,  a  relative  maxi¬ 
mum  of  the  amplitude  curve  is  observed  which  is  missing 
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g.  5.  Amplitude-distance  curves  for  model  PG11  with  0.0,  0.2, 
3  and  5.0  km  of  low-velocity  material  ( inset  upper  right)  overlying 
e  basement 


Fig.  6.  Amplitude-distance  curves  for  model  PG1I  (2.5  Hz)  with 
different  velocity  structures  of  the  low-velocity  cover  of  the  base¬ 
ment  Unset  loner  left) 
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Fig.  7.  Record-sections  for  model  PG11  (4  Hz)  with  and  without 
low-velocity  cover  (0.2  and  0.0  km,  lower  and  upper  sections  re¬ 
spectively).  Note  the  difference  in  the  wave  field  of  secondary  arriv¬ 
als  (see  also  Fig.  1) 


for  the  model  without  sediments.  We  interpret  this  as  an 
effect  of  the  change  of  the  angle  of  incidence  affecting  the 
transmission  coefficient.  On  the  other  hand,  models  with 
a  low-velocity  cover  alone  and  without  a  gradient  in  the 
upper  crust  do  not  have  a  relative  amplitude  maximum. 

Results  for  models  with  different  velocity  structures  in 
the  sediments  are  displayed  in  Fig.  6,  which  show  that  there 
is  no  significant  difference  in  the  overall  shape  of  the  Pg 
amplitude-distance  curves.  However,  the  level  of  the  ampli¬ 
tudes  at  the  local  maximum  around  60  km  relative  to  the 
amplitudes  at  15  km  is  lower  for  sediments  with  a  velocity 
of  5  km/s  than  for  those  of  4  km/s  or  for  sediments  with 
a  velocity  gradient  starting  with  4  km/s  at  the  surface.  This 
again  is  due  to  the  fact  that  the  angle  of  incidence  is  stee¬ 
pened  by  the  lower  velocity  at  the  suiface,  thus  producing 
higher  P-wave  amplitudes  on  the  vertical  component 
seismograms  studied  here. 

As  a  result,  we  can  state  that  any  influence  of  the  thick¬ 
ness  and  velocity  structure  of  the  sedimentary  layer  on  the 
amplitude-distance  behaviour  of  the  first  cycle  of  the  Pg 
phase  is  restricted  to  shorter  distances.  However,  the  total 
wave  fields  for  models  with  and  without  sediments  are  radi¬ 
cally  different.  Reverberations  within  the  sediments  and  P 
to  S  conversions  result  in  conspicuous  seismic  phases  that 
appear  after  the  Pg  phase  (Fig.  7,  see  also  Fig.  1). 


Thickness  of  Gradient  Zone  and  Low-Velocity  Layers 

Significant  differences  in  the  amplitude-distance  curves  for 
models  having  the  same  gradient  but  different  thickness 
of  the  gradient  zone  are  shown  in  Fig.  8.  The  models  include 
structures  with  a  continuous  gradient  between  2  and  1 5  km 
(model  PGll),  2  and  10  km  (model  PGUKK)  and  2  and 
7  km  (model  PG11K)  on  top  of  a  half  space.  A  decrease 
in  the  thickness  of  the  gradient  zone  leads  to  a  faster  drop¬ 
off  of  the  amplitudes  with  distance. 

The  introduction  of  a  low-velocity  layer  below  a  gra¬ 
dient  zone  between  2  and  7  km  and  2  and  10  km  (models 
PG11L  and  PG11LL)  show  another  interesting  effect.  A 
significant  shift  of  the  maximum  and  a  change  in  the  slope 
of  the  amplitude  decay  is  evident  when  a  low-velocity  layer 
is  present  at  the  same  depth  at  which  the  gradient  is  termin¬ 
ated  (compare  PG1  IK  and  PG1 1 L  in  Fig.  8). 

The  influence  of  the  frequency  content  of  the  source 
in  the  presence  of  a  low-velocity  layer  is  shown  in  Fig,  9. 
The  main  effects  are  to  sharpen  the  maximum  and  to  shift 
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Fig.  8.  Amplitude-distance  curves  for  models  PGll,  PG11K, 
PG 1 1 KK  with  variable  gradient  zone  thickness  ( continuous ,  dashed 
and  doited  lines,  respectively  in  the  lower  inset)  and  PG11L  and 
PG11LL  with  low-velocity  layer  at  different  depths  ( dashed  and 
continuous  lines,  respectively,  in  the  upper  inset )  computed  for  a 
frequency  of  4  Hz 


it  to  larger  distances  as  the  frequency  is  increased.  Model 
PG11L  peaks  very  smoothly  for  2.5  Hz  at  about  30  km 
whereas  it  peaks  sharply  at  40  km  for  the  higher  frequen¬ 
cies.  The  waves  behave  as  if  they  “sense"  U.e  low-velocity 
layer  well  before  this  would  be  expected  from  ray  considera¬ 
tions.  To  illustrate  the  effect  of  the  truncated  gradients, 
the  ray  tracing  for  model  PG1  IK  is  shown  in  Fig.  10.  Here, 
even  though  the  gradient  zone  is  terminated  at  a  depth 
of  7  km,  we  observe  rays  emerging  at  distances  out  to 
58  km,  which  is  in  disagreement  with  the  results  from  these 
amplitude  calculations.  As  stated  above,  a  decrease  in  thick¬ 
ness  of  the  gradient  zone  causes  the  Pg-amplitude  drop-off 
to  occur  at  smaller  distances  than  expected.  For  models 
PG11K  and  PGUL,  in  which  the  gradient  zone  is  termin¬ 
ated  at  7  km  depth,  the  amplitudes  start  deviating  from 
those  of  model  PGll,  whose  gradient  extends  down  to 
15  km,  at  distance  as  short  as  45  km.  As  is  shown  in  Fig.  10 
rays  arriving  at  this  distance  reach  a  maximum  depth  of 
only  4  9  km,  which  is  well  above  the  bottom  of  the  gradient 
zone.  In  fact,  the  distance  between  the  turning  point  of 
the  ray  and  the  bottom  of  the  gradient  zone  is  greater  than 
one  wavelength  (about  1.5  km  at  4  Hz).  This  phenomenon, 
which  can  not  be  explained  by  geometrical  ray  theory,  is 
analogous  to  the  Fresnel-zone  effect  in  the  case  of  electro¬ 
magnetic  waves:  the  energy  arriving  at  a  receiver  is  due 
not  only  to  the  ray  propagating  with  minimal  travel-time, 
but.  as  a  consequence  of  Huyghens'  principle,  consists  of 
difractions  interfering  with  each  other,  whose  travel-times 
are  greater  than  that  of  the  direct  wave  (Bom,  1933).  Modi- 


Flg.  9.  Amplitude-distance  Curves  for  models  PGl  1 LL  and  PG1 1 L 
(continuous  and  dotted  lines  respectively,  in  the  inset )  computed 
for  2.S  and  4  Hz 
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Fig.  10.  Ray-trace  for  model  PGl  IK,  with  the  ray  emerging  at 
45  km  and  the  ray  grazing  the  lower  boundary  of  the  gradient 
zone  emerging  at  58  km  (solid  curves).  The  lower  dashed  curve 
shows  the  edge  of  the  first  Fresnel-zone  (corresponding  to  a  phase 
shift  of  T/2)  while  the  upper  dashed  curve  corresponds  to  a  phase 
shift  of  m  for  a  signal  with  4  Hz  dominant  frequency 


fications  of  asymptotic  ray  theory  which  take  this  effect 
into  account  have  been  introduced  by  various  authors  (see 
Spudich  and  Orcutt,  1980  for  a  review)  In  particular. 
Wiggins  (1976)  approximates  the  wave  behaviour  by  en¬ 
visaging  the  energy  transported  by  a  ray  as  being  distributed 
over  a  disk  travelling  with  that  ray,  and  constructs  a  synthet¬ 
ic  seismogram  by  summing  the  contribution  of  neighbour¬ 
ing  rays  as  their  disks  intersect  the  surface  at  the  distance 
of  interest.  As  can  be  seen  from  the  dashed  ray-paths  in 
Fig.  10,  the  first  Fresnel  zone,  which  contains  the  most  sig¬ 
nificant  contributions  to  the  amplitude  of  the  signal, 
extends  to  a  depth  beyond  the  limit  of  the  gradient  zone 
for  energy  arriving  at  a  distance  of  45  km.  Thus,  changes 
in  gradient  will  influence  the  amplitudes  of  waves  with  lon¬ 
ger  periods  at  smaller  distances  than  those  with  shorter 
periods.  This  is  aiso  the  reason  for  the  marked  shift  of 
the  amplitude  peak  in  models  PGl  1 L  and  PGl  1  LL  between 
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Fig.  II.  Amplitude-distance  curves  for  models  PG12  (negative  gra¬ 
dient).  PG1  (no  gradient)  and  PG3,  PG6,  PG11  and  PG22  (positive 
gradients)  for  a  frequency  of  2.5  Hz  (see  Table  1  for  model  specifi¬ 
cations) 


2.5  and  4  Hz,  and  for  the  frequency  dependence  ,  *n 
in  Fig.  4. 

Reflections  originating  from  the  low-velocity  layer  have 
also  been  studied.  Our  results  do  not  differ  from  those  pub¬ 
lished  by  Braile  and  Smith  (1975),  Smith  et  al.  (1975)  and 
further  modified  by  Banda  (1979):  the  reflection  from  the 
top  of  the  low-velocity  layer  will  only  be  seen  as  a  separate 
phase  when  the  upper  crust  has  a  very  small  or  zero  gra¬ 
dient,  the  frequency  is  high  enough  and  the  transition  is 
nearly  a  first  order  discontinuity.  In  this  case  the  interfer¬ 
ence  of  Pg  with  the  reflection  changes  the  Pg  amplitudes 
at  larger  distances  (this  has  not  been  studied  in  detail  in 
this  paper).  If  the  above  mentioned  conditions  are  not  ful¬ 
filled  the  Pg  amplitudes  will  barely  be  affected  by  the  reflec¬ 
tion  from  the  top  of  the  low-velocity  layer. 


Variation  of  Gradients 

It  is  well  known  that  even  a  slight  positive  velocity  gradient 
greatly  influences  the  amplitude  of  refracted  phases 
(Cerveny,  1966;  Hill,  1971).  Thus,  we  have  tried  to  vary 
systematically  the  velocity  gradient  in  the  basement  models 
in  order  to  determine  the  corresponding  amplitude-distance 
curves.  We  have  arbitrarily  chosen  the  gradient  to  be  linear. 
Amplitudes  for  models  with  various  gradients  between 
-0.015  and  0.108  km/s/km  were  computed  and  are  shown 
in  Fig.  11.  From  a  comparison  with  a  calculation  for  model 
GP11  at  higher  frequency  and  from  the  discussion  in  the 
previous  section,  it  can  be  stated  that,  although  these  curves 
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Fig.  12.  Amplitude-distance  curves  for  models  PGl,  PG6  and  PG1 1 
(see  inset  and  Table  1)  for  Q  =  500,  continuous  lines,  Q  =■=  200, 
dotted  lines  and  Q  =  2000  (dashed  line  for  model  PG  1 1 ) 


were  calculated  for  2.5  Hz,  they  are  also  representative  of 
higher  frequencies  as  long  as  the  gradient  zones  extend  to 
sufficient  depth.  The  results  shown  in  Fig  11  illustrate  that 
small  positive  velocity  gradients  in  the  upper  crust  will  be 
resolvable  by  amplitude  measurements  on  reasonably  good 
experimental  data.  For  example,  the  amplilude-d  •— 'oe 
characteristics  of  Models  PGl  and  PG6  are  significantly 
different  although  the  velocity  structure  differs  only  by  the 
presence  of  a  small  (0.038  km/s/km)  gradient  in  Model 
PG6. 

Attenuation  (Q) 

We  have  tested  models  with  Q  values  for  compressional 
waves  of  200,  500  and  2000  for  the  upper  crust  to  investigate 
the  influence  of  attenuation  upon  the  amplitude  of  the  Pg 
wave.  The  curves  shown  in  Fig.  12  are  not  sufficiently  dis¬ 
tinct  for  us  to  infer  apparent  Q  from  Pg  amplitude  calcula¬ 
tions.  As  already  demonstrated  by  Hill  (1971)  a  slight  chan¬ 
ge  in  the  gradient  could  make  up  for  the  differences  shown 
in  Fig.  12  without  taking  the  attenuation  into  account.  To 
estimate  Q  in  the  basement  one  could  use  the  method  pro¬ 
posed  by  Braile  (1977)  which  requires  that  the  Pg  ampli¬ 
tudes  be  modelled  simultaneously  with  the  amplitudes  of 
the  reflection  from  the  bottom  of  the  upper  crust.  Alterna¬ 
tively,  if  the  geometrical  spreading  factors  are  sufficiently 
well  known,  and  the  band-width  of  the  source  is  broad 
enough,  it  is  conceivable  that  spectral  ratio  methods  (see 
B3th,  1974  for  a  review)  might  be  adequate  to  infer  Q  in 
the  upper  crust. 
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Fig.  13.  Bottom:  ray-trace  for  model  Sulz-7.  Note  vertical  exagerauon  of  model  plot.  Middle  record  section  of  Sclz-south  data,  vertical 
component,  bandpass  filtered  8-16  Hz,  trace  normalized.  Travel-time  curve  corresponds  to  ray-trace  below  Top.  synthetic  record  section 
of  flat  layer  approximation  of  model  Sulz-7  (see  Fig.  14  and  Table  1)  vertical  component  velocity  amplitude  multiplied  by  distance 


Poisson's  Ratio  (a) 

Seismic  velocity  studies  of  the  crust  show  that  a  Poisson 
ratio  of  0.25  is  in  general  a  good  average.  However,  in 
some  cases  strong  deviations  from  this  value  have  been 
found.  The  influence  of  a  has  been  studied  for  a  few  models 
with  values  0.2-0.35  The  resulting  amplitude-distance 
curves  do  not  differ  significantly  from  each  other.  Only 
curves  for  models  involving  a  sudden  change  of  a  at  a 
discontinuity  produce  distinct  features  in  the  Pg  versus  re¬ 
flection  amplitude  ratio  as  already  shown  by  Olsen  et  al. 
(1979). 

Examples  of  Data  from  the  Black  Forest  (Germany) 
and  the  Basin  and  Range  (U.S.A.) 

Between  1974  and  1980  several  quarry  blasts  near  Sulz, 
southern  Germany,  were  recorded  along  a  1 1 3  km  long 
profile  running  south  along  the  eastern  margin  of  the  Black 
Forest  into  the  Swiss  Alpine  foreland.  The  first  2  s  of  the 
resulting  seismic  record-section  are  shown  in  Fig.  13  togeth¬ 
er  with  the  travel-time  curve  corresponding  to  the  ray-trace 
model  in  the  lower  part  of  the  figure.  The  ray  tracing  was 
performed  using  the  method  of  Gebrande  (1976). 

The  velocity  structure  of  the  sediments  (including  a 
wedge  of  late  Paleozoic  and  early  Mesozoic  sediments  with 
P-velocity  of  3.8  km  s  between  0  and  40  km)  and  the  shape 
of  the  basement/sediments  boundary  were  modelled  to  fit 
the  available  bore-hole  data  (Buechi  et  al.,  1965;  Boigk  and 


Schoeneich,  1968;  Lemcke  et  al.,  1968).  The  velocity  struc¬ 
ture  of  the  basement  was  then  adjusted  until  the  general 
shape  of  the  travel-time  data  was  matched.  No  attempt 
was  made  to  model  the  local  variations  in  sedimentary 
structure  and  basement  depth  which  produce  the  small  time 
discrepancies,  limited  to  one  or  two  consecutive  records. 
This  applies  particularly  to  the  local  anomaly  between  94 
and  101  km  which,  because  of  the  ooor  signal  to  noise  ratio, 
cannot  be  resolved  with  these  data  alone. 

For  plotting  convenience,  the  seismograms  in  Fig.  13 
were  trace  normalized.  The  peak-to-peak  amplitudes  of  the 
first  cycle  were  multiplied  by  the  scale  factor  marked  above 
each  seismogram  to  obtain  the  true  particle  velocity  values 
in  nm/s.  The  amplitude  data  for  the  Pg  phase  are  plotted 
as  crosses  with  corresponding  shot  and  station  numbers 
in  Fig.  14. 

The  size  of  shot  number  1  was  only  700  kg  while  the 
others  were  around  2,000  kg,  so  that  the  corresponding  am¬ 
plitude  values  were  multiplied  by  a  correction  factor  of  1  4, 
equal  to  the  cube  root  of  the  charge  ratio,  which  seems 
to  be  appropriate  for  this  specific  quarry.  The  lemaimng 
scatter  seems  to  be  independent  of  the  shots,  and  although 
it  amounts  to  about  a  factor  of  two.  the  data  define  an 
amplitude-distance  behavior  characterized  by  a  rapid  de¬ 
crease  in  the  first  20  km,  a  local  maximum  at  about  40  km 
and  a  smooth  decay  out  to  about  100  km. 

In  order  to  be  able  to  apply  the  reflectivity  method  to 
the  data,  the  curved  layers  derived  from  the  rav-tracing 
technique  were  approximated  by  fiat  layers  in  model 
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Fig.  14.  Amplitude  date  of  Sulz-South  record  section  (Fig.  13)  (crosses  with  shot  and  station  numbers)  and  amplitude-distance  curves 
for  the  four  models  shown  in  the  lower  inset  and  listed  in  Table  1.  Top  inset  shows  an  enlarged  view  of  the  gradient  zones  of  the 
four  models 


number  7  (insets  in  Fig.  14  and  Table  1)  and  the  sediments 
were  simplified  to  a  single  layer.  The  resulting  synthetic 
record  section  (vertical  ground  velocity  with  a  dominant 
frequency  of  8  Hz)  is  reproduced  at  the  top  of  Fig.  13.  The 
offset  of  the  calculated  amplitude-distance  curve  as  a  whole 
was  adjusted  to  the  experimental  data  by  a  least-squares 
method  (Fig.  14,  curve  7).  For  three  other  ray-trace  models 
which  satisfy  the  travel-time  data  equally  well,  amplitude- 
distance  curves  were  calculated  in  the  same  way  and  plotted 
in  Fig.  13  together  with  their  velocity-depth  functions.  This 
illustrates  some  of  the  possible  model  variations.  A  compar¬ 
ison  between  models  3  and  4  shows  how  a  second  gradient 
zone  of  sufficient  strength  and  extent  can  significantly  in¬ 
crease  the  distance  at  which  the  amplitude  decay  occurs. 
Models  6  and  7  illustrate  the  effect  of  small  changes  in 
the  extent  of  a  second  zone  with  slight  gradient.  The  loga¬ 
rithmic  standard  deviations  of  data  points  from  these  curves 
range  between  0.23  and  0.27,  which  correspond  to  ampli¬ 
tude  factors  of  1.7  and  1.9  respectively.  The  amplitude- 
distance  behaviour  of  the  models  presented  here  do  not 
differ  from  each  other  sufficiently  to  be  able  to  discriminate 
between  them  on  the  basis  of  this  data  alone.  However, 
as  the  insets  in  Fig.  14  show,  the  velocity-depth  variations 
of  these  models  are  very  small,  so  that  the  results  are  good 
evidence  for  a  strong  gradient  (0.07-0.08  km/s/km)  in'  the 
upper  5-6  km  of  the  basement. 

It  should  be  noted  that  computational  techniques  re¬ 
strict  us  to  consideration  of  fiat-layered  models  for  the  am¬ 
plitude  modelling  even  though  small  lateral  variations  are 
evident  in  the  travel-time  model  for  the  Sulz  data  as  shown 


in  Fig.  13.  Shooting  up-dip  over  the  first  30  km  and  down- 
dip  beyond  that  distance  will  have  a  slight  focusing  effect 
on  the  rays,  thus  contributing  somewhat  to  the  observed, 
relative  amplitude  maximum  at  that  distance.  However,  be¬ 
cause  the  lateral  changes  in  velocity  structure  indicated  in 
Fig.  13  are  small  (note  that  the  model  is  plotted  with  2  x 
vertical  exageration)  and  because  the  positive  velocity  gra¬ 
dients  in  the  fiat-layered  models  which  satisfy  the  amplitude 
data  also  fit  the  travel-time  data  modelled  with  a  curved 
sediment-basement  boundary,  we  anticipate  that  the  ampli¬ 
tude-distance  effects  of  the  two-dimensional  structure  will 
be  negligible. 

A  very  different  amplitude-distance  character  of  the  Pg 
phase  is  observed  for  the  Basin  and  Range  province  of 
the  tectonically  active  western  North  American  continent. 
A  Pg  data-set  from  this  region  was  analyzed  for  this  study. 
The  seismic  records  are  from  the  northern  Basin  and  Range 
(NBR)  (shotpoints  Mountain  City,  Eureka  and  Elko)  in 
Western  United  States  and  were  originally  studied  by  Hill 
and  Pakiser  (1966)  and  presented  by  Prodehl  (1970, 1979). 

Partial  record  sections  emphasizing  the  Pg  arrivals  for 
the  NBR  data  are  shown  in  Fig.  15.  The  record  sections 
are  from  Prodehl  (1979)  and  the  travel-time  curves  shown 
are  the  results  of  ray-trace  modelling  by  Fauria  (1981).  Am¬ 
plitudes  of  the  Pg  phase  were  read  from  the  sections  shown 
in  Fig.  IS  and' corrected  for  plot  scaling  factors  using  the 
calibration  signals  presented  on  the  original  sections 
(Prodehl,  1979).  The  Pg  amplitudes  were  adjusted  to 
account  for  the  differing  shotpoints  and  plotted  as  a  func¬ 
tion  of  distance  in  Fig.  16.  Although  there  is  considerable 
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Fig.  15.  Partial  record  lections  of  the  Pg  phase 
for  the  northern  Basin  and  Ranee  province. 
Travel-time  curves  shown  are  the  result  of  ray- 
trace  modelling  by  Fauna  (1981) 


scatter  in  these  data  (the  range  is  illustrated  by  the  shaded 
area  of  the  plot)  it  appears  that  Pg  amplitudes  for  the  NBR 
decay  rapidly  with  distance  and  do  not  display  -he  local 
maximum  which  was  characteristic  of  the  Sulz  data  and 
which  indicated  a  velocity  gradient.  In  fact,  the  NBR  Pg 
amplitude  data  are  fitted  bat  by  the  PCI  model  consisting 
of  homogeneous  upper  crust  The  scatter  in  the  data  shown 
in  Fig.  16  and  the  fact  that  the  data  represent  a  compilation 
from  several  different  refraction  lina  preclude  more  de¬ 
tailed  analysis.  Note  that  the  very  small  positive  or  negative 
gradients  (leu  than  that  for  model  PG3)  and  virtually  any 
Q  value  could  be  present  in  the  upper  crust  of  the  NBR 
based  on  these  data  alone.  However,  because  no  strong 
gradient  or  low-velocity  effects  are  observed  in  the  ampli¬ 
tude  distance  data,  the  upper  crust  in  the  NBR  appears 
to  be  nearly  a  homogeneous  velocity  layer. 


The  fine  structure  of  the  continental  basement,  i.e.  the  dis¬ 
tribution  of  physical  parameters  with  depth,  can  reveal  sig¬ 
nificant  aspects  of  the  crustal  evolution  and  the  interaction 
of  different  crustal  units  as  well  as  its  thermal  and  composi¬ 
tional  history.  Very  often  the  upper  part  of  the  crust  does 
not  show  pronounced  interfaces  with  discontinuities  in  ve¬ 
locity  and  density.  Vertical  incidence  reflection  cannot 
resolve  velocity  structure  in  the  absence  of  distinct  imped¬ 
ance  boundaries.  On  the  other  hand,  combined  travel-time 
and  amplitude'  analysis  of  refraction  data  can  give  more 
detailed  insight  into  the  velocity  structure.  Therefore  com¬ 
bined  reflection  and  refraction  surveys  should  be  carried 
out  in  crustal  investigations. 

A  representative  velocity-depth  structure  requires  either 
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Fig.  16.  Amplidide'distanoe  data  for  the  Pg  phase  for  the  NBR 
data  shown  in  Fig.  13.  The  shaded  region  represents  the  range  of 
scatter  of  the  data  points.  Theoretical  amplitude-distance  curvet 
calculated  by  the  reflectivity  method  are  shown  for  models  PG1, 
PG3  and  PG12  (Table  1) 


detailed  and  highly  accurate  travel-time  data,  i.e.  a  spacing 
of  observation  points  on  the  order  of  a  few  hundred  meters, 
or  a  combination  of  travel- time  and  reliable  amplitude  data 
obtained  from  an  average  record  spacing  of  a  few  kilome¬ 
ters.  The  following  conclusions  drawn  from  this  study  can 
help  to  achieve  a  better  interpretation  of  the  amplitude  date 
and  hence  lead  to  a  more  reliable  estimate  of  velocity-depth 
structure  of  the  upper  continental  crust 

1.  The  frequency  content  of  the  source  is  a  critical  parame¬ 
ter  when  gradient  zones  and/or  discontinuities  of  velocity 
are  present  Higher  frequencies  produce  more  pronounced 
variations  of  amplitudes  with  distance.  The  amplitude  vari¬ 
ations  are  likely  to  be  resolvable  within  the  scatter  of  the 
data  and  thus  high  frequency  source  signals  lead  to  better 
resolution. 

2.  A  low-velocity  cover  of  the  basement  accentuates  relative 
maxima  in  the  amplitude-distance  curve  when  velocity  gra¬ 
dients  are  present  in  the  crystalline  upper  crust  Layer  thick¬ 
nesses  of  only  0.2  km,  i.e.  even  weathered  layers,  can  cause 
this  effect  In  addition,  multiple  reflections,  refractions  and 
conversions  in  the  low-velocity  cover  contribute  considera¬ 
bly  to  a  complicated  wave  field  following  the  initial  Pg 
arrival  and  probably  make  later  arrivals  from  deeper  parts 
of  the  crust  undetectable. 

3.  Reflections  from  the  top  of  a  low-velocity  layer' within 
the  upper  cruet  will  affect  the  Pg  amplitudes  by  interference 
only  if  they  are  strong  enough.  This  requires  low  gradients 
in  the  basement  above,  a  first  order  discontinuity  at  the 
bottom  of  the  upper  crust  and  a  high  frequency  source 
signal. 


4.  The  increasing  thickness  of  the  gradient  zone  in  the  base¬ 
ment  reduces  the  amplitude  decay  with  distance  and  broad¬ 
ens  the  width  of  the  relative  amplitude  maximum 

5.  The  Pg  amplitudes  decay  with  distance  much  faster  in 
the  presence  of  a  low-velocity  layer  than  would 'be  expected 
from  geometrical  ray  theory. 

6.  Reliable  Pg  amplitude  data,  characterized  by  a  continu¬ 
ous  decrease,  or  a  steady  value  over  some  distance  or  even 
a  local  maximum,  allows  us  to  distinguish  between  homoge¬ 
neous  constant  velocity  layers  and  various  gradient  zones 
in  the  upper  crust 

7.  Apparent  Q  values  cannot  be  extracted  from  the  ampli¬ 
tude  distance  behaviour  of  Pg  data  alone. 

8.  Synthetic  seismogram  calculations  based  on  asymptotic 
ray  theory  can  be  used  only  as  a  first  approximation  to 
determine  detailed  velocity  depth  structures. 

These  results,  of  course,  apply  only  to  a  crust  which 
can  be  approximated  by  latently  homogeneous  models.  Fo¬ 
cusing  effects  due  to  lateral  heterogeneities  may  completely 
mask  the  amplitude-distance  behaviour  due  to  vertical 
structures  and  must  be  studied  by  other  methods. 
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ABSTRACT 


Esplndola,  Juan  Manuel.  Ph.D.,  Purdue  University,  August  1979.  Finite 
Difference  Synthetic  Seismograms  for  Kinematic  Models  of  the 
Earthquake  Source.  Major  Professor:  Lawrence  W.  Bralle. 

The  dynamic  displacement  field  of  kinematic  dislocation  In  two 
dimensions  Is  modeled  by  a  finite-difference  technique  In  skew 
coordinates.  The  scheme  Involves  the  solution  of  the  heterogeneous, 
elastic  vector  wave  equation  subjected  to  dislocation  conditions. 
Algorithms  are  Included  to  model  the  free  surface  and  absorbing 
artificial  boundaries  which  simulate  an  infinite  half-space. 

Specification  of  faults  with  arbitrary  dip  Is  Implemented  by 
means  of  the  skew-coordinate  system.  Other  source  parameters  such  as 
time-history,  slip-function,  and  rupture  velocity  are  user  selected. 

Near-field  synthetic  seismograms  were  calculated  by  this  pro¬ 
cedure  for  selected  source  parameters  for  two  cases  of  geologic  Interest 
(1 .e. ,  burled  basin  and  layer  over  a  half  space).  These  examples 
Illustrate  the  Influence  of  fault  parameters  and  local  geology  on 
ground  motion.  Significant  effects  due  to  the  heterogeneous  media  and 
various  fault  and  source  parameters  suggest  that  the  method  can  provide 
a  new  realistic  modelling  tool  for  the  study  of  the  Interaction  between 
the  earthquake  source  and  the  local  geology  In  the  near-field  for  long 
period  ground  motion. 


ABSTRACT 


Maxxelia.  Frederick  Elwood.  Ph.D. .Purdue  Univaraity,  August 
1979.  Tha  Ganaratlon  of  Synthetic  Seismograms  for  Laterally 
Heterogeneous  Modal a  Using  tha  Pinita  Dlffaranea  Taehniqua. 
Major  Profaaaor:  Lawrence  W.  Braila. 

A  aathod  for  bha  ganarabion  of  aynbhatic  seismograms  for 
two-dimensional  habaroganaoua  mods  Is  specified  only  by  P- 
wav a  (dilatation)  valoeiby  haa  baan  davalopsd.  Tha  finiba 
dlffaranea  approximation  taehniqua  is  applied -to  tha  two- 
dimanaional  vaetor  wava  aquation  for  an  alastie  and 
iaobropie  aodal.  Tha  flnlta  dlffaranea  calculations  ara 
appliad  bo  a  das Ignat ad  portion  of  tha  aodal  at  any 
particular  tine  step  during  tha  lifetime  of  tha 
approximation.  This  portion  of  the  aodal  is  referred  to  as 
tha  window  and  is  prescribed  to  include  tha  waves  of 
interest  which  propagate  through  tha  modal.  A  suite  of 
models  ara  evaluated  with  various  window  paraaabars  which 
provides  insight  into  bha  limitations  and  capabilities  of 
tha  method.  Tha  computational  affieiancy  of  tha  windowing 
taehniqua  allows  evaluation  of  models  which  ara 
approximately  twice  as  large  as  those  which  can  be  evaluated 
by  other  numerical  approximation  techniques  not  employing  a 


windowing  scheme  for  tho  mi*  computing  rosourcoo  (time  and 
memory) . 

Two  typoo  of  boundary  conditions  aro  employed  on  the 
walls  of  the  window.  The  upper  horizontal  surface  is 
treated  with  a  second  order  heterogeneous  free  surface 
boundary  eondition  developed  during  this  study.  All  other 
surfaces  of  the  window  are  treated  with  a  nonreflecting 
boundary  eondition  to  minimize  false  reflections  and  mode 
conversions.  Seismograms  can  be  obtained  for  any  location 
inside  or  on  the  boundaries  of  the  sodel. 

Synthetic  seismograms  are  calculated  for  seven 
heterogeneous  (both  vertical  and  horizontal)  geologically 
realistic  models.  Each  model  is  described  graphically  and  • 
is  accompanied  by  a  displacement  seismogram  record  section 
for  the  vertical  and  horizontal  components.  These  record 
sections  reveal  significant  amplitude  and  time  effects  which 
result  from  the  velocity  structure. 


